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FOREWORD 


This  report  was  prepared  under  USAF  Contract  No.  F33615-70-C-1126 , 
"Unsteady  Aerodynamics  and  Flutter  of  Wing  Horizbntal  Tail  Configurations 
in  Supersonic  Flow",  as  a  part  of  Project  No.  1370,  "Dynamic  Problems  in 
Military  Flight  Vehicles"  and  Task  No.  137003,  "Prevention  of  Dynamic 
Aeroelastic  Instabilities  in  Advanced  Military  Aircraft".  The  work  was 
performed  by  the  Flutter  Research  Group  of  The  Boeing  Company,  Commercial 
Airplane  Group,  Renton,  Washington,  for  the  Aerospace  Dynamics  Branch, 
Vehicle  Dynamics  Division,  Air  Force  Flight  Dynamics  Laboratory,  Wright- 
Patterson  Air  Force  Base,  Ohio.  The  computer  program  was  prepared  and 
documented  by  Mr.  Gordon  D.  Kramer  and  Mr.  George  F.  Keylon  of  Boeing 
Computer  Services,  Inc.  The  overall  program  management  was  provided  by 
Mr.  William  S.  Rowe  of  The  Boeing  Company.  Mr.  Lawrence  J.  Huttsell  of 
the  Aerospace  Dynamics  Branch  (FYS)  monitored  the  contract  for  the 
Air  Force. 

This  document  constitutes  the  first  part  of  a  two-part  final  report. 

Part  II  containing  the  computer  program  listing  is  available  upon  request 
from  the  Air  Force  Flight  Dynamics  Laboratory  (FY) ,  Wright-Patterson  AFB, 
Ohio  45433. 

This  report  covers  work  performed  from  December  1969  to  April  1971. 

The.  manuscript  was  released  by  the  authors-  Dr.  Jack  Morito^Er^  Mr.  Christo¬ 
pher  J.  Borland  and  Mr.  John  R.  Hog.ley  of  the  Flutter  Re^arciv  Group  of 
Commercial  Airplane  Group  for  publication  as  an  AFFDL  technical  report. 

This  report  is  also  known  as  Boeing  Document  D6-24860. 

This  technical  report  has  been  reviewed  and  is  approved. 


WALTER  J.  MYKYTOW 
Assistant  for  Research  &  Technology 
Vehicle  Dynamics  Division 
Air  Force  Flight  Dynamics  Laboratory 


ABSTRACT 


A  three-dimensional  extension  of  the  Mach  Box  technique  has  been  de¬ 
veloped  for  the  unsteady  aerodynamic  analysis  of  non-planar  wings  and  wing- 
tail  configurations  in  supersonic  flow.  Various  refinement  procedures  have 
been  included  to  improve  the  accuracy  of  the  results, 

A  general  purpose  computer  program  has  been  written  for  the  CDC  6600. 
The  program  is  capable  of  treating  ving  tail  combinations  with  or  without 
vertical  separation,  longitudinal  separation  and  dihedral  on  either  surface 
A  non-intersecting  wing-vertical  tail  combination  may  be  examined.  If  a 
wing  alone  is  treated,  perturbation  velocity  components  in  the  flow  field 
may  be  found. 

The  program  has  been  used  to  obtain  pressure  distributions  and  general 
ized  forces  on  several  wing-tail  configurations.  In  addition,  correlations 
with  experimental  flutter  results  have  been  performed  for  several  models 
tested  in  the  low  supersonic  Mach  number  range  at  Cornell  Aeronautical 
Laboratory. 
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^  n  i  b  i  n  «  x 

=Li 

TL1FT 

non- 

dimensional 

Nondimensional  amplitude  of 
total  lift  (  program  output  ) 

—  |  bi  \  *  t  ut'i 

ii  ■  ij  •  -ufh  e 

t 

EL 

non- 

dimensional 

The  l  coordinate  location  of  a  pulse 
c 

receiving  point,  i.e.,  the 
perpendicular  distance  from  the 
sending  plane  to  receiving  point. 

M 

XMACH 

non- 

dimensional 

Mach  number 

KiS 

— 

mass 

Generalized  mass  for  the  mode 

m 

— 

non- 

dimensional 

The  m  coordinate  location  of  a 
c 

pulse  receiving  point. 

m. 

i 

— 

mass 

i*k  lumped  mass 

n,m 

N 

xyz 

ENRUS 

ENRLS, 

EN, 

ENSUBD 

non- 

dimensional 

Normal  wash  at  box  n,m  on  surface 
"xyz"  due  to  local  source  strength, 
where  possible  subscript  values  are; 

x  •[ R-right 
\L-left 

y  «*jU-upper 
(L-lower 

z  »/W-wing 
VT-tall 

e.g.  means  normal  wash  on  the 

right  upper  wing  at  box  n,m 

xlll 


NOMENCLATURE 


Mathe¬ 

FORTRAN 

Dimension 

Definition 

matical 

Symbol 

Symbol 

A  n»m 

N 

ENRULU, 

non- 

Normal  wash  at  box  n,m  on  surface 

xyz 

ENRLLL, 

dimensional 

"xyz"  due  to  remote  source  strengths 

abc 

ENRURW, 

ENRULW 

on  surface  "abc",  where  possible 
subscript  values  for  a,b,  and  c 
are  the  same  as  for  x,y,  and  z, 
respectively,  defined  above; 
a  n»m 

e.g.  N  means  normal 

RUT 

LLW 

wash  at  box  n,m  on  the  right  upper 
tail  due  to  source  strengths  on  the 
left  lower  wing. 

N 

s 

NSUBDV 

non- 

dimensional, 
odd  integer 

No.  of  "sub-boxes"  (chordwise 
and  spanwise)  to  be  used  in  the 
subdivision  improvement  technique. 

n 

— 

non- 

dimensional 

the  n  coordinate  location  of  a 
c 

pulse  receiving  point. 

nc,V  ^c 

— 

— 

Sending  Surface  Coordinate  System 

l*° 

1  B° 

1  c° 

— 

— 

Receiving  Point  Coordinate 

System 

p»p(x,y,t) 

— 

force/area 

local  static  pressure 

P« 

— 

force/area 

Free  stream  static  pressure 

A  p(x,y,0 

force/area 

pressure  difference  between  upper 
and  lower  surfaces  at  point  (x,y) 
at  time  t 

^  p(x,y,t)«*  P  (x,y,  t  )  “p(x,y,i). 

K  »•"  >  r  upper  H  iowei 

NOMENCLATURE 


Mathe¬ 

matical 

Symbol 

FORTRAN 

Symbol 

Dimension 

Definition 

^/>(x,y) 

— 

force /are a 

Amplitude  of  pressure  difference: 

e 

4j>(x,y,t)  =  Apix, 

force /length 

Generalized  for^g  due  to  the  defor¬ 
mation  in  the  l  elastic  mode  and 
loading  for  the  modal  deflections 

— 

force /length 

Amplitude  of  generalized  force 

% 

GENAF 

non- 

dimensional 

Non-dimenBional  generalized  force 
(program  output); 

~  iwt  twf 

QiJ  “  QiJS 

QAGARD 

1/ (length) 2 

Generalized  force  in  the  AGARD 
notation 

/  if 

V  Qu 

GENAFC 

1/ (length) 2 

Real  and  imaginary  parts  of  ,  in  the 
AGARD  definition  (program  output) 

q. 

— 

force /are a 

dynamic  pressure 

length 

Generalized  coordinate  relating  physical 
deflection  to  modal  deflections: 

Z(x,y,t)  =2  fjU;y)  ^  j(t) 

U 

5j 

— 

length 

11. 

Amplitude  of  J  generalized  coordinate 

8 

S 

length 

Wing  semi-span. 

NOMENCLATURE 


Mathe¬ 

matical 

Symbol 

FORTRAN 

Symbol 

Dimension 

Definition 

S 

— 

length^ 

Area  of  integration.  Bounded  by 
edge  of  plan form  plus  diaphragm 
and  lying  inside  the  forward  Mach 
cone  of  the  receiving  point. 

t 

— 

time 

Time 

U 

— 

length/ 

time 

Free  stream  velocity. 

V 

— 

non- 

dimensional 

Velocity  spatial  AIC  for  perturbation 
velocity  parallel  to  the  free  stream. 

U 

— 

length/ 

time 

Perturbation  velocity  in  the  stream- 
wise  direction,  positive  downstream. 

V;rx 

V 

non- 

dimensional 

Velocity  spatial  aerodynamic  influence 
coefficient  (AIC)  for  velocity  normal 
to  the  free  stream  and  parallel  to 
the  sending  surface. 

(xy} 

V 

y/X 

V 

non- 

dimensional 

Spatial  AIC  giving  velocity  normal  to 
the  free  stream  and  parallel  to  sur¬ 
face  "ab"  at  a  point  on  surface  "xy" 
or  in  the  flowfield  (FF),  due  to  con¬ 
stant  normal  wash  over  a  box  on 
surface  "ab" 

V 

length/ 

time 

Perturbation  velocity  in  the  span- 
wise  direction,  positive  right 
(looking  upstream). 

w 

non- 

dimensional 

Velocity  spatial  aerodynamic  influence 
coefficient  (AIC)  for  velocity  normal, 
to  the  sending  plane. 

NOMENCLATURE 


Mathe-  FORTRAN 

matical  Symbol  Dimension  Definition 

Symbol 


W  non-  Spatial  AIC  giving  velocity  normal 

dimensional  to  surface  "ab"  at  a  point  on  sur¬ 
face  "xy"  or  in  the  flow  field  (FF) 
due  to  a  constant  normal  wash  over 
a  box  on  surface  "ab". 


w  length/  Perturbation  velocity  in  the  vertical 

time  direction,  positive  upward. 


X,Y,Z  length  Reference  (globed)  coordinate  system* 

X  positive  aft,  Y  positive  right, 

Z  positive  upward. 


xvil 


NOMENCLATURE 


Mathe-  FORTRAN 

matical  Symbol  Dimension  Definition 

Symbol 


X 

w, 


Y  ,Z 
w 


w 


Wing  local  coordinate  system. 
Tail  local  coordinate  system 


^LE 

XWLE 

length 

XTLE 

XTLE 

length 

‘Vie 

XWTE 

length 

Y 

TTE 

XTTE 

length 

(?) 

f  XEDGE  Y 
\  XCENTR j 

length 

/  xcw] 

l  XCT, 


The  location  of  a  leading  edge 
definition  point  of  the  wing 
planform,  measured  along  the  X 
co-ordinate.  w 


Same  as  above  for  the  tail 
planform,  measured  along  the 
X^  coordinate. 

The  location  of  a  trailing  edge 
definition  point  of  the  wing 
planform. 


Same  as  above  for  the  tail, 
measured  along  the  X,j,  co-ordinate. 

Location  of  the  /edge  "{ 

\ center/ 

of  a  Mach  box  used  for  the 
placement  of  the  box  pattern, 
measured  along  the  X^  coordinate. 


The  location  of  the  most  upstream 
row  of  boxes  on  thei'wing) 

1  tail/ 

measured  along  thefXw-)  co-ordinate. 

w 


NOMENCLATURE 


Mathe¬ 

matical 

Symbol 

FORTRAN 

Symbol 

Dimension 

Definition 

ywle1  | 
ytle  / 

'  YWLE  } 
k  YTLE  / 

length 

The  location  of  a  leading  edge 
definition  point  of  thefwing! 

{tail/ 

planform  measured  along 
coordinate  ' 

Iyt;  * 

y 

YBAR 

non- 

dimensional 

Offset  of  receiving  chord  from 
the  nearest  sending  chord. 

zu' 

)  — 

length 

1  Upper| 

1  Lower J8urface  deflection  at 

(x,y)  as  a  function  of  time 

Zm(x,y,t) 

— 

length 

mean  surface  deflection: 

Zj^Xfy.t)  -  f j (x,y)elW^ 

Zx  (x,y) 

— 

length 

Local  thickness  at  (x,y) 

T*  (x,y) 

TSLFN 

non- 

dimensional 

Local  thickness  slope  at  (x,y) 

~r 

— 

non- 

dimensional 

Thickness  slope  piston  theory 
correction; 

S- 

*  2 

n,m 

X 

ALPHA 

non- 

dimensional 

Edge  box  area  ratio  for  box(n,m) . 

/3 

BETA 

non- 

dimensional 

}J  M2  -  1 
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nomenclature 


Mathe-  FORTRAN 


matical 

Symbol 

Symbol 

Dimension 

Definition 

\/ft 

BIBETA 

length 

Spanwise  dimension  of  the 
Mach  box. 

bl  //3 

s 

BIBTAS 

length 

Spanwise  dimension  of  a  sub¬ 
divided  Mach  box. 

Y 

GAMMA 

non- 

dimensional 

■  Ratio  of  specific  heats,  -  l. 

s 

ZETA 

non- 

dimensional 

Dummy  variable  in  the  Z  or  Z 

w 

coordinate. 

i 

ETABAR 

non- 

dimensional 

Dummy  variable  of  integration 
in  the  mc  coordinate. 

0 

THETA 

radians 

sin  -  sin-1  * 

0iJ 

radians/ 

length 

Torsion  of  ithlumped  mass  in 
mode  j. 

A 

— 

non- 

dimensional 

L  c  coordinate  location  of  a 
pulse  sending  box. 

M 

MU, YMU 

non- 

dimensional 

mc  coordinate  location  of  a 
pulse  sending  box. 

? 

YMUBAR 

non- 

dimensional 

™c  coordinate  location  of  a 

pulse  sending  box. 

XX 
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NOMENCLATURE 


Mathe- 

FORTRAN 

matical 

Symbol 

Symbol 

Dimension 

V 

NU,XNU 

non- 

dimensional 

XNUBAR 

non- 

dimensional 

\ 

XI 

non- 

dimensional 

? 

XIBAR 

non-dim 

dimensional 

f,  r. 

Mass/ 

volume 

X 

TAU 

non- 

dimensional 

0(x,y,t) 

— 

length2/ 

time 

Definition 


nc  coordinate  location  of  a 
pulse  sending  box. 


n£  coordinate  location  of  a 
pulse  sending  box. 


Dummy  variable  of  integration 
in  the  X  or  X  coordinate. 

W  i 

Dummy  variable  of  integration 

in  the  n  coordinate, 
c 

Free  stream  density. 


Disturbance  velocity  potential 
at  point  (x,y)  and  time  t, 
defined  so  that B 0  is  velocity, 

positive  in  positive  x^  direction, 
where 

xi  =  X,  Y,  or  Z 


4  0(x,y,t)  —  length2/ 

time 


Disturbance  velocity  potential 
difference  between  the  top  and 
bottom  side  of  the  surface  at 
point  (x,y)  and  time  t: 


<30(x,y,t) 


0 

upper 


lower 


NOMENCLATURE 


Mathe¬ 

matical 

Symbol 


FORTRAN 

Symbol 


Dimens ion 


Definition 


A  0j (x,y) 


length  /  Amplitude  of  velocity  potential 

time  difference  at  point  (?c,y)  or 

for  box  n,m  for  the  j*-**  mode 


d1j(x,y)]  DELPHI 


A  0 j (XT£ 


"I 


non-  Non-dimensional  velocity  potential 

dimensional  difference  due  to  the  unit 

generalized  coordinate  (program  output) ; 

Irt  n>m  U  /  n'm  ih,t 

**\  6  ' 

non-  Trailing  edge  velocity  potential 

dimensional  difference. 


degrees 


Dihedral  angle  of^ail)  *  radians, 
positive  upwards  from  horizontal. 


radians/ 

time 


Circular  frequency 


radians/ 

time 


Circular  frequency  of  mode  j 


Superscripts 


(n,m)  T 

(F,/0  } 


Box  location 


Subscripts 


Lower  limit  of  Integration;  Left-hand  surface;  Lower  surface 

Right-hand  Surface 

Subdivided 


Upper  limit  of  Integration  Upper  Surface 
Wing 

Flowfield 
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GLOSSARY  OF  TERMINOLOGY 


Aftmost  Box  -  Each  chord  on  each  planform  and  diaphragm  combination  has 
one  such  box.  It  is  the  aftmost  box  on  that  chord  for  which  AIC  arrays 
must  be  calculated  and  may  be  on  the  planform  or  diaphragm. 


AIC  -  Aerodynamic  Influence  Coefficient 


Area  Ratio  -  On-planform  fraction  of  a  box  which  is  cut  by  the  planform 
boundary. 


Apex  Box  -  The  box  on  the  sending  surface  which  encloses  the  apex 
of  the  Mach  hyperbola  associated  with  the  receiving  box. 


Box  Grid  -  Non-dimensionalized  geometric  array  of  boxes  whose  extent  is 
determined  by  the  geometric  properties  of  the  planforms.  The  term  "grid" 
embraces  the  arrays  on  both  surfaces. 


Control  Point  -  The  location  at  which  a  receiving  box  is  deemed  to  be 
influenced  by  other  boxes.  In  general,  the  center  of  the  receiving  box. 


Effective  Area  -  A  concept  which  relates  entirely  to  the  sub-division 
technique.  It  is  composed  of  those  boxes  sufficiently  close  to  the 
receiving  box  that  their  influence  on  it  is  large  enough  for  the  subdivision 
refinement  to  affect  results  significantly.  The  size  has  been  arbitrarily 
set  to  include  the  Nfi0X/Ng  rows  immediately  ahead  of  the  receiving  box. 


Global  Co-ordinate  System  -  An  overall  reference  system  of  co-ordinates. 
For  example,  the  airplane  co-ordinate  system  X  ~  aft,  Y^right,  Z^up. 

Y  -  0  at  centerline  of  airplane. 


Leading  Edge  Diaphragm  -  All  diaphragms  on  whichA0  ■  0. 


Local  Co-ordinate  System  -  A  co-ordinate  system  lying  in  the  plane  of 
the  surface,  x^aft,  y^root  to  right  tip.  y  ■  0  at  center  line  of 
airplane. 


Longitudinal  Separation  -  Streamwise  distance  between  the  trailing  edge 
of  the  wing  and  the  leading  edge  of  the  tail,  measured  along  the  centerline. 


GLOSSARY  OF  TERMINOLOGY 


Mach  Asymptote  -  The  asymptote  of  the  Mach  hyperbola. 


Mach  Hyperbola  -  The  intersection  of  the  sending  plane  and  the 
forward  Mach  cone  of  the  receiving  point.  Since  this  is  always  non- 
dimensionalized,  it  is  a  rectangular  hyperbola. 


Map  -  A  condensed  description  of  a  large  amount  of  data  which  can  be 
used  to  locate  any  desired  data  element.  A  map  of  a  banded  sparse  matrix 
might  consist  of  two  numbers  per  row,  the  first  being  the  first  non-zero 
column  of  that  row  and  the  second  being  the  band  width  for  that  row.  The 
matrix  itself  could  then  be  stored  as  band  elements  only. 


Normal  Offset  -  The  lc  distance  between  the  sending  box  and  the  receiving 
point. 

Parallel  Offset  -  The  m  distance  between  the  sending  box  center  and  the 
receiving  point. 


Partial  Box  -  A  sending  box  which  is  cut  by  the  Mach  hyperbola  but  which 
is  neither  an  apex  box  nor  an  edge  box. 


Planar  A.I.C.-  An  A.I.C  defined  by  the  geometric  relation  between  a 
sending  box  and  receiving  box  which  lies  in  the  same  plane.  Cr—  n 
only. 


Receiving  Box  -  In  defining  the  relationship  between  two  boxes  the 
receiving  box  is  the  box  which  can  be  influenced  by  the  other  box. 


Receiving  Chord  -  Those  receiving  boxes  which  lie  on  the  same  chord. 

The  receiving  chord  is  significant  in  that  all  the  boxes  lying  on  it  use 
AIC  arrays  which  are  a  subset  of  those  for  the  aftmost  box  lying  on  that 
chord. 


Sending  Box  -  In  defining  the  relationship  between  two  boxes,  the  sending 
box  is  the  box  which  influences  the  other  box  (c.f.  Receiving  Box). 


Spatial  A.I.C.  -  An  AIC  defined  by  the  geometrical  relationship  between 


two  boxes  which  do  not  lie  in  the  same  plane. 


V__T 

V/SA 


W_  _■ 
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GLOSSARY  OF  TERMINOLOGY 


Sub -box  -  A  member  of  the  array  of  boxes  formed  when  the  grid  of  sending  boxes 
is  subdivided.  Note  it  refers  to  the  small  box  which  is  a  fraction  of  the 
large  box,  and  not  to  a  large  box  which  has  been  subdivided. 


Tail  -  The  downstream  surface. 


Vertical  Separation  -  The  vertical  distance  between  the  center  lines  of 
the  two  surfaces.  Positive  if  the  second  surface  is  above  the  first. 


Wake  Diaphragm  -  That  vart  of  the  diaphragm  where  /JP  =  0  due  to 
the  influence  of  a  surface. 


Wing  -  Upstream  Surface  -  (E.g.  a  Canard  could  be  referred  to  as  a 
wing) 
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Figure  1  Coordinate  Systems  For  A  Right  Wing 


Symbol 

Transformation 

befinition 

Dimension 

X 

Global  or  Reference  Coordinate  System.  X  positive  Aft, 

Y 

Y  positive  Right,  Z  positive  Upward 

Length 

Z 

X 

1  w 

X  -  WLAX 

Wing  Local  Coordinate  System, 

Y 

w 

Ycos  ♦  +  (Z  -  WLAZ)sin 

used  to  define  wing  leading 

and  trailing  edges.  X  ,Y  ,Z 

Length 

Z 

w 

(Z  -  WLAZ)cos  ♦  -  Ysin  ♦ 

are  similarly  defined  for  the 

tail  local  axes 

n 

c 

<VV/bl  + 1 

Sending  Surface  Coordinate 

m 

Y  /(b./|  )  +  1/2 

System  used  to  define  box  grid. 

The  (n  ,m  )  plane  lies  within 

Non-dimen- 

Jtc 

1— _ 

W*> 

c  c 

the  plane  of  the  sending  surface, 
in  this  case  the  right  wing. 

sional 

|C  |  6  |H 


Symbol 


Transformation 

Definition 

Dimension 

Receiving  Point  Coordinate 

-(m  -m ) 

System  parallel  to  the 

-Uc-D 

n  ,m  ,1  coordinates  but 

opposite  in  sign  and  having 

non-dimen 

their  origin  at  the  pulse 

receiving  point  (n,m,l  in 

the  n  ,m  ,1  coordinates), 
c  c  c 

sional 

t 

WTz 


m 


* 


SYMBOL 

TRANSFORMATION 

DEFINITION 

DIMENSION 

*T 

yt 

X-TLAX 

Tail  Local  Coordinate 

Ycob  \j>T+(Z-TLAZ)sin^T 

System  used  to  define 

length 

ZT 

(Z-TLAZ)cos  r|>T-Ysin4)T 

tail  leading  and 

trailing  edges. 

n 

c 

{  Xt+1'LaX-(WLAX+Xc)  (  /b1+l 

Sending  Surface  Co- 

in 

V(V®  }  +  1/2 

ordinate  System.  In 

non- 

c 

this  case  the  right 

y  (y  0 ) 

tail  is  shown  as 

the  sending  surface. 

dimensional 
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SECTION  I 


INTRODUCTION 


This  report  describes  the  development  of  a  method  for  predicting  the  unsteady- 
aerodynamics  of  flexible  aircraft  with  non-planar  wing  and  tail  surfaces. 

The  aerodynamic  interference  between  the  wing  and  tail  has  been  taken  into 
account.  The  method  is  based  on  a  three-dimension  extension  and  refinements 
of  the  well-known  "Mach  Box"  method,  first  developed  in  Ref.  (l).  The  refine¬ 
ments  consist  of  a  subdivision  technique  for  improving  the  accuracy  of  the 
pressure  distributions,  a  velocity  potential  smoothing  technique,  and  a  thick¬ 
ness-correction  procedure  based  on  second-order  piston  theory.  The  theo¬ 
retical  developments  are  described  in  Section  II. 

A  general  purpose  computer  program  for  implementing  the  theory  has  been 
written  for  the  Control  Data  Corporation  6600.  All  information  necessary 
for  operation  of  the  program  is  contained  in  Section  III  of  the  report. 
Detailed  information  on  the  computer  program,  including  flow  charts  and 
listings,  is  contained  in  Part  II. 

The  program  has  been  tested  by  obtaining  results  on  s.everal  typical  con¬ 
figurations.  Unsteady  aerodynamic  results,  in  the  forms  of  pressure  distri¬ 
butions,  generalized  forces,  and  perturbation  velocities  in  the  flow  field, 
have  been  obtained  for  a  standard  AGARD  wing-tail  configuration  for  several 
Mach  numbers,  reduced  frequencies,  and  geometrical  arrangements.  The  pro¬ 
gram  has  also  been  used  for  comparison  of  experimental  and  analytical 
flutter  results  of  several  wing-horizontal  tail  configurations  tested 
in  the  low  supersonic  Mach  number  range  at  Cornell  Aeronautical  Laboratory 
(Ref.  (2)). 


The  possible  configurations  which  may  be  analysed  by  the  use  of  the  program 
described  herein  may  be  summarized  as  follows: 

a)  Single  Planar  Wings  (flow  field  sampling  permitted); 

b)  Single  Wings  with  Dihedral  up  to  ^45°  (flow  field  sampling 
permitted) ; 

c)  Coplanar  wing  and  tail; 

d)  Planar  wing  and  tail  with  vertical  separation; 

e)  Non  planar  wing  and  tail,  dihedral  on  either  surface  up  to  ^45°, 
with  or  without  vertical  separation; 

f)  Wing,  single  vertical  tail,  non-intersecting 

The  planforms  of  either  surface  are  entirely  general,  being  specified  by 
linear  boundaries  between  points,  for  up  to  ten  (x,y)  coordinates  on  either 
leading  or  trailing  edge.  Segments  of  the  leading  edge  may  not  be  swept 
forward.  Either  subsonic  or  supersonic  leading  or  trailing  edges  are  per¬ 
mitted.  Tips  may  be  streamwise  or  raked.  The  surfaces  are  not  permitted 
to  have  ’’mutual  interaction”  i.e.,  where  the  tail  directly  affects  the  wing 
as  well  as  vice-versa,  except  in  the  case  where  the  dihedral  angles  of  the 
two  surfaces  are  equal  and  there  is  no  vertical  separation.  This  includes 
the  coplanar  case. 

Structural  modal  or  rigid  body  deflections  may  be  input  at  arbitrary  points 
on  the  lifting  surfaces,  a$  polynomial  coefficients,  or  at  the  control 
points  of  a  previously  determined  Mach  box  grid.  The  program  is  capable  of 
recycling  on  Mach  number  and  reduced  frequency  without  reinput  of  mode  shapes 
or  geometry.  New  mode  shapes  may  be  input  without  reinput  of  geometry 
information.  The  reverse  may  also  be  accomplished  if  it  is  desired,  for 
example,  to  change  the  fineness  of  the  Mach  box  grid. 

The  applicable  Mach  number  range  of  the  program  is  greater  than  1.0  to  less 
than  5.0.  However,  because  of  the  limitations  of  linearized  theory,  results 
obtained  for  Mach  numbers  of  less  than  1.2  or  greater  than  3.0  should  be 
used  with  extreme  caution. 


SECTION  II 


THEORETICAL  DEVELOPMENT 


1.  A  BRIEF  SUMMARY  OF  THE  THEORY 

The  motion  of  the  configuration  under  consideration  Is  assumed  to 
be  simple  harmonic  motion  of  sufficiently  small  amplitude  that  the 
perturbation  velocity  potential  $(x,y,z,t)  satisfies  the  linearized 
differential  equation  of  invicid  fluid  motion  (Ref.  3): 
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wnere  U  is  the  free  stream  velocity  and  is  parallel  to  the  x-direction 
and  M  is  the  free  stream  Mach  Number. 

In  order  to  determine  the  pressure  distributions  on  the  planform, 
and  from  thence  the  generalized  forces,  the  method  of  Aerodynamic 
Influence  Coefficients  is  used.  As  is  described  in  References  3  and 
4  the  disturbances  of  the  flow  are  reproduced  by  distributing  pul¬ 
sating  sources  related  to  the  amplitudes  of  the  modal  displacements 
over  each  surface  of  each  planform  and  over  artificial  "diaphragm 
regions"  where  necessary  to  isolate  the  upper  and  lower  surfaces  of 
the  planform.  There  results  an  integral  expression  for  $(x,y,z) 
which  is  dependent  upon  the  distribution  of  sources  and  the  location 
of  the  point  (x,y,z).  This  equation  is  developed  in  sub-section 
2  below. 

This  integral  expression  cannot,  in  general,  be  evaluated  analytically. 
The  numerical  approach  used  divides  the  surface  into  a  grid  of 
rectangular  Mach  Boxes  whose  diagonals  lie  along  Mach  Rays,  and 
which  are  sufficiently  small  that  the  source  distribution  may  be 
considered  constant  within  them.  The  expression  for  <f>  then  becomes 
a  summation  over  the  Mach  Boxes.  Each  element  of  the  summation  is 
the  source  strength  at  the  box  center  multiplied  by  an  Aerodynamic 
Coefficient  which  is  dependent  only  on  the  relative  locations  of 
the  box  and  the  point  (x,y,z).  This  expression  is  derived  in  sub¬ 
sections  3  and  4. 

The  velocity  potentials  and,  where  unknown,  the  source  distributions 
are  then  evaluated  at  each  surface  of  each  planform  and  diaphragm, 
taking  into  account  the  boundary  conditions  at  the  surfaces  (section 
5).  The  pressure  distributions  and  generalized  aerodynamic  forces 
are  then  calculated  from  the  difference  between  the  velocity  potential 
distributions  on  the  upper  and  lower  surfaces  of  the  planform,  as 
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described  in  sub-section  6. 


It  is  well  known  that  the  Mach  Box  method  produces  pressure  distri¬ 
butions  which  contain  inaccuracies  in  the  fora  of  irregular  fluctu¬ 
ations.  Two  refining  techniques  to  reduce  the  fluctuations  are 
used  in  this  computer  program.  Of  prime  importance  is  the  sub¬ 
division  refinement  in  which  the  sending  Mach  Box  grid  is  divided 
into  an  even  finer  grid  while  the  receiving  grid  of  control  points 
remains  unchanged.  This  refinement  reduces  the  magnitude  of  the 
pressure  distribution  fluctuations  substantially.  A  second  refine¬ 
ment  is  the  application  of  a  least  squares  surface  or  chordwise  polynomial, 
fitted  to  the  velocity  potential  distribution,  which  smoothes  the 
pressure  distributions. 

The  pressure  distributions  are  further  refined  by  applying  to  them 
a  thickness  correction  based  on  second  order  piston  theory.  All 
these  optional  refinements  are  described  in  sub-section  7. 

THIN  AIRFOIL  THEORY  IN  SUPERSONIC  FLOW 

Thin  airfoil  or  small  perturbation  theory  is  applied  to  describe 
the  flow  patterns  that  result  when  3mall  disturbances  are  superposed 
on  parallel  uniform  flow.  When  coordinate  axes  are  fixed  on  a  thin 
body  placed  in  uniform  flow,  the  perturbation  velocity  potential 
satisfies  the  differential  equation. 

a  -  m2)  *„„  +  *  +  ♦„  -r1  *xt +  -V  *tt  (1) 
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The  linearized  equation  (1)  can  be  transformed  into  the  canonical 
wave  equation  through  the  mixed  Lorentz-Galilean  transformations, 
which  does  not  involve  any  spacial  rotation  but  is  concerned  only 
with  a  uniformly  moving  system,  and  which  preserves  the  speed  of 
sound  in  both  coordinate  systems  (Refs.  4  and  5), 

Then  equation  (1)  will  be  transformed  into 

*x'x'  +  Vy'  +  "  “V  *t't'  "  0  (2) 
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One  form  of  the  general  solution  of  equation  (2)  is 

♦Cx'.yV.f)  »(f  -JL!)  +  £(t.  +  i), 
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which  signifies  a  pulsing  source  singularity  moving  with  a  speed  of 
U  (Ref.  6),  and  where  F(x',y,,z')  is  independent  of  time  and  is 
related  only  to  the  spacial  coordinates,  and 
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Consider  a  lifting  surface  whose  leading  edge  is  supersonic  (for  conve¬ 
nience  the  y  -axis  is  placed  on  the  leading  edge) ,  so  that  there  is  no 
communication  between  the  upper  and  lower  surfaces  (see  Figure  3) . 

The  potential  at  P(x,y,z)  and  at  time  t  is  computed  from  a  source 
singularity  distribution  over  the  x,y  plane  and  is 

x  -  8z  y  +  n 
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The  non-time  dependent  function  F  is  determined  by  letting  z-»*-0+ 
and  ?*-0+.  The  result  is 

FU,n,0+)  - - ~~  W  (C,n,0+)  (8] 

2ne 

where  Wq  is  the  amplitude  of  the  upwash  w. 

If  the  time  variation  f(t)  of  the  source  singularity  is  harmonic, 
that  is, 

r  /  \  iU)t 

f(t)  -  e 

then,  by  applying  the  identities  (6)  the  following  can  be  obtained: 
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Therefore,  substitution  of  Equations  (8)  and  (9)  into  (7)  yields 
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where 


R^=  Hyperbolic  radius  =  B^r  =  [(x-O^-  B^((y-n)^  +  z^}]  ^ 
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and  0+  indicates  the  upper  surface  of  the  pulse-sending  box. 

From  the  previous  discussion,  Equation  (10)  is  only  valid  for  harmon 
ically  oscillating  source  singularities  in  which  the  wing  leading 
edges  are  supersonic. 

The  time-independent  part  of  Equation  (10)  may  be  written  a  little 
more  generally  as 


For  a  planform  with  a  supersonic  leading  edge,  S  indicates  that 
portion  of  the  planform  which  lies  inside  the  forward  Mach  cone 
emanating  from  (x,y,z).  If  the  planform  has  a  subsonic  leading 
edge,  or  trailing  edge,  or  a  streamwise  edge  at  its  tip,  "diaphragms 
are  introduced  using  the  "Eward  concept"  as  described  in  Ref.  7. 

The  diaphragms  are  artificial  "surfaces"  which  present  no  barrier 
to  the  flow,  across  which  there  can  be  no  pressure  discontinuity, 
and  through  which  continuity  of  mass  is  maintained.  They  are 
introduced  so  that  sources  may  be  distributed  upon  them  to  account 
for  the  effects  on  the  upper  surface  of  the  planform  induced  by 
the  source  distribution  on  the  lower  surface.  The  upper  and  lower 
surfaces  of  the  planform  plus  diaphragm  can  now  be  considered  non- 
communicative.  Thus,  S  in  Equation  (11)  is  presumed  to  include 
any  diaphragm  regions  needed  to  isolate  the  upper  and  lower  surfaces 
but  is  again  restricted  to  that  part  of  the  planform  and  diaphragm 
contained  inside  the  forward  Mach  cone  emanating  from  (x,y,z). 

The  source  distribution  that  is  placed  on  the  diaphragm  is  initially 
unknown,  and  is  evaluated  using  the  boundary  condition  Ap  *  0  at 
ail  points  not  on  the  wing  surface.  This  condition  may  be  expressed 
as 

3  3 

U-M-0  (12) 

Thus  A<f>  **  0  for  an  off-wing  diaphragm  and 

A<j>(x,y ,0, t)  -  A$(xTE,y,0,  t - jj - ) 

for  a  diaphragm  that  lies  in  the  wake  of  the  wing.  For  simple 
harmonic  motion,  this  latter  expression  may  be  stated  thus: 
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The  wake  may  have  an  effect  on  the  wing  (if  the  trailing  edge  is 
subsonic)  or  on  the  tail.  Considerations  have  so  far  been  based 
upon  purely  linearized  theory,  where  the  wake  coincided  with  the 
mean  plane  of  the  lifting  surface.  However,  the  actual  wake  deforms 
in  the  direction  of  disturbed  free  flow  and,  as  is  known,  tends  to 
roll  up.  The  distortion  depends  on  the  spanwise  station,  initial 
angle  of  attack,  flow  Mach  number,  reduced  frequency,  etc.  The 
interaction  effects  on  a  tail  surface  with  elevated  vertical  position 
is  exaggerated  by  the  wing  with  a  positive  initial  angle  of  attack 
due  to  this  distortion  phenomenon.  This  effect  can  approximately 
be  taken  into  account  by  assuming  the  actual  position  of  the  tail 
surface  relative  to  the  steady  wake  position  to  be  the  effective 
vertical  position  to  which  the  unsteady  wake  will  be  superimposed, 
(Ref.  8). 

The  pressure  due  to  the  j  1  mode  of  oscillation  is  found  from  the 
usual  linearized  expression  in  the  x,y,z,  coordinate  system. 
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Thus,  the  airloads  are  evaluated  by  ar.  integration  of  the  pressure 

over  the  area  of  interest.  The  generalized  force  over  the  area  A 
due  to  deformation  in  the  i-th  elastic  or  rigid  body  mode  and  loadings 
of  the  j-th  mode  is 
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where  Ap 


j 


Pressure  difference  between  upper  and  lower  surfaces  at 
point  (x,y)  at  time  t  due  to  a  unit  generalized  de¬ 
flection  of  j-th  mode. 


Further  discussion  of  the  computation  of  pressures  and  generalized 
forces  will  be  found  in  Section  6, 
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3.  NUMERICAL  COMPUTATION  IN  THE  MACH  BOX  GRID  SYSTEM 

If  the  integration  area  may  be  split  into  sufficiently  small  elementary 
areas,  one  may  make  approximations  in  the  evaluation  of  the  double 
integral  in  Equation  (10) . 

The  form  of  elementary  areas  presently  available  for  practical  comp¬ 
utation  are  squares,  the  so-called  Mach  Boxes  whose  diagonals  are 
parallel  to  the  Mach  lines,  and  the  so-called  characteristic  boxes 
which  are  rhombuses  with  sides  parallel  to  the  Mach  lines.  Due  to  the 
fact  that  the  Mach  Box  appears  to  be  the  most  useful  and  versatile  of 
the  elementary  shapes  proposed  to  date,  the  numerical  scheme  for 
evaluating  the  aerodynamic  forces  used  in  this  report  is  based  on 
the  Mach  Box  method  (Refs.  1,  9). 

In  order  to  establish  a  convenient  counting  system  for  the  grid  of 
sending  boxes  a  transformation  is  made  from  the  dimensional  coordinate 
system  to  a  non-dimensional  coordinate  system  (n  ,  m  ,  i  )  in  which 
the  sending  boxes  are  located  at  (v,  y,  0)  whereCv, and  yCare 
integers.  This  transformation,  with  the  wing  as  sending  surface  is: 


“c  ■  «W  -  XCV>'bl  +  1 

“c  “  6Vbl  +  1/2 

*c  "  BVbl 

where  X„,,  is  the  location  of  the  first  box  on  the  sending  surface 
Cw 

(see  Figure  2)  andCXy,  Yy,  Z y)  in  the  wing  local  coordinate  system. 

A  similar  transformation  is  used  when  the  tail  is  the  sending  surface, 
with  Xy,  Yw,  Zw,  Xyc  replaced  by  Xj,  YT ,  Zrj.,  XCT  respectively. 

The  location  of  the  receiving  point,  (x,y,z)  in  the  local  coordinate 
system),  is  designated  (n,m,£)  in  the  (nc»  mc,  i £)  coordinate  system, 
as  shown  in  Figure  4. 

It  is  assumed  that  each  box  in  the  grid  is  sufficiently  small  that 
the  upwash  within  it,  w0^  is  constant.  Thus,  Equation  (11)  may 

be  expressed  in  the  form 
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Receiving  box 


Figure  4  Receiving  Point  Coordinate  Syitem  Correiponding  To 
A  Right  Wing  Sending  Surface 


where 


A(v,u)  denotes  integration  over  box(v,p)  and  b.  is  the  streamwi3e 
dimension  of  a  Mach  box 


A  final  transformation,  which  is 


establishes  a  non-dimensional  coordinate  system  centered  on  the  receiving 
point,  (Figure  A).  The  location  of  a  sending  box  is  (v,p,"X)  in  this 
receiving  coordinate  system. 

The  expression  for  the  perturbation  velocity  potential  (Equation  1A)  may  now 
be  expressed  as 
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(15) 
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and  is  known  as  an  "Aerodynamic  Influence  Coefficient"  or  AIC. 
Also,  £,  n  are  dummy  variables  of  integration 


£  ,  £,  »  n  ,  nT  are  the  limits  of  integration  associated  with 

U  la  11  la 

the  Mach  Box  whose  center  is  at  v,p,X  and  are  in  general  functions 
of  v,p,X 


Also, 


v  *■  n-v 


and 


p  ■  ra-p 

X  =•  £-X  =  £  since  X  must  be  0. 

£  -  Cx  -  c) /b1  n  -  B(y  -  n)/^ 


The  AIC  for  velocity  potential,  c~s“,'y»  la  Equations  (15)  and  (16) 

may  be  interpreted  as  the  velocity  potential  induced  at  an  origin 
of  the  nc>  m  ,  coordinate  system  due  to  a  unit  upwash  distributed 

over  unit  "pulse-sending"  Mach  Box  whose  center  is  located  at 

(v,p,X).  The  "pulse-sending"  Mach  Box  lies  wholly  or  partially 
inside  the  forward  Mach  cone  associated  with  the  point  (0,0,0). 


Three  more  forms  of  Aerodynamic  Influence  Coefficients  are  required 
by  this  three-dimensional  analysis.  They  are  the  AIC's  associated 
with  the  upwash,  sidewash,  and  longitudinal  wash  generated  at  th^ 
receiving  point  and  are  related  to  C~,— ,-r  as  shown  below. 

If  the  perturbation  velocity  potential  at  a  point  (x,y,z)  in  the 
local  coordinate  system  of  a  surface  is  <K*»y,z)»  the  upwash  is  f 

defined  as  the  z  -  component  of  this  velocity  and  is  therefore 
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w(x,y, z)  (x,y,z) 

This_ma^  be  expresaed  in  the  numerical  computation  scheme,  in  the 
(nc>mc,ic)  coordinate  system  using  Equation  (15) 


IJL  „  i.  L±  m  ±  22  w 
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The  AIC  for  upwash  is  therefore  defined  as 

3 

w  -  -  C-  -r* 

vyX  3  -  vyX 


(18) 


which  may  be  defined  as  the  upwash  induced  at  the  pulse-receiving  point 
due  to  a  unit  upwash  distributed  over  the  pulse-sending  box. 

Similarly  for  the  sidewash, 


v(O.O.O) 
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°VP 
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vyX 


and  the  AIC  for  sidewash  is 

v  s  — —  c — t- 
vyX  3  7  Wk 

For  completesess,  the  longitudinal  wash  is 

.  3 

u (0)0,o)  -  -  E  Z  »0  —  c-5I 

1  v  p  vy  3v 

and  the  AIC  for  longitudinal  wash  is 


(19) 


(20) 


(21) 
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U--r  = - C  — -T- 

vyX  3  -  vyX 


(22) 


The  symmetry  of  the  general  wing-tail  configuration  allows  the  assumption 
that  the  receiving  point  is  on  the  right  hand  side  of  the  configuration. 

Furu  ciiisore  need  only  be  evaluated  at  a  planform  or  diaphragm,  and  conse¬ 
quently  the  receiving  points  of  interest  will  be  only  the  centers  of  Mach  boxes 
on  the  right  wing  and  right  tail.  The  possible  combinations  of  receiving  point- 
and  sending  boxes  may  therefore  be  summarized  as  follows: 
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Sending  box  lies  on: 


Receiving  point  lies  on: 


1) 

Right 

Wing 

2) 

Right 

wing 

3) 

Right 

tail 

4) 

Right 

tail 

5) 

Right 

tail 

6) 

Right 

tail 

Right  wing 
Left  wing 
Right  wing 
Left  wing 
Right  tail 

Left  tail 


Each  of  these  combinations  gives  rise  to  a  distinct  coordinate  system, 
(n^,  mc,  &c),  with  its  origin  at  the  receiving  point,  and  with  the 
(nc,  mc)  plane  lying  paralled  to  the  sending  box  plane.  Thus,  in 
Equations  (15)  to  (22)  v  is  always  an  integer,  A  is  the  same  for 
every  sending  box,  and  p  is  a  real  number,  In  cases  1)  and  5) 
above  A  *  0. 


Consider  now  case  3),  the  influence  of  the  right  wing  on  the  right 
tail.  The  wing  and  its  diaphragm  are  the  sending  surface;  the 
receiving  points  are  the  centers  of  the  Mach  Boxes  on  the  tail. 
Figure  5  shows  the  general  relationship  between  the  surfaces. 
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If  P  is  the  center  of  a  receiving  box  on  the  tail,  (nc,mc,£.c)  forms 
the  coordinate  system  for  the  analysis  where  nc  disappears  into  the 
paper  in  Figure  5.  Figure  6  shows  how  the  wing  looks  from  P. 

A  box  is  considered  to  be  a  sending  box  if  its  center  lies  on  the 
wing,  and  if  the  box  lies  v?holly  or  partially  inside  the  Mach 
—2  —2  —2 

hyperbola  K  =  p  +  A  ,  which  is  defined  as  the  intersection  of  the 
wing  plane  with  the  forward  Mach  cone  emanating  from  P, 


Modified  Wing  l.E 


^  (Supersonic) 


Modified  Wing  L.E. 
(Subionic) 

.  DIAPHRAGM 
rft.  BOX 


'.]  MACH  HYPERBOLA 


m 


IBES^Si 


F 


The  size  of  each  AIC  array  is  determined  by  hew  many  boxes  on  the  sending 
surface  can  influence  the  receiving  box.  Receiving  boxes  at  different 
spanwise  locations  on  the  receiving  surface  will,  in  general,  require 
different  AIC  arrays.  However,  receiving  boxes  which  lie  on  the  same  chord 
of  the  receiving  surface  use  the  same  AIC’s,  and  consequently,  an  AIC  array 
that  is  large  enough  to  satisfy  the  requirements  of  the  aftmost  box  on  that 
surface  can  be  used  for  all  boxes  on  that  chord.  Note  again  that  "surface" 
includes  any  diaphragm  areas  needed  in  the  solution  of  the  problem. 

When  the  receiving  surface  lies  in  the  plane  of  the  sending  surface,  the 
AIC  arrays  V--^  and  W--^  are  not  required. 

4.  EVALUATION  OF  AERODYNAMIC  INFLUENCE  COEFFICIENTS 


a.  Evaluation  of  C — r 
vpA 

The  analysis  in  the  previous  section  has  formulated  the  aerodynamic 
influence  coefficients  in  the  Mach  Box  grid  system  with  reference  to 
the  receiving  point  coordinate  system.  Since  the  integrand  of  C — r  in 
Equation  (16)  becomes  singular  at  the  Mach  hyperbola  some  manipulation 
of  C-~  is  needed  to  render  it  amenable  to  numerical  integration.  To 

this  end,  the  function  G(£,n,A)  is  introduced,  where 


G(S,n,A) 


cos [  k  /M  (I- n-A2)1/2) 

— - : — - 

<c2-S2-  j2>1/2 


where  F  is  the  Bessel  function  series  expressed  by 


F(T,n)  -  Jq  (kjT/M)  6  +  £  J2r  (kiT/M)  sin(2r0> 


r  -  1 


where 


t  -  (  c2-  A2)1/2  0 


0  *  sin  ^  n  /  t 


The  significance  oft  and  0  is  illustrated  in  Figure  7. 


Figure  7.  ValuwOf  (r  ,0)  Correiponding To  (I  ,i?) 


Substitution  of  Equations  (23)  and  (24)  into  Equation  (16) 
yields  ~ 
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a  ikl^  t  F(T,nu)  -  F(T,nL)l  dl 


C — r  can  be  evaluated  from  this  expression  using  numerical 


vpX 


integration  techniques,  since  F  is  closely  approximated  by 
the  first  few  terms  of  the  infinite  summation  in  Equation  (24). 


The  limits  of  integration  £u»  nu,  nL  are  determined  by  the 


location  of  the  box  with  respect  to  the  Mach  hyperbola  (see 
Figure  7).  If  the  box  is  cut  by  the  Mach  hyperbola,  only  that 
part  of  the  box  upstream  of  the  Mach  hyperbola  contributes  to 
the  vaiue  of  the  AIC. 


In  general  the  integration  limits  may  be  summarized  in  the 
following  way : 


(25) 


16 


.  Development  of  W-^  ^  ,  and  U--? 


As  written,  Equations  (18), (20)  and  (22)  imply  that  W—  ,V—  and 
u---x  are  obtained  by  differentiating  (>--  .  c--  hwever^is  a 

e!^luated  approximation,  and  is  not  suitable  for  differen- 
tiation  furthermore,  differentiation  tends  to  increase  the  effects 

8a^drurS  and  appr°xi“ations*  Therefore,  the  expressions  for 
grH?ionVpA  vPr  are  raanipulated  t0  a  f°rm  which  requires  only  inte- 


Normal  Wash 


Substituting  Equation_(25)  into  Equation  (18)  yields: 


!i_f  /lki{ 


vSA  "  “  TT  a!  J  e  [  f(t,  nu)  -  f(t,  n  )  ]  dt 


(27) 


Applying  Leibnitz  Theorem  for  differentiation  of  an  integral 


“c;r  -  -  iff  ’  -1 


-ik^ 


e  [F(T,n(J)  -  F(t,n  )])  dl 


-ik  ? 

3  5 

+  f  «  [F(T,0y) 

-  F(t ,n  ) — ; 

3  A 

Now,  T  ■[  I2-  A2 

jl/2  3  j 

Therefore , - 

A  3  T 

= - and - 

3  A 

T  3  5 

3  I 

A  3  t 

Thus  ,  -  =  — 

—  - - 

3  A 

T  3  T 

In  Equation  (28)  , 

since  and  nT  are 

function  of  t. 

-K 


(28) 


3  F(t,7Tu) 

i3 

3  A 

T3  X 

3  F(t,1Tl) 

_  A  8  F<T>\> 

3  A 

1  3  l 

(29) 


18 


i  n -i  infWirf 


Equation  (28)  can  be  integrated  by  parts  as  follows! 


Sidewash 


Combining  Equations  (20),  (24)  and  (25)  yields 


nU  ?U 


f  f 

*  3uJ-  J 


_ -lk  C 

G(5,n,A)  e  dK  dn 


L  F*L 


Since  and  ^  can  not  be  functions  of  p  ,  the  inner  integral  is  not  a 
function  of  p  .  Therefore  by  applying  Leibnitz's  Theorem,  we  have, 


1  3  nu  r  -  -  -  -1Kv 
~  U  J  G(?»nu(A)  e 


“ik,  5 


3nL  f 


}  G(^,nL>A)  e 


Since 


s  k_  999  k- 

~sin (-i  [  r-  n  -  a  ]  )  5G(j;,n,A) 

3£  M  M 


Equation  (33)  integrated  by  parts  gives: 

-ik  F  -  -  _ 

XV  an„  k  _7  _7  _7  1/2.  3nT  k  1/2 

^vpA  -  "  - (^~  V  "  °U  ‘  *  1  ~  Sin7U  "n  "X  1 

-  --  ____ 

.  f  ,,9nU  ,  klr-2  -2  r2,l/2  3nL  ,  klr72  -2  r21l/2\1  +  kl5  ~ikl  \  .7 
+  I  l(r—  sin  —  (5  -n  -A  ]  -r=-  sin  —  [£  -n  -A  ]  ) — — e  )dr. 


The  limit  of  Integration  for  this  expression  are  summarized  in  Equation  (26). 
They  will  yield3riy/3p  =  1  when  does  not  lie  on  the  Mach  hyperbola,  and 
3n.,/3p  a  0  when  n„  lies  on  the  Mach  hyperbola.  It  is  similar  for3nT /3p  . 

It  will  be  seen  that  evaluation  of  Equation  (35)  will  be  troublesome  when  k^=0. 
Consequently,  Equation  (35)  is  remanipulated  in  the  limit  as  k..-*~0. 
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1 


i;.,  ■, 


lim  V-- - =  lim 

rI“°  V"X  *! 


Urn  [-  — 

ol  ’ 


exp(-lk,5)  s1"  — <  5  S  'X  > 


1  ,  72  _“2  -2  ,1/2 

M  “  8  n 


U 


8  T 


a.  1  r,  2  -2  —2.1/2 
sin  ^  (£  **  n^-  ^  ) 


9  n. 


v  f 


5U 


f  1  +  ik  ,C 

-  +|-  t - Zo exp  (-ik.  5) 

C  '*  C  rZ  1 


5l  '  5l  t 


(36) 


{ 


s1-— <  52-  V  »2)1/2  81f  sin  A  (l2.  iiV 

M  d  nU  M  L 


3  p 


k3 


3  n 
3  U 


*~] 

II 


Applying  L 'Hospital's  rule  yields 

-[< 


lim  V--- 

k^O  vyA 


(g2-ng-I2)1/2  3nu  (i2-n2-X2  )1/2  3nL 
C  3u  5 


3k  L 


.(^tJ?.f:t2)1/2±JhL-  _g2_-n2-X2)1/2  3  nL  “I 

j  1  -2  -  -2  _  f  d  - 

7  c  3  V  C  3  w  J 


(37) 


Noticing  that  30y/3w  and  3n^/3p  are  constant  and  by  applying: 


,  2  2,1/2 

f  (  x  -  a  ) 


dx 


,  2  2,1/2 

(  x  -  a  ) 


+  log 


(x2-a2)1/2 


3n 


|l+  a2-n2-72l1/2  ap 
lim  Vvux  “  ~  lo* - 


u 


kj-**  0 


7 1  /?2  -2  -r2. 1/2 1 

C+  (C  -nL-T  )  3- 


(38) 


Equation  (38)  is  used  to  calculate  V-^jfor  the  case  when  k  is  very  small. 
Longitudinal  Wash 

Combining  Equations_(22) , (24)  and  (25)  gives 


3  iu  -nU 


a  u  ^  o 


7  7  T\  ?  jr  jr 


3^l 


X)  e  Is  dn  dc 
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(39) 


ft-n  jr  ft 


-  — . dmcfidtetgteaa  jaaMtiattaasaJiiiiiiiyg 


1 


Since  n  and  n,  are  not  a  function  of  v,  the  inner  integral  is 

U  Li 

not  a  function  of  v.  Thus,  applying  Leibnitz's  Theorem  yields 


35l  r  _ 

— j  6(5  ,n,X)  e 
3v  JnT 


-ikl^L 


dn  ] 


Using  G(F,n,A)  «-r=-F(T,n)  from  Equation  (23) 


Va  ‘  -  V  I  "7  <«’•"„>  -  «A,nL))  .  1  U 

_ 5-  5„  3v 

-ik  5  u 

x{F(T,nn)  -  F(T,n.)} _ e  lL]  ( 

5=  5 

It  should  be  noted  that  use  of  Equation  (41)  in  Equation  (21) 
yields  the  perturbation  longitudinal  wash.  The  free  stream 
velocity  should  be  added  to  it  to  give  the  total  longitudinal 
wash. 


5.  VELOCITY  POTENTIAL  AND  PERTURBATION  VELOCITY  CALCULATIONS 

a.  Zones  of  Influence  in  Wing-Wing  and  Wing-Tail  Interaction 

The  basic  process  of  solving  an  interference  flow  problem  as 
described  in  Ref.  13,  is  one  of  calculating  the  flow  field 
induced  at  one  surface  by  a  source  distribution  on  another 
surface,  cancelling  this  flow  field  by  modifying  the  boundary 
condition  of  the  "interfered"  surface,  and  solving  for  the 
unknown  quantities  by  treating  the  interfered  surface  as 
though  it  were  in  an  undisturbed  uniform  flow. 

This  process  tray  be  more  fully  explained  by  consideration 
the  various  "nor.-comraunicative  zones"  or  disturbed  regions 
of  the  flow  field  and  their  bounding  surfaces,  and  the  various 
physical  conditions  to  be  satisfied  by  source  distributions 
on  these  surfaces.  Figure  9  shows  a  wing  and  tail  combination 
with  subsonic  leading  edges,  both  having  positive  dihedral 
with  the  tail  centerline  located  above  the  wing  centerline. 

(In  our  analysis,  the  tail  may  not  be  located  directly  above 
the  wing,  but  must  be  located  far  enough  aft  to  prevent  mutual 
interference) . 
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Figure  9.  Zones  In  Wing-Tail  Configuration 


The  region  of  disturbed  flow  may  be  divided  into  three  zones e 
as  shown  in  Figure  9.  Zone  (1)  lies  below  the  wing  and  its 
diaphragm,  and  inside  the  Mach  Cone  emanating  from  the  apex 
of  the  wing.  It  is  disturbed  only  by  the  sources  placed  on 
the  lower  surface  of  the  wing  and  wing-diaphragm. 

Zone  (2)  lies  between  the  wing  and  the  tail,  is  inside  the 
wing  Mach  Cone.  It  is  disturbed  by  the  sources  placed  upon 
the  upper  surface  of  the  wing  and  wing  diaphragm  and  the 
lower  surface  of  the  tail  and  tail  diaphragm. 

Zone  (3)  lies  above  the  tail,  and  inside  the  tail  Mach  Cone. 

The  complete  solution  of  the  flow  field,  then,  would  consist 
of  separately  finding  the  velocity  potential  in  these  three 
zones.  It  may  be  seen  however,  that  the  tail  Mach  Cone  lies 
at  least  partially  inside  the  wing  Mach  Cone.  This  has  the 
important  implication  that  the  flow  conditions  ahead  of  the 
tail  Mach  Cone  are  non-uniform.  Hence,  the  usual  methods  of 
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linearized  theory  cannot  be  used  to  find  the  flow  field  inside 
zone  (3),  which  includes  the  upper  surface  of  the  tail. 

An  alternate  formulation  of,  the  problem,  using  modified 
boundary  conditions,  is  thus  used  to  allow  for  this  non- 
uniformity. 

The  boundary  conditions  used  for  the  various  surfaces  and 
diaphragms  may  be  described  as  follows; 


a) 

Wing 

surface: 

The  tangency  condition  of  flow  parallel 
to  the  surface; 

b) 

Wing 

diaphragm: 

Zero  pressure  difference  and  continuity 
of  normal  wash  across  the  diaphragm; 

c) 

Tail 

diaphragm: 

Zero  pressure  difference  and  continuity 
of  normal  wash  across  the  diaphragm; 

d) 

Tail 

surface : 

A  modified  form  of  the  boundary  condition 
which  permits  the  tail  to  be  examined 
as  if  It  were  in  a  uniform  flow. 

1 

We  may  note  that  since  the  non-uniform  flow  affects  both  the 
upper  and  lower  tail  surfaces,  the  modified  boundary  conditions 
are  to  be  applied  to  both  surfaces,  whereas  in  the  zonal 
picture  the  upper  tail  surface  is  isolated  from  the  upper  wing 
surface, 

A  further  illustration  of  this  replacement  of  non—? iniform  flow 
effects  by  modified  boundary  conditions,  showing  how  the  tail 
diaphragms  should  be  treated,  is  given  in  Figure  10. 


The  upper  part  of  Figure  10  shows  a  wing  which  Is  pitched  at  an 
angle  a,  and  a  tail  which  is  parallel  to  the  undisturbed  flow. 
The  influence  of  the  wing  creates  a  pressure  distribution  on 
the  tail  because  the  tail  is  now  placed  in  a  non-uniform  flow. 
This  pressure  distribution  could  be  reproduced  on  the  tail  by 
removing  the  wing  and  deforming  the  tail  surface  appropriately. 
The  appropriate  deformation  could  then  be  represented  by  a 
source  distribution  on  the  tail  planform  upper  and  lower  sur¬ 
faces  which  would  cancel  exactly  the  flow  disturbance  caused 
by  the  wing,  thus  satisfying  the  boundary  condition  on  the  tail. 
This  source  distribution  would  in  turn  give  rise  to  a  distri¬ 
bution  on  the  tail  diaphragm  so  that  the  upper  and  lower  tail 
surfaces  can  be  considered  non-eomraunicative.  There  would, 
however,  be  no  direct  influence  on  the  tai j.  diaphragm  from 
the  wing. 
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TAIL  IN  UNIFORM  FlOW  FIELD 


Figure  10.  Intereetsd  Flow  Field  in  Wing-Tail  Configuration  In  Steady  Ccse 


Once  the  appropriate  regions  of  influence  have  been  determined, 
and  the  necessary  AIC's  calculated  as  described  in  the  previous 
section,  the  strength  of  the  source  distribution  and  the  local 
velocity  potential  differences  may  be  determined.  Since  the 
application  of  this  process  to  the  general  configuration  is 
rather  complex  it  will  be  described  for  three  distinct  problems. 

1)  The  velocity  potential  distribution  will  be  evaluated 
for  a  single  wing  with  dihedral. 

2)  The  distributions  will  then  be  evaluated  for  a  wing  and 
tail  where  neither  surface  has  dihedral,  and  the  tail  does 
not  lie  in  the  plane  of  the  wing. 

3)  These  solutions  will  be  combined  to  give  distributions  in 
the  case  of  the  general  configuration. 
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b.  Solution  for  Single  Wing  with  Dihedral 

This  problem  requires  only  one  receiving  surface  -  the  right 
wing,  and  two  sending  surfaces  -  the  left  wing  and  the  right 
wing.  The  AIC  arrays  required  are  consequently: 


1) 

2) 


3) 

4) 


The  velocity  potential  array  representing  the  influence 
of  the  right  wing  upon  itself,  C- 


Velocity  potential  arrays  representing  the  influence  of  the 
left  wing  upon  the  right  wing,  C-  -  r.  A  separate  array 

V  ,U  ,  A 

is  required  for  each  chordwise  row  of  receiving  boxes  (a 
receiving  chord)  on  the  right  wing; 

IW 

Sidewash  arrays,  left  wing  upon  right  wing  -  - 
one  for  each  receiving  chord;  *  *  * 

(28) 

Upwash  arrays,  left  wing  upon  right  wing  W-  -  r  one  for 

v ,p , A , 

each  receiving  chord. 


The  solution  process  begins  by  considering  the  innermost  box 
on  the  most  upstream  row  of  the  right  wing,  box  (1,1)  in  the 
(r<c,  rac>  £c)  coordinate  system  (see  Figure  11).  Since  the  wing 

leading  edge  may  not  be  swept  forward,  this  box  must  be  on  the 
wing,  and  is  not  influenced  by  any  other  box  on  the  wing  or  its 
diaphragm.  Consequently  the  source  strength  on  the  upper  and 
lower  surfaces  of  box  (1,1)  are  given  by 


p.l 

RUW 


Df}’1 


D  t 


D  i1*1 

NM  =  -_LL 

RLW  D  t 


(42) 


where  if ,n)  and  Np,n)  denotes  the  complex  outward  normal  wash  due  to 

KUW  KLW 

the  source  distribution  on  the  upper  and  lower  surfaces  respectively 
of  box  (n,m)  on  the  right  wing. 

f^,m  is  the  scaled  modal  deflection  at  box  (n,m)  in  mode  j  and 


D  fn,m 


D  t 


1  =_1 


bx  0  f“’m 
if  D  t 


3  f 


n,m 


i  +  blT-J-  ) 
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If  the  leading  edge  of  the  wing  is  at  least  partially  super¬ 
sonic  there  may  be  more  than  one  box  on  the  first  row; 
Equation  (42)  holds  for  all  such  boxes. 

Since  the  configurations  are  symmetric  about  the  center  line, 
and  the  modes  are  either  symmetric  or  antisymmetric,  the 
source  strengths  on  the  left  wing  can  be  found  as 


Here  the  positive  sign  is  taken  for  symmetric  modes,  the  negative 
for  antisymmetric  modes,  and  the  zero  factor  for  the  cases  where 
the  left  hand  surface  contributions  are  to  be  ignored. 
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When  has  been  calculated  for  all  boxes  in  the  first  row, 

the  solution  process  proceeds  to  the  second  row,  then  the  third 
and  subsequent  rows.  Each  row  is  influenced  only  by  boxes  on 
the  rows  upstream  of  it.  Any  row  other  than  the  first  may  in¬ 
clude  diaphragm  boxes,  and  all  rows  except  the  first  will  be 
influenced  by  the  left  hand  side  of  the  planform  and  diaphragm. 

Let  us  suppose  that  box  (n,m)  lies  on  the  upper  surface  of  the 
right  wing  planform.  The  source  strength  placed  on  box  (n,m) 
must  include  the  effect  of  the  modal  deflection  at  (n,m),  and 
also  the  influence  of  the  sources  on  the  upper  surface  of  the 
left  wing  (planform  and  diaphragm)  at  that  box. 


The  tangenty  condition  must  still  be  satisfied  on  box  (n,m), 
i.e.  the  complex  normal  wash  due  to  modal  deflection  and  slope 
must  equal  the  sum  of  normal  washes  induced  by  both  local  and 
remote  source  strength.  (The  only  remote  sources  inducing 
normal  wash  are  on  the  left  hand  surface.)  The  boundary  condition 
may  thus  be  expressed  as 


r\  *  m 

Df ,  n#m  An,ra 

— 1  *  (N  +  N  ) 
Dt  ' HRUW  RUW  ' 

LUW 


(45) 


*n,ra 

where  NRUW  is  the  complex  outward  normal  wash  at  (n,m)  induced 
LUW 

by  the  source  distribution  over  the  upper  surface  of  the  left 
wing,  given  by: 


N. 


,n,i 

RUW 

LUW 


i  _ . 

left  wing 
+  diaphragm 


RW 

W— — r  cos 2<h. 
vpX  rW 


RW 

m  ur 

Vvi^  ain2*W 


The  local  source  strength  is  then  found  by  solution  of  Equation 
(45). 


n,m 

1 

Dfn,m 

ao  |  ro 

nruw 

RUW 

“  Dt 

LUW 

(46) 


It  should  be  noted  that  for  any  given  box  (n,m)  all  the  quantities 
on  the  right  hand  side  of  Equation  (45)  and  hence  of  Equation 
(46)  are  known,  since  the  sura  of  Equation  (45)  includes  only 
boxes  in  rows  ahead  of  the  row  under  consideration,  where  the 
source  strengths  have  been  found  previously. 

A  similar  argument  is  used  to  show  that  if  (n,m)  is  a  box  on 
the  lower  surface  of  the  right  wing 
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Figure  12.  Normal  Warfi  On  Upper  And  Lower  Surface*  Of  Wing 


(47) 


(48) 


Figure  12  illustrates  the  reason  for  the  negative  signs  appearing  in 
in  Equations  (47)  and  (48).  It  must  be  noted,  that  is  here 

defined  by  motion  of  the  upper  surface. 

If  the  box  (n,m)  is  located  off  the  planform  i.e  in  a  leading  or  trail¬ 
ing  edge  diaphragm,  a  different  boundary  condition  must  ’*•>  applied  in 
order  to  find  the  unknown  source  strength.  Since  the  dlarhragm  is  in¬ 
capable  of  supporting  a  pressure  difference; 

ApJ,m  -  0  (49) 

for  boxes  in  any  diaphragm,  where  Ap?,m  is  the  pressure  difference  at  box 
(n,m)  resulting  from  the  motion  of  mode  j. 

The  relation  between  pressure  and  velocity  potential  difference, ^,m, 
is  J 
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.  n DA<f> . 


n,ra 


Dt 


p  ( iw+U— )  A4>j 


n,m 


-  Cilf-ii,  *Tn*ra  i  v  3  .7*1, m  \  iwt 
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where  A^j,ra  is  the  non-dimensionalized  velocity  potential 


(50) 


difference  at  box  (n,m)  due  to  a  unit  generalized  deflection 
of  the  jth  mode.  The  actual  value  of  the  perturbation  velocity 
potential  is  then 


A  4- 


n,m 

j 


£i> 

8  J 


iwt 

e 


(51) 


Setting  the  pressure  difference  to  zero,  and  solving  Equation 
(44)  as  an  ordinary  differential  equation  subject  to  initial 
conditions  at  x  : 

°  X  -  X 

r  -ik  ■ n  °-i 

=n  =,  .  1  b, 

A4.  -  M(xo)[e  1  (52) 

For  leading  edge  diaphragms,  the  initial  value  for  velocity 
potential  is  that  at  the  foremost  Mach  line,  i.e.  zero. 

Thus  for  leading  edge  diaphragm  boxes 

A|t?,m  =  0  (53) 

J 

For  trailing  edge  diaphragm  boxes,  the  initial  velocity  potential 
is  the  value  at  the  trailing  edge  of  the  wing  for  that  span- 
wise  location. 

Thus : 
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(54) 


=.  n,m  s  m 

A  <?  .  *  A  *  TE  exp  [  -ik^x^  -  xTE  Wbj^  ] 

J  m 


The  value  at  the  trailing  edge  is  found  by  extrapolation  from  previously 
calculated  values  as 


4  *  TE  '£<XTE,»  -  V>/bl*4  5  -  4  *  ”  ‘1’”  )  +  4* 


Where  n'  is  the  index  of  the  last  complete  box  on  chord  m.  If  there  is 
only  a  single  box  on  the  chord,  as  in  the  case  of  a  pointed  wing  tip  the 
extrapolation  is  performed  along  a  Mach  line. 

The  non-dimensional  velocity  potential  distribution  is  related  to  the 
local  normal  washes  (source  strength)  through 


=  n,m 


cvp°  (nruw  "  NmP  + 


right  wing 
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.  vp  vp  „ 
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left  wing 
+  diaphragm 


The  source  strength  difference  at  box  (n,m)  on  the  diaphragm  can  now  be 
found  by  solving  Equation  (56).  All  the  source  strengths  ahead  of  box 
(n,m)  are  known,  and  the  velocity  potential  is  given  by  either  Equation 
(53)  or  Equation  (54). 
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An  additional  condition  to  be  solved  is  the  continuity  of 
total  normal  wash  across  the  diaphragm.  This  may  be  expressed 
as 


n,m  ^  n ,  m 
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or  rearranging: 
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Equations  (57)  and  (59)  can  now  be  solved  simultaneously  to 
give  the  source  strength  on  the  upper  and  lower  surfaces  of 
box  (n,m). 

The  source  strengths  found  on  the  upper  and  lower  surfaces 
are  in  general,  different  in  magnitude.  (For  the  case  of 
zero  dihedral,  they  will  differ  in  sign  but  not  in  magnitude, 
in  which  case  only  the  upper  surface  need  be  considered.) 

For  boxes  on  the  planform,  the  velocity  potential  difference 
is  found  by  application  of  Equation  (56). 

c.  Solution  for  Flat  Wing-Flat  Tail  with  Vertical  Separation 

Discussion  of  the  case  of  a  flat  wing  and  tail  with  positive 
vertical  separation  will  lend  some  insight  into  the  process 
of  performing  the  solution  with  two  separated  lifting  surfaces, 
without  the  complicating  features  of  the  dihedral  solution 
discussed  above. 

The  solution  may  be  divided  conceptually  into  three  parts: 

1)  Performing  the  solution  for  the  isolated  wing  and  its  wake; 

2)  finding  the  interference  field  over  the  tail  and  cancelling 
it;  and  3)  performing  the  solution  for  the  isolated,  but 
modified  tail  surface.  Actually,  however,  the  second  and  third 
parts  are  performed  simultaneously. 

The  complete  solution  for  the  wing  is  performed  first.  Since 
there  is  no  dihedral,  the  velocity  potential  difference  is 
just  given  by: 
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For  tail  diaphragm  boxes,  the  zero  pressure  difference  and 
continuity  conditions  are  applied,  and  the  unknown  source 
strength  found  as  for  the  wing.  It  must  be  noted  here,  that 
the  sources  on  the  tail  diaphragm  do  not  include  effects 
directly  from  the  wing,  as  these  are  accounted  for  by  the 
influence  of  the  wing  on  the  tail  surface  itself.  In  the 
absence  cf  effective  motion  on  the  tail  (i.e.  in  the  case 
where  the  interference  normal  wash  just  cancels  the  modal 
induced  normal  wash)  the  source  strength  on  the  tail  would  be 
zero.  Hence  there  would  be  no  tail  diaphragm,  since  it  is 
necessary  only  to  isolate  the  upper  and  lower  tail  surfaces 
in  the  presence  of  a  finite  source  strength.  The  case  of  wing 
dihedral  is  dissimilar  since  the  interference  normal  washes 
on  the  upper  and  lower  surfaces  are  generally  of  different 
magnitudes  while  the  modal  deflections  are  the  same.  Thus 
a  diaphragm  is  always  necessary  to  isolate  the  upper  and  lower 
surfaces,  and  it  is  affected  by  the  opposite  surface. 

Solution  for  Wing  and  Tail  with  Dihedral  and  Vertical  Separation 

The  general  case  is  found  by  combining  the  wing  with  dihedral 
and  wing  tail  solutions  described  above.  First  the  complete 
wing  with  dihedral  solution  is  performed  as  described  in 
Section  (a),  complete  with  a  sufficient  wake  region  to  permit 
tail  interference  effects  to  be  considered.  The  source 
strength  difference  used  to  satisfy  the  boundary  condition 
of  tangential  flow  for  boxes  on  the  tail  surface  is  given 
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Than  the  velocity  potential  difference  is  found  in  the  usu j1  manner, 


35 


fr  a!  in-  - 


,m 


l  C — 

L  vyo 


(NVP 

V  RUT 


rigjtit  tail 
+diaphragm 


-  NVy  ) 
RLT‘ 


,RT} 

VyA 
left  tail 
+&iaphragra 


l  Or 


ID 


(N^ 

U RUT 


-  NVP  ) 
RLT; 


(75) 


For  diaphragm  boxes : 
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which  may  be  solved  for  the  unknown  source  strength  on  the  tail. 
It  may  be  seen  that  the  tail  surface  boundary  condition  (Equation 
67)  includes  direct  wing  contributions,  but  the  diaphragm 
boundary  condition  (76)  does  not. 

e.  Special  Cases 

When  the  dihedral  angle  of  either  surface  is  zero,  or  when  the 
dihedral  angles  of  both  surfaces  are  equal,  or  when  the  tail 
lies  in  the  plane  of  the  wing,  some  simplifications  of  the 
procedure  can  be  made.  There  are  special  cases  which  occur 
quite  often  in  practice  and  are  discussed  briefly  below. 

Flat  Coplanar  Wj.ng-Tail 

In  this  case,  the  two  lifting  surfaces  may  almost  be  treated 
as  one  surface  and  the  source  strength  and  velocity  potential 
found  in  the  same  way  as  for  a  flat  wing  alone.  The  sole 
difference  is  that  there  exists  a  diaphragm  between  the  wing 
and  the  tail  which  is  handled  in  the  usual  manner  of  a  wake. 

It  may  be  theoretically  demonstrated  that  the  three-dimensional 
solution  just  reaches  this  limit  as  L  ■>  0  and  the  tail  surface 
merges  into  the  wing  plane.  One  additional  feature  of  this 
case  is  that  the  "mutual  interaction"  condition,  i.e.  where 
there  is  direct  influence  of  the  tail  on  the  wing,  may  be 
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handled,  since  the  computational  scheme  works  with  a  single 
box  pattern  and  can  conveniently  switch  from  wing  to  tail  and 
back  again  as  required.  This  case  would  occur  with  highly 
swept  coplanar  wing-tail  combinations  in  close  proximity  at 
low  Mach  number. 


Coplanar  Wing-Tail  with  Dihedral 

This  case  is  similar  to  that  described  above,  but  both  the  wing 
and  tail  have  dihedrals,  which  are  equal.  In  this  case,  tbe  two 
surfaces  are  treated  as  a  single  wing  with  dihedral,  and  the 
above  diaphragm  between  the  two  planforms  apply. 

Wing  Tail  with  Equal  Dihedral  and  Vertical  Separation 

This  case  is  similar  to  the  general  case.  Since  the  dihedrals 
are  equal,  some  saving  can  be  effected  by  using  the  same  AIC 
arrays  for  the  effect  of  the  wing  on  itself  and  the  tail  on 
itself.  Separate  arrays  for  the  wing-tail  interaction  effect 
must  be  computed,  however,  Mutual  interaction  is  not  permitted 
in  this  case. 


"■'ng  with  Vertical  Tail 

This  case  is  another  variant  of  the  general  case.  The  tail 
dihedral  is  set  equal  to  90°  and  the  contribution  of  the  left 
tail  surface  is  ignored.  This  case  is  meaningful  only  for  anti¬ 
symmetric  mode  shapes.  Overlap  between  wing  and  tail  surfaces 
is  not  permitted,  so  that  a  T-tail  cannot'  be  analyzed. 


f.  Wing  -  Flowfield  Sampling  Calculations 


The  method  of  interference  flow  fields  described  above  can  be 
used  to  sample  the  perturbation  velocity  components  in  the  flow- 
field,  resulting  from  the  motions  of  the  wing.  The  upwash  and 
sidewash  perturbation  velocity  components  at  a  point  (n,m)  are 
found  by: 
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where  N^j,  N^J  are  defined  above  for  the  general  case,  and  the 
AIC’s  are  the  same  as  those  for  the  wing-tail  effect  witn  ij>_  =  0. 

The  longitudinal  wash  is  found  frotn  the  velocity  potential  in 
the  flow  field 
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Note  that  it  is  most  convenient  to  find  the  values  of  the  pertur¬ 
bation  velocities  at  the  centers  of  "pseudo-boxes"  located  in 
the  flow  field,  rather  than  at  artibrary  points.  Interpolation 
to  different  points  could  be  performed,  of  course. 
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6.  PRESSURE  DISTRIBUTIONS,  GENERALIZED  FORCES  AND  GENERALIZED  AERODYNAMIC 
COEFFICIENTS 

a.  Pressure  Distributions 

Once  the  velocity  potential  difference  is  evaluated  over  a  surface, 
the  pressure  distribution  for  the  jth  mode  is  found  from  the  usual 
linearized  expression  (ref,  3) 
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The  complex  amplitude  of  the  pressure  for  simple  harmonic 
motio.i  is  then 
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In  terms  of  the  nondimensionalized  velocity  potential  difference 
for  box  (n,m),  the  local  pressure  coefficient  difference  is 
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where  the  velocity  potential  difference  A$j(x,y)  is  evaluated 

al  the  leading  and  trailing  edges  of  box  (n,ra)  by  interpolation 
between  the  value  for  box  (n,ra)  and  the  values  for  boxes 
directly  ahead  or  behind.  The  on-planfcrra  area  of  each  box  is 
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where  an,m  is  an  area  factor  used  to  adjust  the  basic  box  area 
along  planform  edges  and  equals  the  proportion  of  box  (n,m) 
that  lies  on  the  planform.  For  interior  boxes,  a  '  is  always 
unity. 

The  lift  on  each  box  due  to  a  unit  generalized  deflection  of  j4 
mode  is  given  by 
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where  the  L^*m  is  the  machine  output  and  is 


(87) 


r  n,m 
1 


[XTE  — — 

[AfjCx.y)]  +  i  an*m  kj?  ^raJ 

XLE 


(88) 


The  lift  on  a  chordwise  section  (chordwise  row  of  boxes)  is  given  by 
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The  total  lift  for  half  the  airplane  is 
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if  the  section  or  total  lift  for  the  wing  pr  tail  alone  is  desired,  the  sums 
of  Equation  (89)  or  (90)  are  taken  only  over  the  appropriate  on-planform  boxes. 

The  lift  and  moment  coefficients  per  unit  span  C  and  C  ,  which  take  account 
of  three  dimensional  interaction  effects,  may  be  round  mdirectly  from  the 
sectional  lift  and  generaliaed  force  by  the  following  relationships: 
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where  c  is  the  local  chord  and  Q  ^  is  the  sectional  generalized 

force  due  to  i-th  pitching  mode  about  a  specified  axis  and  loading  for 

1-th  mode. (see  the  following  section).  ,  =m  =m 

j  uii  &  The  vaiues  L  and  q  are  optional 

machine  outputs.  j  <ij  r 

b. Generalized  Forces 


The  presence  of  the  differentiation  with  respect  to  x  in  the  pressure 
expression  (Equation  85)  makes  it  preferable  to  calculate  the  generalized 
forces  directly  from  the  velocity  potential,  rather  than  by  integration  of 
the  pressure  distributions  weighted  by  the  various  modal  functions. 

This  may  be  accomplished  as  follows.  If  f  (x,y)  is  the  i-th  modal  deflec¬ 
tion  in  rion-dimensional  form,  the-  a&sociatea  generalized  fotce  due  to 
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the  pressure  for  the  j  mode  is 


Q  =  Q 


Qij  “  H  fiCx.yjApjCx.y)  dxdy 


The  applicable  sign  convention  is  that  a  positive  generalized 
force  is  in  the  direction  opposing  the  motion  (restoring  force). 

The  generalized  force  associated  with  each  box  (n,m)  is  therefore 
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Or,  substituting  Equation  (85)  into  Equation  (94) 
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In  terms  of  the  nondimensionalized  velocity  potential  difference 


T  n»m  =  £l- 

3b, 


JJ  £i<*. 

A(n,m) 


3^(x.y) 


y)  [ — - +  ik..A$,(x,y)]  dxdy 

Hb)  j 


Integration  by  parts  and  rearrangement  of  terms  yields 
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Equation  (97)  is  integrated  by  assuming  at  most  linear  variation 
Of  the  parameters  oyer  each  box.  The  box  generalized  load  is  thus 
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The  total  generalized  force  in  the  ith  mode  due  to  a  unit 
generalized  deflection  of  the  jth  mode  is  thus 


(99) 


where 


'ij 


In 

m  n 


(fn,mA=n,m 


lTE  n,m 

»n,o  n,m51  . 


3  -b]A<J> 

XLE 


j 


8x 


ik1an,mf"»mAij,ra  (100) 


C.  Generalized  Aerodynamic  Coefficients 

In  order  to  be  able  to  compare  results  of  the  methods  of  this 
report  to  experimental  values  or  other  theoretical  methods,  it 
is  desirable  to  have  a  standard  notation  and  nondimensionalization 
for  deflections  and  generalized  forces.  Such  a  standard  has 
been  established  by  the  Advisory  Group  for  Aeronautical  Research 
and  Development  (AGARD), 

This  standard  is  based  primarily  on  a  nondimensionalization  by 
semi -span, s  .  In  the  AGARD  notation  system  the  total  general¬ 
ized  force  is  given  by 
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where 
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7.  REFINEMENTS  OF  THE  MACH  BOX  METHOD 
a.  Subdivision 

It  is  well  known  that  the  Mach  Box  Method  may  proauce  inaccurate 
pressures.  The  inaccuracies  take  the  form  of  irregular  fluct¬ 
uations  in  the  chordwise  and  spanwise  pressure  distributions. 

These  pressure  fluctuations  are  more  pronounced  for  a  planform 
which  has  subsonic  leading  edges,  due  to  the  inaccurate  repre¬ 
sentation  of  the  singular  upwas'n  distribution  on  the  diaphragm 
region.  It  has  appeared  in  the  past  that,  in  spite  of  these 
pressure  fluctuations,  generalized  airforces  have  been  acceptable. 

A3  discussed  in  subsection  3  above,  the  true  leading  edge  of  rhe 
planform  i9  replaced  by  the  edges  of  boxes  in  the  Mach  Box  Grid, 
to  form  an  artificial  leading  edges  Unless  the  true  leading 
edge  is  parallel  or  perpendicular  to  the  stream  flow  the  arti¬ 
ficial  leading  edge  will  be  "jagged",  as  in  Figure  13. 
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Figure  13.  Irregularities  in  Distorted  Jagged  Leading  Edges 
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As  may  be  seen  from  Figure  13  ,  leading  edges  which  are  not  parallel 

to  the  Mach  rays  are  represented  by  an  irregular  pattern  of  Mach 
boxes.  It  is  this  irregularity  which  is  the  primary  cause  of  the 
pressure  fluctuation. 

The  concept  of  "distortion"  may  be  introduced  to  quantify  the  ir¬ 
regularity  in  the  relationship  between  the  true  and  artificial  lead¬ 
ing  edges.  This  distortion  is  defined  as  the  maximum  distance 
between  the  true  and  artificial  leading  edges,  measured  either 
spanwise  or  chordvise  (whichever  is  the  larger),  divided  by  the  box 
dimension,  b^or  b^/0,  respectively.  If  the  distortion  does  not 
exceed  1,  the  planform  is  "non-distorted".  Clearly  if  a  finer  box 
grid  is  chosen,  the  distortion  is  unchanged;  this  corresponds  to  the 
finding  that  a  finer  box  grid  merely  increases  the  spatial  frequency 
of  the  pressure  fluctuations ,  but  does  not  reduce  their  amplitudes 
(  Figure  lh  ). 

This  phenomenom  has  been  reported  on  extensively  in  References  10  &  11 
which  found  that  the  pressure  fluctuations  are  principally  caused  by 
(a)thesingular  nature  of  the  AIC's  along  the  forward  Mach  lines 
emanating  from  the  receiving  point, (b)  the  strong  upwash  which  is 
discontinuous  at  a  subsonic  leading  edge  and  in  the  opposite  direction 
on  the  diaphragm  ahead  of  the  edge,  and  (c)  the  irregular  points  * 
produced  by  the  jagged  leading  edge  of  a  grid  of  boxes.  These  effects 
combine  to  cause  noticeable  fluctuations  in  the  velocity  potential 
distributions  which  when  differentiated  cause  the  pressure  distri¬ 
bution  to  fluctuate  unrealistically  over  the  wing  planform. 

It  is  the  aim  of  the  subdivision  technique  to  reduce  the  distortion 
of  the  leading  edge.  This  is  achieved  by  decreasing  the  minimum 
unit  scale  of  the  sending  box  grid  system,  i.e.  the  box  dimension, 
without  changing  the  dimension  of  the  receiving  box,  as  shown  in 
Figure  lh. 


The  "irregular  points"  are  defined  as  the  down-stream  comers 
of  regions  where  the  leading  edge  cuts  more  than  one  box  on 
the  chordwise  or  spanwise  strips. 
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Figure  14  Refinement  Of  Jagged  Leading  Edges 
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Each  box  in  the  sending  box  grid  is  divided  into  Ng  "sub-boxes" 

where  Ng  is  an  odd  integer.  The  receiving  boxes  remain  unchanged, 

their  centers  forming  the  "control  points"  of  the  system.  Clearly, 
each  control  point  will  be  at  the  center  of  the  central  "sub-box". 
Figure  15  shows  that  the  distortion  is  reduced  by  a  factor  Ng, 
and  that  the  effects  of  the  irregular  points  on  the  control  points 
are  effectively  reduced. 
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Figure  15  Avoidance  Of  Irresuler  Point  Effects  By  Subdivision 
In  Odd  Number  Increments 


Principles  for  Application  of  the  Subdivision  Scheme 

The  subdivision  refinement  is  applied  using  the  following  assumptions 

and  techniques: 

1)  Subdivision  may  be  used  only  when  the  receiving  and  sending  boxes 
lie  in  the  same  plane.  This  means  that,  in  the  general  con¬ 
figuration,  subdivision  may  be  applied  only  to  the  calculation 
of  the  effect  of  the  right  wing  on  itself  and  of  the  left  wing 

on  itself.  If  the  wing  dihedral  is  zero,  subdi. vision  may  be  used 
in  calculating  the  effect  of  the  left  wing  on  the  right  wing. 

2)  Subdivision  is  not  applied  to  the  whole  surface.  A  receiving 
box  is  influenced  primarily  by  sending  boxes  which  are  close  to 
it.  Therefore ,  there  is  introduced  the  concept  of  "effective 
area"  shown  in  Figure  lb.  The  effective  area  may  Contain  more  or 
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fewer  subdivided  boxes  than  the  total  number  of  ordinary  boxes 
required  in  the  solution.  The  influence  of  sending  boxes  which 
lie  inside  the  effective  area  on  the  receiving  point  will  be 
determined  using  subdivision;  the  influence  of  those  boxes  on 
the  planform  which  lie  outside  the  effective  area  will  be  cal¬ 
culated  in  the  usual  way.  For  off  planform  boxes  outside  the 
effective  area,  ap  average  value  of  upwash  will  be  used. 


t  • -  CONTROL  POINT  LOCATION 


Figure  16  Subdivision  Of  Box  Grid  System  tn  Effective  Area  Zone 


3)  The  modal  deflection  and  slope  is  calculated  at  the  center  of 
each  unsubdivided  box  that  lies  on  the  planform.  The  deflection 
at  the  center  of  each  sub-box  that  lies  on  the  planform  is 
determined  by  point  slope  extrapolation  from  the  center  of  the 
unsubdivided  box  within  which  it  lies. 

It)  The  upvash  at  the  center  of  each  sub-box  which  does  not  lie  on 
the  planform  is  calculated  using  the  boundary  condition  of  zero 
pressure  difference.  The  upwash  associated  with  any  unsubdivided 
box  that  includes  a  diaphragm  sub-box  is  calculated  as  the  average 
of  the  upwashes  of  the  sub-boxes  which  lie  within  the  unsub¬ 
divided  box. 

5)  Only  the  velocity  potentials  at  the  unsubdivided  box  centers  are 
used  to  calculate  the  generalized  farces  and  pressure  distributions. 

Application  of  Subdivision  to  Flat  Wings-or-Coplanar  Wing  Tails 


The  subdivision  refinement  may  be  applied  only  in  cases  where  the 

sending  and  receiving  boxes  lie  in  the  same  plane.  If  such  a  case 

were  analized  without  subdivision  an  array  of  AIC's  for  velocity 

potential,  C-  -  ,  would  be  used  where  C-  -  is  given 

v,u,o  v,y,o  ° 


C, 


U,M,0 


d  | 


(109) 


for  v  =  o,  1,  ...  N.  where  N.  =  number  of  rows  in  the  array 

and  5  =  0,  1,  2,.,.v  rt 

When  subdivision  is  used  a  second  array  of  AIC's  is  needed.  The 

subdivided  AIC  array  for  velocity  potential  differs  from  the  un¬ 

subdivided  array  only  in  the  value  of  frequency,  k.^  being  replaced 
by  k^/Ng.  Thus  some  saving  can  be  affected  when  multiple  frequencies 
are  considered  by  choosing  them  in  multiples  of  Ng.  Some  differences 
in  the  necessary  sizes  of  the  array  may  be  found,  however,  due  to 
the  finer  breakdown  of  the  planform  geometry  into  a  box  pattern. 

The  nondimensional  velocity  potential  difference  at  the  receiving 
box,  when  subdivision  is  applied,  is 
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This  equation  is  used  to  evaluate  the  source  strength  N  on  the 

diaphragm  by  applying  the  boundary  conditions  Eq.  (  53)  or  (5*4  ) 
for  each  of  the  subdivided  boxes  on  the  diaphragm. 

Application  of  Subdivision  to  a  Wing  with  Dihedral 


If  the  small  (subdivided)  box  center  corresponds  to  a  large  (un¬ 
subdivided)  box  center,  the  unknown  source  strength  for  the  snail 
box  may  be  found  directly  from  the  value  for  the  large  box; 


(n.m) 

RUW 


s 


(N 


n,nu 

RUW' 


(111) 


where  Ns  is  the  subdivision  factor. 

It  can  be  shown  that,  since  both  the  mode  shapes  and  interference 
normal  washes  are  obtained  for  large  box  centers  only,  that  the  un¬ 
known  source  strength  for  small  boxes  on  the  planforra  rfhich  do  not 
correspond  to  large  box  centers  may  be  found  directly  from  the  value:, 
at  large  box  centers  by  a  point -slope  extrapolation.  Assuming  the 
complex  source  strength  represents  the  normal  wash  produced  by  a 
pseudo  mode  shape, 


(n,m) ' 
RUW 


Ns  |  Re  tNRuvj]+  1  [Im  ^Nruw^  +  '  x)  Re  [nr^]_  I 

;  (lid; 
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where  xg  is  the  X  coordinate  of  the  sub-box  center.  A  similar 

expression  can  be  deduced  for  N^m)  on  the  lower  surface. 

If  the  small  box  is  not  contained  within  a  large  box  on  the 
planform: 

*")S  ■  ?,  tC-”l 

- 1 11”  +  -  *>■>*  i]  ttuj 

where  n  +  p  is  the  index  of  the  nearest  on-planform  large  box, 
upstream  or  downstream. 

The  velocity  potential  for  on-planform  unsubdivided  boxes  is  then 
found  as : 


n,m 


<N 
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+  ^  cvpo  ^nruw  "  nrlw  ^  (114) 

right  wing  +  diaphragm 
outside  the  effective  area 


left  wing  +  diaphragm 
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For  subdivided  boxes  in  leading  or  trailing  edge  diaphragms,  the  local 
source  strength  is  found  by  satisfying  the  zero  pressure  difference 
and  continuity  conditions  explicitly  at  every  point  (rather  than 
extrapolating  from  large  box  values).  Thus  in  the  diaphragm  region 


N(n,m)  N(n,m)  _1_ 

RUW  "  wRLW  s 

000 
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Here  the  interference  normal  wash  values  are  again  found  by  the  point- 
slope  extrapolation  method  from  large  box  center  values. 
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The  velocity  potential  difference  for  subdivided  diaphragm  boxes  is 
found  in  a  similar  manner.  For  leading  edge  diaphragm  boxes: 
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for  trailing  edge  diaphragm  boxes : 


rf-xf 


Af["-n)S  =  A$"”* 


(120) 


where 


)(  A?*”'- A?  ?«'")  +AQ*  ■"  (1?1) 


where  n  is  the  index  of  the  last  complete  box  on  chord  m. 

Once  all  the  on-planform  and  diaphragm  subdivided  source  strengths 
are  determined,  the*  velocity  potential  for  unsubdivided  boxes  is 
found  from  Eq.  (l"4  '. 


b.  Velocity  Potential  Smoothing 

The  subdivision  technique  described  above  reduces  substantially  the 
fluctuation  in  the  velocity  potential  and  pressure  distributions. 

Some  fluctuations  still  occur  however,  for  the  subdivision  technique 
is  still  a  finite  element  approach.  These  smaller  fluctuations  may 
be  eliminated  through  the  application  of  a  least  squares  surface 
fitting  technique  to  the  velocity  potential  distribution.  The  method 
is  described  as  follows. 

The  function  chosen  to  represent  the  velocity  potential  is 
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where 


x  = 


X 


y  = 


y 

V* 


n,m  =  integers  ;  n  +  m  ^  5 


f  ,  f  =  sets  of  least  squares  polynomial  coefficients 
used  to  fit  the  real  and  imaginary  parts, 
respectively,  of  the  velocity  potential  function. 


Application  of  the  least-square  surface  fit  to  the  velocity  potential 
functions  is  performed  separately  over  each  lifting  surface.  The 
number  of  control  points  required  (including  the  leading  edge  points 
where  the  velocity  potential  is  zero  or  determined  from  the  wake)  is 


Number  of 

Control  Points  >  S  i  for  a  k-th  order  polynomial 

i-1 


In  practice,  however,  the  control  points  used  for  a  given  order 
polynomial  should  exceed  these  theoretical  limits  by  at  least 


k+1 

1/2  Si  in  order  to  obtrin  a  smooth  fit. 

ial 
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c.  Thickness  Corrections  to  the  Pressure  Distribution 

Tn  the  realm  of  linearized  flow  solutions  without  interference 
effects,  the  pressure  distribution  due  to  thickness  may  be  accounted 
for  separately  and  combined  by  superposition.,  The  thickness  effect 
is  more  critical  to  flutter  stability  when  the  flow  Mach  number  is 
close  to  one. 

When  interference  effects  are  included,  the  thickness  problems  may 
not  simply  be  decoupled.  For  simple  harmonic  motion,  however,  thick¬ 
ness  alone  would  generate  an  oscillatory  loading  only  when  the  upper 
and  lower  surface  modal  displacements  have  the  opposite  sign.  Since 
this  is  an  unlikely  case,  thickness  has  been  ignored  in  the  methods 
for  interference  calculations  used  in  this  report.  However,  a 
correction  on  the  pressure  distribution  for  thickness  can  be  per¬ 
formed  locally  by  applying  second  order  piston  theory  (Ref.  12). 

This  means  that  the  local  pressure  is  uniquely  determined  by  the 
local  velocity  normal  to  the  lifting  surface. 

The  second  order  piston  theory  may  be  expressed 


P-  Po» 


?  a2 

O*  a> 
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III 

4 


(123) 


In  this  expression  the  velocity  w  is  normal  to  the  oscillating 
surface  and  p  is  the  local  static  pressure. 

For  non-planar  wings,  the  shape  of  lifting  surface  is  expressed  by 

ZU  (x,y,t)  =  ZM  (x,y)  elwt 

\  (x.y.t)  =  ZM  (x,y)  elwt 

where  Z  (x,y)  is  the  thickness  distribution  about  the  mean  surface 
which  is  represented  by  the  mean  surface  distribution,  Z^  (x,y)  as 
shown  in  Figure  18. 


4  \  (x,y) 
-  \  (x,y) 
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Z,  (x,y) 


Figure  18  Thickness  Distribution  About  Mean  Surface 


Since  surface  pulsing  is  not  allowed,  Z^x.y)  is  not  time-dependent. 
Therefore  the  normal  velocity  at  the  surface  Z  (x,y,t)  is 


W(x,y,'t)  = 


3  2 

W  +  u 


=  .  iW?„<x,y  )  ]  e”'  ± 


where  (+)  ana  (-)  correspond  to  the  upper  and  lower  surface  respect¬ 
ively  and  Wq  (x,y)  is  the  amplitude  of  the  mean  surface  motion. 

By  substituting  the  non-time  dependent  term  of  Equation  (125)  into 
Equation  (123),  the  pressure  over  the  upper  surface  is 
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Similarly  the  pressure  over  the  lower  surface  is 

[p-  p„]u  =  -?«».[ »- u  It1] 


(127) 


'Hie  net  pressure  difference  between  the  upper  and  lower  surfaces 
is  thus 
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The  first  term  in  Equation  (128) is  equivalent  to  the  1st  order 
piston  theory  which  is  computed  over  a  planar  wing.  Therefore, 
the  thickness  correction,  Z  (x,y),  to  be  applied  to  the  pressure 
difference  is 

w  If  +  1  Z  -r* 

-j— m  ^ 


The  expression  for  pressure  at  a  point  is  thus  modified  by  the 
thickness  correct  ion  factor  and  is 
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The  local  pressure  coefficient  difference  becomes 
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and  the  non-dimensionalized  box  lift  becomes 
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The  non-dimensionalized  generalized  £grce  in  the  i  mode  due  to  a 
unit  generalized  deflection  of  the  j  mode  is  thus 
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SECTION  III 


COMPUTER  PROGRAM  USAGE 


MACHINE  AND  OPERATING  SYSTEM  REQUIREMENTS 

The  program  is  written  in  FORTRAN  IV  for  execution  on  the  CDC  6600 
under  the  SCOPE  operating  system.  Title  FORTRAN  coding  is  compatible 
with  the  RUN  compiler  and  with  FORTRAN  Extended.  Using  the  overlay 
feature  of  the  SCOPE  operating  system,  a  field  length  of  about 
132000g  is  required  to  load  the  program,  and  about  12H000q  for 
execution.  Disk  space  or  tape  drives  are  required  for  five  internal 
binary  scratch  files,  in  addition  to  the  standard  INPUT  and  OUTPUT 
files,  and  further  optional  binary  input  and  output  files  are  allowed 
for . 

TIMING  AND  OUTPUT  ESTIMATES 
a.  TIMING 

Three  measures  of  time  are  available  for  most  CDC  operating 
systems;  Central  Processor  (CP)  time,  Peripheral  Processor 
(PP)  time,  and  Occupancy  time.  Of  the  three,  CP  time  is  most 
consistent,  since  the  program's  usage  of  the  Central  Processor 
is  not  dependent  upon  the  job  mix  or  system  configuration.  The 
PP  times  for  two  runs  of  the  same  job  may  vary  widely,  since 
Periph „ral  Processors  may  be  held  up  while  waiting  for  other 
10  devices  used  and  the  operating  system  driving  them.  Occupancy 
(clock)  time  is  the  least  consistent  of  the  three,  since  it  is 
also  dependent  on  the  priority  levels  of  other  jobs  in  the  system. 
With  this  in  mind,  CP  time  usage  is  presented  in  this  report  in 
some  detail,  and  rough  PP  and  Occupancy  times  are  given  only  to 
present  an  idea  of  order  of  magnitude . 

CP  time  in  seconds  can  be  estimated  from: 

T  =  60  to  70  for  UPDATE,  compilation  and  load 

+2  to  7  for  Geometry  and  Modal  input 

ft  of  K- values 

+  (time  for  AIC  calculations,  see  following  description) 
i  =  1 

+  (ti  of  modes)  x  {#  of  K-values )  x  (time  for  Velocity  Potentials 

+  time  for  Generalized 
Forces,  see  following 
descriptions ) 
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The  time  necessary  to  do  a  full  UPDATE,  compilation  and  load 
of  the  program  is  roughly  60-70  seconds  CP,  200  -  300  seconds 
PP,  and  U  to  7  minutes  occupancy.  Most  of  this  can  be  avoided 
if  the  program  is  maintained  in  absolute  form  on  tape  (or  disk). 

Geometry  setup  and  modal  input  for  a  typical  case  ranges  from 
two  to  seven  seconds  CP  time. 

The  time  involved  in  calculating  AIC's  is  variable,  depending 
upon  Mach  number,  reduced  frequency,  geometry,  and  desired 
integration  accuracy.  In  general,  time  increases  with  K-value 
and  decreases  with  increased  Mach  number,  for  a  constant  number 
of  rows  on  the  box  pattern.  A  completely  flat  configuration 
requires  one  AIC  array  if  no  subdivision  is  desired,  and  two  with 
subdivision.  A  flat  wing  and  flat  tail  with  vertical  separation 
require  two  sets  of  arrays.  Dihedral  angles  require  in  addition, 
one  set  per  chord  of  the  tail  for  the  wing/tail  interaction. 

Thus,  a  typical  configuration  having  30  chords  on  the  box  pattern 
requires  one  AIC  when  completely  planar,  2  AIC  arrays  if  vertical 
separation  is  specified,  and  between  50  and  70  if  either  surface 
has  dihedral.  The  timing  involved  is  not  directly  proportional 
to  the  number  of  AIC  arrays  needed,  since  spatial  AIC's  are 
typically  sparse  and  various  sizes,  so  a  rough  rule  of  thumb  is 
to  multiply  the  time  estimated  for  the  required  planar  AIC  by  .5 
for  the  time  needed  per  spatial  AIC  set.  The  desired  integration 
accuracy  parameter  plays  an  important  role,  since  an  accuracy  of 
.01%  (EX AIC  =  .TRUEj  takes  roughly  5  1/2  times  as  much  CP  time  as 
the  default  accuracy  of  1.055.  Figure  19  indicates  some  timing 
results  from  the  application  phase  executions. 
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By  far,  the  most  time  consuming  part  of  the  program  is  the  section 
which  calculates  normal  washes  and  velocity  potentials.  Without 
subdivision,  CP  seconds  per  mode  can  be  given  roughly  as 


1  2 

■  Y^'ooo —  X  (#  of  boxes)  ,  see  Figure  20  .  The  number  of  boxes 

includes  all  diaphragm  boxes  for  both  wing  and  tail.  If  subdivision 
is  specified  (SUBDV  =  .TRUE,)  the  question  becomes  more  complex. 

The  size  of  the  effective  area  controls  the  timing  for  planform 
box  computations,  and  diaphragm  boxes  take  roughly  nine  times 
as  much  time  as  planform  boxes.  A  few  examples  are  given  in  the 
table  below: 


BOX  PATTERN 

Effective 

Area, 

Rows 

CP  Seconds  per 
Mode 

ft  Rows 

H  Chords 

H  of  Boxes 

Planform 

Total 

Planfom 

Diaphragm 

37 

10 

2k 

2U2 

28U 

8 

268 

3  6 

6 

21 

1* 

39.9 

37 

6 

22 

150 

300 

default* 

175. 1»71 

12 

0.1*13 

*  Default  is  1/3  (H  of  rows)  or  600  sub-boxes,  whichever  is  smaller. 

FIGURE  20  SAMPLE  TIMING  OF  VELOCITY  POTENTIAL  CALCULATIONS 
WITH  SUBDIVISION 


The  time  spent  calculating  Generalized  Forces  is  a  function  of  box 
pattern  size,  number  of  modes,  and  number  of  reduced  frequencies. 
The  range  of  CP  time  is  from  0,7  seconds  to  2.1  seconds  per  mode 
for  each  reduced  fre^UGncy. 

b .  Output 

Output  can  range  from  a  minimum  of  about  two  pages  per  reduced 
frequency,  for  a  small  cuse  with  no  print  options,  to  around 
100  pages  per  mode  for  each  reduced  frequency,  for  a  large  case 
with  all  print  options  turned  on. 


KIwUKj:  .  1 
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PROGRAM  LIMITS  AMD  RESTRICTIONS 


The  primary  limit  on  the  program  is  the  number  of  Mach  boxes  which  can  be 
handled.  This  is: 

(l)  If  a  single  surface  (  or  a  co-planar  wing-tail)  without  dihedral  is 
being  analyzed  without  subdivision,  the  box  pattern,  including  tip 
diaphragm,  may  extend  50  boxes  fore-to-aft  and  50  boxes  root-to- 
tip.  (See  Figure  22) 


FIGURE  22  Determination  of  Box  Pattern  Limits 

(All  measurements  are  in  number  of  boxes) 

(2)  For  any  spatial  configuration  or  sampling  case,  a  hO  X  40  box  pattern, 
for  wing  and  tail  combined  (with  diaphragms)  is  the  maximum. 

( j)  Total  number  of  on-planform  boxes  (control  points)  £  1000. 

(It)  Total  number  of  boxes  considered,  planfom  and  diaphragm,  <  l.’V1'. 
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Geometric  configurations  that  fit  the  program  analysis  capabilities  are 
shown  in  figures  25,  26,  27  and  29.  Figures  28  and  30  illustrate  some 
of  the  configurations  not  allowed.  The  planforms  must  be  defined  in 
terms  of  straight  line  segments  and  may  take  on  various  shapes  as 
indicated  in  figure  25.  The  wing  and  tail  may  have  dihedral,  vertical 
separation  and  horizontal  separations.  Aerodynamic  interaction  of  the 
tail  affecting  the  wing  may  be  evaluated  only  for  the  coplanar  case. 

The  tail  may  physically  overlap  the  wing  for  the  spatial  case'  however 
the  tail  aerodynamic  "Mach  cone"  must  not  intersect  the  wing  as  shown 
in  figure  29. 

Each  planform  edge  must  be  defined  by  10  points,  located  by  (x,  y) 
in  the  local  coordinate  system  (fig.  25).  Each  edge  must  start  with 
y  =  0,  and  the  y  value  for  the  outmost  leading  edge  point  and  trail¬ 
ing  edge  point  must  be  identical  for  each  planform.  The  y-values 
increase  outward,  and  x-values  increase  aft. 

The  leading  edge  of  any  planform  cannot  have  any  forward  swept  sections. 
The  trailinp  edge  may.  (Figures  25,  oO) 

Planform  edge  definitions  must  be  monotonic  in  y;  i.e.,  no  planform 
edge  may  "double  back"  in  the  y-direction.  Segments  of  an  edge  may  be 
parallel  to  the  x-axis,  however,  other  considerations  may  cause  bad 
answers  for  planforms  of  this  nature.  (Figure  25) 

The  dihedral  angles  and  ^  must  lie  in  the  range  0°  ^  ^  <  45°  or 
T  =  90° .  (Figures  26,  27) 

A  maximum  of  20  modes  is  allowed. 


Number  of  reduced  frequencies  20. 


Mach  number  >1.0 
Mach  number  <5.0 


(recommended  range  from  1.2  to  3.0) 


For  printing  purposes  it  is  assumed  that  reduced  frequencies  based  on  b. 
will  lie  within  the  range  lO-5^^^  1.0  or  k  =  0.0,  and  reduced 
frequencies  based  on  semispan  s  will  lie  within  the  range 

.  5  . _  i .  ^  i  nnn  ft  _ _  ’  ""  ft 

10  >  t\  v.  i uw.u  or  k  =  u.u. 

N  S  'v  s 
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DECK  SETUP 


a.  Deck  Setup  for  running  the  program  from  a  SCOPE  3  UPDATE  tape. 
(  Job  Card  CM132000. 


Control  / 
Cards  | 


REQUEST  TAPE.  (Tape  containing  OLDPL) 
REWIND (TAPE) 

COPYBF( TAPE, OLDPL) 

UNLOAD(TAPE) 

UPDATE(F) 


RUN(S,, 1022, COMPILE) 
LGO. 


End-of-Record  Card 

(Optional  Correction  Sets) 
End-of-Record  Card 

(Card  Data  for  the  program) 
End-of-File  Card 


V 


b. 


Deck  Setup  for  running  the  program  ns  above  with  previously  generated 
AIC  arrays  saved  on  tape. 


I  Job  Card  CM132000. 


Control 

Cards 


REQUEST  TAPE.  (Tape  containing  OLDPL) 

REWIND (TAPE) 

C0PYBF( TAPE , OLDPL ) 

UN LOAD (TAPE) 

UPDATE(F) 

RUN(S,, 1022, COMPILE) 

REQUEST  PAICS.  (Tape  containing  planar  AICs ) 

REQUEST  SAICS.  (Tape  containing  spatial  AICs) 

REWIND( PAICS) 

REWIIID(  SAICS) 

C0PYBF( PAICS ,TAPE2 , 3 ) 

OOPYBF( SAICS ,TAPE3  >  3 } 

UN LOAD (PAICS) 

UNLOAD (SAICS) 

LGO. 


End-of-Record  Card 

(Optional  Correction  Sets) 

End-of-Record  Card 

(Card  Data  for  the  program  ,  0AIC=2,  0SAIC=3) 
End-of-File  Card 
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5.  CARD  INPUT  DATA 


a.  Namelist  Format 

Several  of  the  card  sets  defined  in  this  input  data  description 
are  in  NAMELIST  format.  The  following  few  paragraphs  explain  this 
format  as  used  here.  Beginning  in  column  2  is  the  quantity  $LNAME, 
where  LNAME  is  a  name  associated  with  a  particular  card  set.  The 
LNAME  for  card  set  C  is  "CARDC11  therefore  Card  Set  C  must  begin  with 
$CARDC  in  column  2. 

$LNAME  is  followed  by  a  blank  and  a  string  of  specifications, 
separated  by  commas  and  terminated  by  a  $.  Each  specification  is 
one  of  the  following  forms : 

(1)  VNAME  =  VALUE 

(2)  VNAME(l)  =  VALUE1,  VNAME(2)  =  VALUE2 ,  VNAME(3)  =  VALUES,  etc. 

where  VNAME  is  one  of  the  legal  variable  (or  array)  names  for  the 
card  set,  and  VALUE  is  the  value  (or  values)  associated  with  VNAME . 

In  form  ( 1) :  (a)  VALUE  may  be  of  the  form  .TRUE.,  .T.,  .FALSE., 

or  .F.  when  VNAME  is  a  logical  variable. 

(b)  VALUE  should  be  an  integer  if  VNAME  is  an 
integer  variable. 

(c)  VALUE  should  be  a  real  number  if  VNAME  is  a  real 
variable.  Examples:  -1.0,0.07,  3.5E2. 

In  form  (2):  (a)  VNAME  must  be  an  ar  .y  name  and  the  values  given 

will  be  stored  in  ti.e  specified  array  locations. 

The  field  on  a  card  consists  of  column?  ?  through  80.  The  string  may 
be  continued  from  card  to  card  with  i  .  •  restriction  that  a  specifi¬ 
cation  may  not  be  split  between  cards.  Embedded  blanks  are  allowed 
except  within  $LNAME,  within  VNAME  and  *.  .thin  VALUE.  The  order  of 
appearance  on  a  card  is  not  important.  <V  variable  for  which  a  de¬ 
fault  value  is  desired  need  not  be  specj "ied  unless  a  previous  pas, 
through  the  program  has  replaced  the  default  value. 


b.  Data  Description  Conventions 


For  conciseness,  if  two  or  more  parameters  have  similar  meanings, 
they  are  bracketed  together  in  the  description.  Norma], ly  each 


bracketed  parameter  may  be  independently  specified, 


PRPAI^l 

prsai;/ 


In  a  few  cases,  two  methods  of  specifying  the  same  information  are 
provided.  In  these  cases,  the  parameters  are  bracketed  with  the 


word  "or" ,  as  in 


0AIC\ 

NAIC  J •  ^  only  one  of 


the  bracketed  paramet'-r.. 


should  be  specified. 
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Data  Descriptions 


Linkage  Card  to  Execute  the  Program  Format-(A6,l*X,2ilO) 

6  20  30 

- i - 1  — i — ■  ~ 

AFMBOX  i  0 

AFMBOX  =  Program  Overlay  File  Name 
1  =  Primary  level  number 

0  =  Secondary  level  number 

Card  A  -  Title  Card  (Once  per  cycle)  Format- (8A10) 

(  TITLE 

TITLE  -  Columns  1  to  80  are  printed  as  page  header 
information. 

Card  B  -  Mach  Number  (Once  per  cycle)  Format- (NAMELIST) 

( $CARDB  XMACH  = $ 

XMACH  =  Mach  number  for  the  airforces  of  this  cycle.  All 

calculations  for  a  cycle  are  done  at  one  Mach  number. 
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Card  Set  C  -  Control  Parameters 

All  of  the  following  control  parameters  have  default  values  (underlined). 
Any  parameter  not  specified  will  take  its  default  value.  For  a  recycle, 
see  parameter  DEFAULT  below. 


Format (NAMELIST) 


'  $CARDC 

•  • 

$ 

Parameter 

Type 

Optional 
Values , 
Default 
Underlined 

Meaning 

DEFAULT 

Logical 

.TRUE. 

All  parameters  for  this  card  set.  are  set 
to  default  values.  Any  other  parameter 
on  this  card  set  will  be  ignored. 

.FALSE. 

All  parameters  listed  on  this  card  set 
will  be  used.  If  this  is  the  first  cycle 
through  the  program,  any  unspecified 
parameters  will  take  their  default  values. 

If  this  is  not  the  first  cycle,  unspecified 
parameters  will  take  the  values  they  had  in 
the  previous  cycle. 

PRVGEOM 

Logical 

.TRUE. 

Use  geometry  defined  in  the  previous  cycle, 
this  cycle  being  a  change  in  Mach  numbers 
and/or  modes.  Cards  F-L  are  not  needed. 

Not  allowed  in  the  first  cycle  of  the 
program. 

.FALSE. 

New  planform  geometry  is  to  be  speoifi  1 
for  this  cycle. 

PRVMODE 

Logical 

.TRUE. 

Mode  shapes  from  the  previous  cycle  are  to 
be  used.  Not  allowed  in  the  first  cycle 
of  the  program,  or  if  geometry  or  Maci 
number  has  changed. 

.FALSE. 

New  mode  shapes  are  to  be  specified  : or 
this  cycle. 

SYM 

Real 

1.0 

Symmetric  analysis 

-1.0 

Antisymmetric  analysis 

0.0 

No  left  hand  surface  contributions  are 
to  be  calculated. 
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Card  Set  C  -  continued 


Parameter 

Type 

Optional 
Values , 
Default 
Underline  d 

Meaning 

NSURF 

Integer 

1  or  2 

Number  of  surfaces  to  be  analyzed. 

fMTYPEW’l 

(mtypetj 

Integer 

Type  of  modal  input  for  the 

fwingl 

|tailf* 

1 

Polynomial  coefficients 

2 

Arbitrary  locations 

3 

Box  center  values  (Mach  number 
dependent) . 

SUBDV 

Logical 

.TRUE. 

Apply  subdivision  to  the  normal-wash 
calculations. 

.FALSE. 

Do  not  apply  subdivision,  (Basic 
analysis) . 

N ROWE A 

Integer 

0  to  12 

If  subdivision  is  .TRUE. ,  this  parameter 
allows  the  user  to  specify  the  maximum 
number  of  unsubdivided  rows  in  the  effec¬ 
tive  area.  If  0,  the  program  uses  1/3  the 
total  number  of  rows  specified.  (See 

Figure  16) 

EXAIC 

Logical 

.TRUE. 

Calculate  AIC's  with  an  integration 
tolerance  of  0.01%. 

.FALSE. 

Calculate  AIC's  with  an  integration 
tolerance  of  1.0$. 

fDIHW'l 

(dirt/ 

Logical 

/ 

.TRUE. 

Include  dihedral  in  calculating  the 
influence  of  the  (Wing)  on  itBelf, 

(Tail) 

.FALSE. 

Use  any  dihedral  only  in  the  calculation 
of  interaction  between  the  planforms. 
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Card  Set  C  -  Continued 


ISMPLW 


’SMOOI'U  "\ 

- - “***  I 


NDEG 


PLYWOOD 


PRCM 


Integer 


T  A  r>  1 

AA/  gXUCU. 


Integer 


Logical 


[Logical 


Optional 
Values , 
Default 
Underlines 


0  to  10 


.TRUE. 


.FALSE. 


0  to  10 


.TRUE. 


.FALSE. 


.TRUE. 


.FALSE. 


Meaning 


Number  of  chordB  for  which  sampling 
of  upwash,  sidewash  and/or  longitudinal 
wash  is  desired  off-planfora.  (See 
Card  G) . 

Smooth  the  velocity  potential  differences, 

■rn  in  Cswrf&cc 

4^  *  ,  by  a  ^  or  V  least  squares 

ichordwisej 

polynomial  fit  before  they  are  used  in 
calculation  of  generalized  forces. 
Generalized  forces  for  smoothed  and  un¬ 
smoothed  are  calculated. 

Do  not  smooth  the  values;  calculate 

generalized  forces  for  only  the  unsmoothed 
1$  values. 

The  maximum  degree  of  least  squares 
polynomial  to  apply  for  smoothing 
(0  causes  the  program  to  determing  the 
degree  as  a  function  of  the  number  of 

f  surface')  . 

boxes  on  the  ^  or  f  '  • 
chord  J 


Calculation  of  lifts  and  generalized 
forces  is  to  be  done  using  full  box 
areas  everywhere. 

Calculation  of  lifts  and  generalized 
forces  is  to  be  done  using  planform 
defined  box  areas. 

Print  Q  ,  sectional  generalized  force  ror 
mode  1.  JIf  first  mode  is  a  pitch  mode,  then 
sectional  moment  coefficients  may  be  o\  ’.ain- 

ed  from  Cm.  =  Q^./c^.  (See  Sect.  II. 6,  eqn. 
91) .  mj  ij 

Do  not  print  Sectional  generalized  forces 
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E> 


Card  Set  C  -  Continued 


Optional 

Parameter 

Type 

Values , 
Default 
Underlines 

Meaning 

PRBOX 

Logical 

.TRUE. 

Print  box  code  pattern(s). 

.FALSE. 

Do  not  print  box  pattern(s). 

CprcoefI 

\PRM0DSj[ 

Logical 

.TRUE. 

Print  "Odaljjefiertiona*  ]  used-  lf 
and  slopes 

available. 

.FALSE. 

e  coefficients  1  , 

Do  not  print  f  deflections  j  used,  ’j 

and  slopes  } 

K 

| PRPAIC^ 
^rRSAICj 

Logical 

TRUE. 

Prlnt  AIC  arraJr(s)'  f 

.FALSE. 

Do  not  print  AIC  array(s).  * 

(PRUw' 

/PRSW 

PRLW 

Logical 

.TRUE. 

/'upvashes,  v/U  "J  ; 

,  J  sidewashes,  v/U  V  -  i 

Print  |  longitudinal  washes,  u/Uj 

4 

sample  chords. 

.FALSE. 

( upwashes  \ 

_  ,  ,  j  sidewashes  V 

Do  not  Prin  | longitudinal  washes J 

PRNW 

Logical 

.TRUE. 

Print  normal  vashes(N^y^,  etc.)  due  to  ( 

local  source  strength. 

.FALSE. 

r 

Do  not  print  normal  washes.  t 

PRVP 

Logical 

.TRUE. 

Print  velocity  potentials,  A$n,m. 

.FALSE. 

Do  not  print  velocity  potentials. 

} 
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Card  Set  C  -  Continued 


Parameter 

Type 

Optional 
Value , 
Default 
Underlined 

Meaning 

PRBL 

Logical 

.TRUE. 

Print  box  lifts,  L^’1”  , 

•FALSE. 

Do  not  print  box  lifts. 

PRDCP 

Logical 

.TRUE. 

Print  pressure  difference 

-n,m 

coefficients  Ac 

.FALSE. 

Do  not  print  pressure  difference 
coefficients. 

PRSL 

Logical 

.TRUE. 

Print  section  lifts,!)1” 

.FALSE. 

Do  not  print  section  lifts. 

PRGNAF 

Logical 

.TRUE. 

Print  generalized  air  forces. 

% 

.FALSE. 

Do  not  print  generalized  air  forces. 

PRGNAO 

Logical 

.TRUE. 

Print  generalized  aerodynamic 
coefficients,  and  Q ,  if 

PRGNAF  =  .TRUE. 

***WARNING  -  To  obtain  generalized  aerodynamic  coefficients  in  the  correct 
AGARD  notation  the  planforra  and  mode  shapes  must  be  normalized.  The  planform 
must  have  a  semi-span  of  1.0  and  the  mode  shapes  must  correspond  to  the 
normalized  planform.*** 


.FALSE. 


Do  not  print  generalized  aerodynamic 
coefficients 


WTGNAF 


Logical 


.TRUE. 


Write  generalize!  air  forces  on  tape. 


.FALSE. 


Mo  generalized  air  forces  are  to  be 
written  on  tape. 


WTBL 


Logical 


.TRUE. 


Write  box  lifts  on  tape. 


•FALSE. 


Box  lifts  are  not  to  be  written  or.  taj  . 
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Card  Set  C  -  Continued 


The  following  parameters  may  be  included  on  Card  Set  C:  DPPCPR,  GEOCPR, 
MODCPR ,  AICCPR,  NWSCPR,  SMCPR,  GAFCPR.  They  are  intended  for  checkout 
only,  ar.d  are  described  in  Part  II,  Section  II,  Subsection  11,  Common 
Block/ CHECKPR/. 


Card  Set  D  -  Tape  Parameters,  all  are  Integers. 


(  $CARDD 


Format  -(NAMELIST) 


Variable 


Default 


Meaning 


'OAIC| 
or  > 

^Aicj 


0 


Logical  tape  number  of 


AIC  tape .  An 


"old"  tape  is  one  which  was  generated  as  a  new 
tape  by  a  previous  execution  of  the  program, 
value  of  "zero"  means  no  tape  is  to  be  us#ed  or 
saved. 


A 


'OSAIC') 
<  or  > 
.NSAICj 


0 


Logical  tape  number  of 


Spatial  AIC  tape. 


An  "old"  tape  is  one  which  was  generated  as  a  new 
tape  by  a  previous  execution  of  the  program.  A 
value  of  "zero"  means  no  tape  is  to  be  used  or 
saved. 


INTAPE 


NOUTP 


0 


1 


Logical  tape  number  of  tape  containing  geometry  and 
mode  shape  information.  If  this  v<Lue  is  0  or  5» 
card  input  will  be  used,  requiring  cards  N  to  0. 
Tape  is  always  rewound  before  use. 

Logical  tape  number  for  output  requested  under 

WT -  options.  This  parameter  is  net  needed 

if  no  WT -  option  is  specified.  Tape  is  re¬ 

wound  first  cycle,  used  from  current  position 
in  subsequent  cycles. 


NOTE:  The  above  four  tape  numbers  must  be  independent, 
.  if  specified. 


INFSP 


0 


Initial  file  spacing  on  INTAPE. 


IOUFSP 


0 


Initial  file  spacing  on  NOUTP. 
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-4 

i 


i 


-  Reduced  Frequencies,  Real 


Format  -  (NAMELIST) 


XKS(l)  = 
or, 

XKl(l)  = 


XKS(2)  = 
XK1(2)  = 


$ 


Array  of  reduced  frequency  values,  based  on 


f  semi span ,s 
or 

box  length 


for  which  airforces  will  be  calculated.  Limit  20.  If  no 
reduced  frequencies  fire  specified,  the  program  will  determine 
the  geometry  und  then  terminate. 


Geometry  Definition 

(None  of  these  cards,  F  through  L,  are  needed  if  PRVGEOM  has  "been  specified.) 


Card  F  Frame  of  Reference  Definition. 

All  planform  description  and  mode  definition  is  based  on 
coordinate  systems  local  to  each  surface.  The  values  on 
this  card  relate  the  local  systems  to  a  global  system. 
(See  Figures  1  &  2) 


/  SCARDF 


Format  -  (NAMELIST) 

$ 


WLAX  =  X  location  of  the  wing  local  axis  system  origin  relative  to 
the  global  coordinate  system.  Positive  aft.  Real, default  0. 


WLAZ  =  Z  location  of  the  wing  local  axis  system  origin  relative  to  the 
global  coordinate  system.  Positive  upward.  Real, default  0. 


PSIW  = 


Wing  dihedral  angle,  degrees,  positive  upward  for  the  right 
hand  side.  Real, default  0. 


TLAX 

TLAZ 

PSIT 


Same  as  above,  for  the  tail  (second  surface)  if  needed. 


Card  G  Box  pattern  size  and  location  Format  -  (NAMELIST) 


r 


$CARDG  N CURDS  = 


XCENTR 

or 

XEDGE 


$ 


NCHRDS  =  Number  of  boxes  in  the  spanwise 


direction  on  the  wing.  Integer, 
no  default.  See  subsection  3  of  Section  III,  Program  Limits  and 
Restrictions. 


XCENTR 
or 
i  XEDGE 


X  location  of  a  box 


center 

edge 


in  the  local  coordinates  of  the  wing, 


from  which  X  is  determined.  XCENTR  or  XEDGE  may  be  specified  anywheri 
on  or  near  tHe  wing.  If  XEDGE  is  specified,  the  program  computes 
a  corresponding  XCENTR  =  XEDGE  +  b^/2 .  Then  the  program  adds  or 
subtracts  increments  of  b^  until  the  first  on-planform  wing  location 
is  reached.  The  resulting  value  is  XCv  the  location  of  row  1  of  the 
box  pattern.  The  same  grid  is  extended  to  the  tail  as  needed. 

Real,  no  default. 
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Card  G  -  Continued 


ICHORD  =  Array  of  chord  numbers  for  sample  of  off-planform  upvash,  sidewash 
or  longitudinal  wash.  Integer,  default  is  no  sample.  ISiOPLW 
(Card  C)  gives  the  maximum  number  expected. 


IBOXF  \ 
IBOXL) 


Array  of  ^ast  pox  samP^e<*  on  eac^  ch°**<3,  Integer,  default 

IBOXF  *  IBOXL  *  last  box  of  the  planfora.  Measured  along  coordinate 
nc .  0  <  IBOXF  (I )  *  IB0XL(  I )  *  kO 


ZLOC  =  Array  of  Z  -  locations  of  the  sample  chords  in  the  global  coordinate 
system.  Real,  default  0.  (The  yertical  separation  between  wing  and 
sampling  chord  must  be  non-zero, } 

NOTE:  ICHORD,  IBOXF,  IBOXL,  ZLOC  will  be  ignored  unless  ISMPLW  on  Card  C 

is>0,  and  NSURF  e  1.  ICHORD,  IBOXF,  IBOXL, ZLOC  are  arrays,  each 
one  requiring  ISMPLW  numbers. 


FIGURE  23  Definition  of  Sampling  Locations 

(Assuming  =  0,  ZLOC(I)  not  shown) 
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Modal  input  is  not  required  if  no  reduced  frequencies  were  listed  on 
Card  E  or  if  PRVMODE  was  specified  on  Card  C.  Otherwise,  Card  M  must 
be  included,  and  the  remainder  are  required  if  INTAPE  =  0  or  5. 


Card  M  -  Number  of  Modes  and  Thickness  slope  distributions 
This  card  is  read  once  per  cycle  of  the  program. 

- - 1 

(  $CARDM  NMODES  =  _  ...  $ 


LISL) _ 


NMODES  =  Total  number  of  modes  to  be  input.  A  "mode"  as  used  here  means 

the  combined  distortion  of  the  wing  and  tail.  Integer,  no  default. 

NTSLOP  =  Number  of  thickness  slope  sets  to  be  input  on  Card  Sets  P  and  Q. 

If  omitted  the  program  takes  a  default  of  0  and  generates  an 
array  of  zero  thickness  slope  everywhere.  Integer. 


There  are  three  options  for  the  form  in  which  mode  shapes  may  be  input  from 
cards.  The  form  indicated  by  MTYPEW  will  be  repeated  NMODES  times,  getting 
all  modes  for  the  wing.  Then,  if  NSURF  =  2,  cards  N  and  0  are  repeated 
NMODES  times  again,  this  time  in  the  form  indicated  by  MTYPET,  to  obtain 
the  corresponding  mode  shapes  for  the  tail.  If  there  is  no  motion  of  the 
wing  (or  tail)  for  a  particular  mode,  the  corresponding  modal  vslues  must 
be  input  as  zero. 

fMTYPEW) 

Option  1:  Polynomial  Input  -  Cards  Nl,  01  i^j^pyPETj  ^ '  RePeated  NMODES  times. 


Card  Nl 


Format-  (15) 


IDEG  =  Degree  of  polynomial,  «  5*  This  controls  the  number  of  coefficients 
on  card  Set  01. 


Card  Set  01 


Format  -  (7E10.0) 


As  many  coefficients  as  are  needed  for  the  polynomial  of  degree  IDEG: 

O  o  o 

fj(X’Y)  =  aoo  +  a10X  +  a01Y  +  a20X  +  allXY  +  a02Y  +  a30X  +  ‘',ao,IDEGY 
where  point  (X,Y)  is  in  the  (or  XT,  YT)  coordinate  system,  and 


f  (X,Y)  is  the  modal  deflection. 


Card  L  -  Tail  trailing  edge  definition  points.  (Omitted  if  NSURF  =  1 
on  Card  C). 


f 


XTTE{1)' 


"ytteTiT  XTTE( 2 )'  ytte(2)v 


Format  -  ( bK10 . 0 


YTTE( 


t)\ 

i) 


XTTE(i)l  -  X  and  Y  coordinates  of  the  trailing  edge  definition  points  for 


the  tail.  There  are  NTTE  pairs  of  coordinates,  measured  in  the 


tail  local  coordinate  system,  (X^,Yfj,). 
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Modal  Input 


Modal  input  is  not  required  if  no  reduced  frequencies  were  listed  on 
Card  E  or  if  PRVMODE  was  specified  on  Card  C.  Otherwise,  Card  M  must 
be  included,  and  the  remainder  are  required  if  INTAPE  =  0  or  5. 


Card  M  -  Number  of  Modes  and  Thickness  slope  distributions 
This  card  is  read  once  per  cycle  of  the  program. 


r 


$CARDM 


NMODES  = 


Format  -  ( NAMELIST ) 


NMODES  =  Total  number  of  modes  to  be  input.  A  "mode"  as  used  here  means 

the  combined  distortion  of  the  wing  and  tail.  Integer,  no  default. 


NTS LOP  =  Number  of  thickness  slope  sets  to  be  input  on  Card  Sets  P  and  Q. 
If  omitted  the  program  takes  a  default  of  0  and  generates  an 
array  of  zero  thickness  slope  everywhere.  Integer. 


There  are  three  options  for  the  form  in  which  mode  shapes  may  be  input  from 
cards.  The  form  indicated  by  MTYPEW  will  be  repeated  NMODES  times,  getting 
all  modes  for  the  wing.  Then,  if  NSURF  =  2,  cards  N  and  0  are  repeated 
NMODES  times  again,  this  time  in  the  form  indicated  by  MTYPET,  to  obtain 
the  corresponding  mode  shapes  for  the  tail.  If  there  is  no  motion  of  the 
wing  (or  tail)  for  a  particular  mode,  the  corresponding  modal  values  must 
be  input  as  zero. 

fMTYPEWi 

Option  1:  Polynomial  Input  -  Cards  Nl,  01  1*  Repeated  NMODES  times. 

Card  N1  Format-  (15) 


5 


IDEG 


IDEG  =  Degree  of  polynomial,  <  5.  This  controls  the  number  of  coefficients 
on  card  Set  01. 

Card  Set  01  Format  -  (7E10.0) 


r 


00 


10 


01 


*20 


‘11 


02 


30 


r 


21 


12 


03 


.etc . 


As  many  coefficients  as  are  needed  for  the  polynomial  of  degree  IDEG: 

=  %  +  a10X  +  a01Y  +  a20x2  +  allXY  *  V*2  +  V3  +  ••■ao,IDE0yIDEG 
where  point  (X,Y)  is  in  the  X^.Y^  (or  X^,,  Y^,)  coordinate  system,  and 

fj (X,Y )  is  the  modal  deflection. 
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Modal  Input  -  Continued 


Option  2: 


Interpolation  -  Cards  N2, 


02  if 


MTYPEW 

MTYPET 


*  2. 


Repeated  MODES  times. 


Card  H2 


Format  -  (215) 


(  IDEG  NPTS 

IDEG  =  Degree  of  polynomial  to  be  least-squares  fit  to  the  surface,  4  5. 
NPTS  =  Number  of  points  at  which  deflections  are  known,  *100. 


Card  Set  02  Locations  and  Deflections,  NPTS/2  cards  Format  -  (6E10.0) 

XJT  '  yy"  '  ZZ "  '  XXL  1  YY_  1  zzT" 

I  -*•  1  -L  2  c,  2 

(  XX^  YY3  1  ZZ3  '  . etc. 

YY 

Y  I  =  The  X  and  Y  location  of  the  ith  point  in  the  planforro  local  coordinate 
system,  Xy,  Yy  (or  Xy,  Yy) 

ZZ^  =  The  deflection  at  point  i,  perpendicular  to  the  surface. 


Option  3 
N3  and  03  if 


Deflections  and  Slopes  at  Box  Centers  -  Card  Sets 
MTYPEW 


MTYPET | 


=  3.  Repeated  llMODES  times. 


Card  Set  H3  -  Deflections  array  for  the  planform. 


_  _  _ Format  -  (7E10.0) _ 

^defl(i)  r  DEFL(2T  ~ 

DEFL(I)  -  Deflection  values  must  be  given  for  the  control  points  of  nl ’ 
the  Mach  boxes  on  the  right  half  surface.  These  values  start 
at  the  leading  edge  root  box  and  proceed  sequentially  fore  to 
aft  on  each  chord,  one  chord  at  a  time  until  the  tip  chord  has 
been  completed.  Each  chord  starts  on  a  new  card. 

Card  Set  03  -  Slopes  array  for  the  planform. 

_ _ Format  -  (7E10.0) _ 

(  slopeTTT  slope  (2  y  rn 

SLOPE(I)  -  The  slope  array  is  given  for  the  same  points  and  ir.  the  same 
order  as  the  deflections  array. 
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Thickness  Slopes 

The  optional  thickness  slope  values  follow  all  modal  input.  The  integer 
variable  NTSLOP  on  Card  M  determines  the  presence  or  omission  of  Card 
Sets  P  and  Q.  If  NTSLOP  *=  0  these  card  sets  are  omitted,  otherwise  there 
will  be  NTSLOP  sets  of  then.  The  thickness  slope  arrays  have  the  same 
order  of  input  as  the  Mode  Shapes,  Option  3.  All  the  card  sets  for  the 
wing  (card  sets  P)  come  first,  followed  by  all  the  card  sets  for  the  tail 
(card  sets  Q). 

Card  Set  P  -  Thickness  slopes  at  Box  Centers  for  the  tail. 

Format  -  (7E10.0) _ 

(  TSLFW(l)'  TSLFW(2) 

iz^ 

TSLFW  =  Array  of  — — — —  values  of  the  wing  at  box  centers,  input  ir.  the 

Sx 

same  order  as  deflections  on  Card  Set  N3.  See  Figure  18,  Section 
II,  subsection  7. 

Card  Set  Q  -  Thickness  slopes  at  Box  Centers  for  the  tail. 

_ _ *  _ Format  -  (7E10.0) _ 

(  TSLFT(l)  TSLFT(2) 


Termination  Control  Card: 


After  completion  of  all  calculation  defined  by  the  previous  data,  the 
following  card  is  read.  If  anothe?  cycle  is  indicated,  the  program  will 
expect  a  new  set  of  data,  starting  with  Card  A  (Title),  Otherwise  the 
program  will  terminate  as  indicated: 


r 


NAME 

NAME 


LI  > 


L2  = 


_ fl _ _ Format  (A6.to.2I10) 

’ -  —y  — 

File  name  of  next  overlay  to  be  executed,  if  L1>0. 

Ignored  if  L1«0. 

0,  Primary  level  number  of  next  overlay  to  be  executed. 

0,  Call  EXIT,  returning  to  control  cards. 

1,  Recycle  to  Card  A 

2,  Return  to  calling  overlay. 

Secondary  level  overlay  number  (ignored  if  L1<0), 
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DIAGNOSTICS 


The  program  prints  warning  and  fatal  diagnostics  when  error  conditions 
are  detected.  A  warning  diagnostic  is  used  to  attract  attention  to  a 
condition  which  is  not  fatal  to  the  basic  program  execution,  but  which 
may  be  symptomatic  of  a  data  error.  It  will  usually  start  with  "*** 
WARNING — After  a  warning  diagnostic  execution  continues.  A  fatal 
diagnostic  is  followed  by  a  call  to  the  subroutine  FLUSH,  which  may  be 
written  to  interface  with  the  error  exit  capability  of  the  user's 
operating  system.  A  fatal  diagnostic  starts  with  n***ERR0R - ." 


a.  General  Disk  Diagnostic 


***ERR0R  -  WHILE 


(READINGl 
WRITING! 


GEOMETRY 
THICKNESS  SLOPE 
MODES 
PLANAR  AIC 
SPATIAL  AIC 
VELOCITY  POTENTIAL/ 


>  SCRATCH  FILE  nnnnn*** 


ERROR  CODE  nnnn  WHILE 


READING 

WRITING 


THE  FOLLOWING  MATRIX 


MATRIX  ID 

PARAMETERS  n.nnnEnn  n.nnnEnn 

(INTEGER)  nn  nn 

FILE  SPACING  =  nnn 
MATRIX  TYPE  = 


n . nnnEnn  n . nnnEnn 

nn  nn 

MATRIX  SPACING  =  nnn 
DIMENSIONED  (nn  x  nn) 


The  error  code  indicates  the  type  of  error  encountered. 

1  Matrix  spacing  is  negative 

2  File  spacing  is  negative 
Matrix  dimensions  are  illegal 

5  Row  dimension  of  matrix  exceeds  the  row  dimension  of 
BUFFER  ARRAY. 

1500+1  Encountered  EOF  after  matrix  I  while  skipping  matrices. 


b.  Diagnostics  from  the  Data  Preprocessor 


Four  warning  diagnostics  may  occur: 


***WARNING  -  MACH  NUMBER 


(  LESS  THAN  1.2 
1  GREATER  THAN  3.0 


IS  BEING  USED. 


The  Mach  range  for  which  reasonable  answers  can  be  expect*;') 
is  1.2  ^  3.0. 


***WARNING  -  NO  K1  VALUES  GIVEN.  PROGRAM  WILL  TERMINATE  AFTER 
GEOMETRY  SECTION*** 


No  values  were  specified  on  Card  K,  non  fatal. 

***WARNING  -  ORDER  FOR  VELOCITY  POTENTIAL  SMOOTHING  TOO  LA' V  . 
IT  HAS  BEEN  REDUCED  TO  10*** 

Too  large  a  value  was  read. 
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"•"WARNING  —  NO  OUTPUT  TAPE  WAS  REQUESTED  FOR 


GENERALIZED  FORCES 
BOX  LIFTS 


The  parameters  WTGNAF  or  WTBL  were  specified,  but  NOUTP  =  0 
on  Card  D. 

One  fatal  diagnostic  may  occur: 

*«*MACH  NUMBER  { .  0  }  CAMOT  BE  USED*** 

This  would  normally  be  a  data  error.  The  theory  used  in  the 
program  is  not  valid  for  Mach  numbers  outside  this  range. 

Diagnostics  from  the  Geometry  Section 

"""WARNING  -  XEDGE  AND  XCENTR  WERE  BOTH  SPECIFIED.  XEDGE  WILL 
BE  IGNORED""" 

The  grid  location  may  be  specified  at  a  box  edge  or  a  box 
center,  Card  G. 

"""WARNING  -  SAMPLING  OF  UPWASHES  CANNOT  BE  DONE  IF  A  TAIL  HAS 
BEEN  DEFINED.  ISMPLW  =  nn  WILL  BE  IGNORED""" 

Self  explanatory. 

"""WARNING  -  SAMPLE  WASH  SPECIFICATION  SET  nn  IS  IN  ERROR.  ONLY 
THE  PRECEDING  WILL  BE  CALCULATED""* 

On  card  G,  ICHORD(nn)  »>NCHRDS,  or  I  CHORD  (nn)  <0,  or 
IBOXF(nn)  ^IBOXL(nn),  or  IBOXF(nn)  -cl,  or  IBOXL(nn)  >40. 

"""ERROR  -  PARAMETER  NCHRDS  WAS  NOT  SPECIFIED.  IT  MUST  ALWAYS 
BE  GIVEN**" 

Card  G,  required  data. 

"""ERROR  -  EITHER  XEDGE  OR  XCENTR  MUST  BE  SPECIFIED*** 

Card  G,  required  data. 

"""ERROR  -  THE  TAIL  AND  THE  WING,  OR  THEIR  DIAPHRAGMS,  CROSS- 
ABOVE  TO  BELOW*** 

If  any  portion  of  the  tail  is  above  the  wing  plane,  all  of  it 
and  its  diaphragm  must  be. 
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***ERROR  -  XCENTR  NOT  WITHIN  50  BOX  LENGTHS  (B1  =  n.nnnnnnnnEnn) 
OF  THE  WING  L.E.  (n.nnnnnnnnEnn)*** 


It  is  assumed  that  the  user  can  place  his  box  grid  fairly 
accurately,  and  any  value  for  XCENTR  that  is  wildly  off 
must  be  a  data  error.  If  XEDGE  was  specified,  the  program 
calculates  XCENTR  =  XEDGE  +  Bl/2. 


***ERR0R  - 


PSIW 

PSIT 

NCHRDS 

NWLE 

NWTE 

NTLE 

NTTW 


IS  OUTSIDE  ALLOWED  RANGE*** 


See  Section  III. 3,  Program  Limits  and  Restrictions. 


WING  L.E. 

***ERROR  -  ]WING  T*E- 
|  TAIL  L.E. 

TAIL  T.E. 


DEFINITION  POINTS  ERROR  nn,  A  COMBINATION 


( 


OF  1,  NONMONOTONIC  Y-VALUES 
k,  FIRST  Y-VALUE  NONZERO 


2,  NONMONOTONIC  X-VALUES 
8,  TIP  T.E.  Y-VALUE  DISAGREES 
WITH  TIP  L.E.  VALUE 


The  error  code  nn  may  equal  the  sum  of  two  or  more  conditions, 
if  more  than  one  error  was  detected.  The  input  values  being 
checked  are  from  Cards  I,  J,  K  and  L. 


***ERR0R  PROCESSING 


WING 

TAIL 


GEOMETRY, 


LEADING  } 
TRAILING  f 


EDGE*** 


This  general  header  may  occur  during  calculation  of  box  codes, 
and  precedes  one  of  the  following  2  descriptors: 


SECTION  nn  IS  BEYOND  THOSE  DEFINED 


Should  only  occur  with  a  machine  failure  or  coding  error. 
SECTION  nn  OF  THE  EDGE  DOUBLES  BACK  TOWARD  THE  CENTER  LINE 


Nonmonotonic ally  increasing  Y-values  have  been  encounter 

WING 


***ERR0R  -  SECTION  nn  OF  THE  TRAILING  EDGE  OF  THE 


CHORD  J  TO  GO  TO  ROW  i ,  WHICH  EXCEEDS  THE  LIMIT*** 


TAIL 


CAUSES 


Too  many  rows  are  being  requested  as  a  function  of  Mach  R>unb< 
and  number  of  chords. 
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***ERROR  -  TOO  MANY  CHORDS  FOR  BOX  CODE  ARRAY*** 


This  occurs  when  additional  tip  diaphragm  chords  are  being 
determined.  The  resulting  box  pattern  exceeds  the  program 
limits,  subsection  3. 

The  above  three  diagnostics  are  followed  by  the  box  code  pattern 
currently  being  developed,  to  aid  in  determining  the  necessary 
data  changes.  The  program  then  terminates. 

d.  Diagnostics  from  the  Modes  Section 


***WARNING  -  PREVIOUS  MODE  SHAPES  HAVE  BEEN  SPECIFIED,  BUT  GEOMETRY 
HAS  CHANGED.*** 

The  mode  shapes  requested  were  input  as  coefficients  on  the 
previous  cycle.  The  planform  has  changed  and  the  mode  shapes 
must  be  calculated  at  the  new  box  centers. 

***ERROR  -  PREVIOUS  MODE  SHAPES  HAVE  BEEN  SPECIFIED,  BUT  THE 
GEOMETRY  HAS  CHANGED.  PREVIOUS  MODE  SHAPES  WERE  AT  BOX 
CENTERS.*** 

The  mode  shapes  requested  were  input  at  box  centers  on  a 
previous  cycle.  The  planform  has  changed  and  the  mode  shapes 
are  no  longer  correct  for  the  box  centers. 

e.  Diagnostics  from  the  AIC  Section 

***ERROR  -  THE  SIZE  OF  THE  AIC  ARRAY  FOR  THIS  PLANFORM  IS  NNNN , 

THE  MAXIMUM  SIZE  ALLOWED  IS  MMMM. *** 

The  planform  for  this  run  is  such  that  it  creates  a  requirement 
for  an  AIC  array  of  size  NNNN.  The  largest  AIC  array  that  can 
be  used  is  MMMM.  The  planform  should  be  checked  to  determine 
if  size  restrictions  have  been  violated. 

***ERR0R  -  THE  ARGUMENT  FOR  A  BESSEL  FUNCTION  IS  OUT  OF  RANGE.*** 

ARGUMENT  =  n.nnnnnn 
LOWER  LIMIT  =  n.nnnnnn 
UPPER  LIMIT  =  n.nnnnnn 

Coding  error  or  machine  failure  during  integration  of  the  AIC 
equations. 

f .  Diagnostics  from  the  Normal  Wash  and  Velocity  Potential  Section 

***WARNING  -  NO  PLANFORM  CONTROL  POINT  FOUND  FOR  SUBDIVIDED  BOX 
DURING  VELOCITY  POTENTIAL  CALCULATIONS*** 

SUBDIVIDED  BOX  (i.J),  CONTROL  POINT  (ii,jj) 
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.  *7 


A  high-swept  leading  or  trailing  edge  segment  may  yield 
subdivided  on-planform  boxes  some  distance  from  the  nearest 
(chordwise)  control  point.  Since  the  substantial  derivative 
for  any  planform  subdivided  box  is  the  chordwise  extrapolated 
value  of  the  nearest  planform  control  point,  no  extrapolation 
beyond  two  box  lengths  will  be  attempted.  The  subdivided  box 
in  question  will  be  ignored. 

***ERR0R  -  THE  TIP  BOX  PATTERN  DOES  HOT  ALLOW  TRAILING  VELOCITY 
POTENTIALS  TO  BE  COMPUTED*** 

The  tip  chord  must  have  two  or  more  on-planform  control 
points,  or  the  inboard  chords  must  have  on-planform  control 
points  as  indicated  in  Figure  2k. 


FIGURE  2k  BOXES  NEEDED  FOR  TRAILING  EDGE  VELOCITY  POTENTIAL 
CALCULATIONS 

***ERROR  -  FAILURE  IN  STORING  (  cOTTRO^POINT )  N0RMAL“WASHES*** 
SUBDIVIDED  BOX  (i.j),  CONTROL  POINT  (ii,jj) 

This  is  symptomatic  of  a  machine  failure  or  a  coding  error. 

Diagnostics  for  the  Velocity  Potential  Smoothing  Section 

***ERROR  -  NO  TIP  TRAILING  EDGE  VELOCITY  POTENTIAL  CAN  BE  COMPUTED.*** 

Same  diagnostic  as  in  Normal  Wash  and  Velocity  Potential  section. 
The  tip  chord  must  have  two  or  more  on-planform  control  points, 
or  the  inboard  chords  must  have  on-planform  control  joints  as 
indicated  in  Figure  2k, 

Diagnostics  for  the  Forces  Section 

***P0INTER  ARRAY  EXCEEDED  FOR  BOX  (I,J)**» 

This  is  symtomatic  of  a  machine  failure  or  a  coding  error. 
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7.  OUTPUT  DESCRIPTION 


a.  Printed  Output 

The  first  page  of  the  printed  output  consists  of  a  program  title 
block,  the  user's  run  title,  and  a  list  of  the  parameters  and 
options  requested.  This  includes  Mach  number,  symmetric  or  anti¬ 
symmetric  ,  basic  or  subdivided  analysis ,  AIC  integration  accuracy 
indicator,  and  all  print  options.  All  tape  variables  (input,  out¬ 
put,  and  AIC  tapes)  are  indicated,  and  the  list  of  reduced  fre¬ 
quencies  is  printed. 

The  geometry  section  prints  all  card  data  pertaining  to  geometric 
definitions.  Any  values  not  specified  in  the  input  data  are 
printed  starting  with  an  asterisk  (*),  and  the  defaults  are  used. 

If  printing  of  the  box  pattern(s)  is  requested  (PRBOX),  these 
follow.  There  is  one  pattern  printed  for  the  control  points  of 
each  surface,  and  if  subdivision  is  requested  there  is  a  sub¬ 
divided  box  code  pattern  printed  for  each  surface. 

The  modes  section  prints  (optionally)  the  modal  data.  This  may 
be  polynomial  coefficients  (PRCOEF)  and/or  deflections  and  slopes 
at  box  centers  (PRMODS). 

The  AIC  arrays  to  be  used  are  printed  (optionally)  as  needed.  For 
each  reduced  frequency,  any  spatial  AIC  arrays  are  printed  first 
(PRSAIC),  followed  by  the  planar  AIC  array  (PRPAIC),  for  which  W^jj, 
and  V-jjQ  are  0.0.  AIC  arrays  needed  for  sampling  are  printed  as 
spatial  arrays. 

The  normal-washes  may  be  printed  under  option  PRNW.  They  are 
printed  in  the  order:  N 

N  Lin-  If  option  PRVP  is  _  _  _  = 

ferences  are  next  printed,  for  the  wing  followed  by  A$  for 
the  tail.  If  sampling  has  been  requested  (ISjpLW^O),  the  upwash,  jj, 
the  side  wash  jj  ,  and  the  longitudinal  wash,  y  ,  may  be  printed  for 
all  boxes  defined  on  each  sampling  chord  requested,  under  options, 
PRUW,  PRSW,  and  PRLW.  All  of  the  printing  requested  in  this  section 
is  repeated  for  each  mode  shape. 

If  either  of  the  velocity  potential  smoothing  options  was  specified, 
a  second  set  of  (yf  arrays  for  all  mode  shapes  will  be  printed  by 
the  smoothing  routine  under  the  control  of  print  option  PRVP. 

The  air  forces  section  piints  (optionally)  box_lifts  Lj’m,  pressure 
difference  coefficients  ACp^m,  section  lifts  3?j,  generalized  air 
forces  and  generalized  aerodynamic  coefficients  Q'  and  QM. 

The  options  controlling  this  are  PRBL,  PRDCP,  PRSL,  PRGNAF,  and 
PRGNAC.  If  a  velocity  potential  smoothing  option  was  specified,  two 
sets  of  box  lifts,  pressures,  etc.,  are  computed  and  printed.  The 
first  set  corresponds  to  the  smoothed  A*  arrays,  and  the  second  set 
corresponds  to  the  unsmoothed  A$'s*  This  section  terminates  the 
printing  for  each  reduced  frequency.  Printing  for  subsequent  reduced 
frequencies  starts  again  at  the  optional  printing  of  AIC  arrays. 


RUW’ 

spec 


Nrt„,  and  if  a  tail  is  defined,  N 
iried,  the  velocity  potential  dif- 
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Tape  Output 


Tape  Output  is  optionally  available  for  box  lifts  Ln,m  and  generalized 
air  forces  Q^.  The  planar  AIC  arrays,  Og,£,o,and  the  spatial  AIC 

arrays,  CvjlX,  Wvpl,  and  V\jjj^,  may  be  saved  for  reuse  by  the  program 
during  subsequent  execution.  Tape  maps  for  all  three  are  in  Part  II. 
TABLE  1  ORDER  OF  ALL  OPTIONAL  OUTPUT 


Option  List 

Reduced  Frequency  List 

Geometry  and  Box  Patterns 

Mode  Shapes 

Spatial  AICS 

Planar  AIC 

Subdivided  AIC 

1 

Wing  Upper  Normal  Wash 

Wing  Lower  Normal  Wash 

Tail  Upper  Normal  Wash 

Repeated 

for 

Tail  Lower  Normal  Wash 

each  mode 

Wing  Velocity  Potential  Difference 
Tail  Velocity  Potential  Difference 

shape 

i 

Repeated  for 
►  Each  Reduced 

Smoothed  Wing  Velocity  Potential  Differences)  Repeated  for 
Smoothed  Tail  Velocity  Potential  Differences/  each  mode  shape 

Frequency 

Wing  Box  Lifts 

Tail  Box  Lifts 

Repeated 

Wing  Pressure  Differences 

for  each 

Tail  Pressure  Differences 

Section  Lifts 

mode  shape 

Total  Lift  „ 

<y 

Sectional  Generalized  Forces  (c  C  ) 

m 

JZ 

u 

o 

Generalized  Air  Forces 

o 

B 

C/3 

Generalized  Aerodynamic  Coefficients 

Wing  Box  Lifts 

Tail  Box  Lifts 

Repeated 

Wing  Pressure  Differences  | 

for  each 

Tail  Pressure  Differences 

Section  Lifts 

mode  3hape 

0) 

Total  Lift  2 

rC 

u 

Sectional  Generalized  Forces  (c  C  ) 

m 

o 

o 

B 

Generalized  Air  Forces 

(0 

e 

Generalized  Aerodynamic  Coefficients 

i 

t 


i 
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SECTION  IV 


APPLICATIONS 


1 .  INTRODUCTION 

The  computer  program  described  in  the  previous  sections  has  been 
developed  to  evaluate  the  unsteady  airforces  existing  on  wing-tail 
lifting  surfaces  in  supersonic  flow.  This  portion  of  the  report 
describes  some  of  the  operational  aspects  of  the  program  and  presents 
an  application  of  the  numerical  procedures  used  in  a  flutter  analysis 
of  wing-tail  models  previously  tested  at  Cornell  Aeronautical 
Laboratory. 

Pressure  distributions,  perturbation  velocities,  and  generalized 
forces  are  presented  for  various  planform  configurations  identified 
as  AGARD  configurations. 


2.  PROGRAM  LIMITATIONS 

Configurations  that  fit  the  program  analysis  capabilities  are  shown 
in  figure  25,  figure  26,  and  figure  27.  The  planforms  must  be  de¬ 
fined  in  terms  of  straight  line  segments  and  may  take  on  various 
shapes  as  indicated  in  figure  25.  The  wing  and  tail  may  have 
dihedral,  vertical  separation  and  horizontal  separations.  Aerodynamic 
interaction  of  the  tail  affecting  the  wing  may  be  evaluated  only  for 
the  coplanar  case.  The  tail  may  physically  overlap  the  wing  for  the 
spatial  cases  however  the  tail  aerodynamic  "Mach  cone"  must  not  in¬ 
tersect  the  wing  as  shown  in  figures  28  and  29.  A  complete  set  of  other  pro¬ 
gram  limitations  are  presented  in  section  III,  Computer  Program  Usage. 

3.  PRELIMINARY  COMMENTS  ON  PROGRAM  CHARACTERISTICS 

In  applying  the  numerical  technique  known  as  the  "Mach  Box"  method 
to  the  solution  of  supersonic  flow  problems  one  is  bound  by  the  res¬ 
trictions  of  the  method  and  must  accept  the  errors  incurred  by  the 
built-in  approximations.  Probably  the  most  critical  approximation 
made  within  the  present  method  is  the  assumption  that  the  normal  washes 
are  of  constant  value  over  the  individual  box  areas.  This  approxima¬ 
tion  may  apply  without  too  much  error  within  regions  on  the  lifting 
surface  where  small  normal  wash  variations  exist  and  also  provided 
that  the  box  sizes  are  relatively  small.  However,  relatively  large 
errors  may  result  in  predicting  the  normal  wash  distribution  over  the 
"off-wing"  diaphragm  region  where  large  variations  in  the  normal  wash 
exist  that  are  inversely  proportional  to  the  square  root  of  the  dis¬ 
tance  from  the  planform  edge. 


91 


In  order  to  allow  a  variation  in  the  basic  box  size  to  reduce  the 
errors  in  the  downwash  calculation,  the  program  has  been  devised  such 
that  two  basic  sizes  of  boxes  may  be  used  within  the  solution  process. 
The  basic  program  uses  a  uniform  size  box  over  the  entire  wing  and 
"off -wing"  diaphragm  regions.  The  application  of  the  optional  data 
item  SUBDV  =  .TRUE,  allows  a  three  to  one  subdivision  of  the  basic 
box  lengths  over  an  "effective  area"  of  the  wing  and  "off-wing"  dia¬ 
phragm  region.  Figure  31  shows  the  distribution  of  uniform  size  boxes 
over  the  wing  and  the  "off -wing"  diaphragm  used  in  the  basic  analysis 
case.  Figure  32  displays  the  "effective  area"  concept  used  in  the 
subdivision  case.  The  streamwise  width  of  the  "effective-area"  may 
be  maintained  at  a  constant  value  or  may  be  variable  upon  the  users 
request.  The  "effective  area"  slides  from  row-to-row  in  a  streamwise 
direction  and  is  used  to  obtain  a  more  accurate  definition  of  the 
velocity  potentials  at  the  box  centers  of  the  basic  case. 

It  should  be  noted  that  once  the  box  pattern  has  been  selected  to 
approximate  the  original  planform  shape  then  the  original  shape  no 
longer  may  be  used  in  the  analysis.  There  is  one  modification  of  this 
condition  in  that  in  using  the  subdivision  option  the  subdivided  boxes 
are  used  to  refine  the  velocity  potentials  at  the  box  center  of  the 
basic  box  pattern  and  do  not  enter  into  the  box  load  or  ?enex-alized 
force  calculations.  However,  the  leading  edge  and  trailing  edge  box 
loads  may  include  the  effect  of  the  additional  subdivided  boxes  lo¬ 
cated  at  the  leading  and  trailing  edges  provided  that  the  program 
option  PLYWOOD  is  not  used.  The  option,  PLYWOOD  =  .TRUE. ,  forces  all 
boxes  to  have  the  same  area  as  ^iven  in  the  basic  analysis  case  by 
setting  the  area  ratio  factor  a*  =1.0  for  all  box  areas.  If  the 
program  option  PLYWOOD  is  not  used,  then  the  program  will  calculate 
the  on*m  values  used  to  modify  the  box  areas  contained  within  the 
leading  and  trailing  edge  regions  for  use  within  the  loads  and  general¬ 
ized  force  calculations.  Figure  33  displays  the  method  used  to  eval¬ 
uate  the  an,m  factors  used  in  the  subdivision  case.  It  should  be  noted 
that  the  option  PLYWOOD  =  .TRUE.  (an,m  =  1,0)  is  to  be  applied  for  all 
cases  except  when  subdivision  is  being  used,  otherwise  the  original 
mathematical  definition  of  the  analysis  planform  will  be  violated  and 
will  produce  inconsistent  results. 

The  purpose  of  subdividing  the  box  sizes  into  smaller  boxes  is  to 
smooth  out  the  streamwise  variations  in  the  velocity  potentials  since 
it  is  estaDlished  that  the  chordwise  loadings  are  a  function  of  the 
velocity  potentials  and  the  first  derivatives  of  the  velocity  poten¬ 
tials.  Small  variations  in  the  velocity  potentials  will  cause  large 
variations  in  the  chordwise  loadings  and  may  appear  to  be  unrealistic 
from  a  physical  standpoint.  The  variations  in  the  chordwise  defini¬ 
tion  of  the  velocity  potentials  are  caused  by  the  irregularities  of 
the  planform  defined  by  the  Mach  box  pattern  and  also  is  affected  by 
the  amount  of  diaphragm  existing  between  the  foremost  Mach  line  and 
the  leading  edge  of  the  planform.  An  example  of  how  well  the  chord- 
wise  loading  irregularities  are  reduced  through  use  of  the  subdivision 
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process  is  shown  in  figure  3^  (reproduced  from  reference  10) •  The  re¬ 
duction  in  the  chordwise  load  variation  is  due  to  reducing  uhe  irregu¬ 
larities  of  the  velocity  potentials  shown  in  figure  35  (reproduced 
from  reference  10).  I-c  can  be  shown  that  subdividing  the  boxes  even 
further  will  produce  a  smooth  converged  velocity  potential  definition 
providing  smooth  chordwise  loadings  that  more  readily  simulate  the 
physical  f'ows.  However,  further  subdivision  would  result  in  an  ex¬ 
orbitant  amount  of  computer  time  required  to  achieve  the  smooth  load¬ 
ing  results.  Consequently  velocity  potential  smoothing  techniques 
were  investigated  in  the  hope  that  a  single  smoothing  process  could 
be  developed  and  applied  to  the  basic  box  pattern  velocity  potentials 
that  would  produce  results  similar  to  the  highly  subdivided  case. 

One  of  the  smoothing  techniques  investigated  and  available  within  the 
program  is  the  least-squares-error  surface  fitting  function  used  in 
the  modal  interpolation  portion  of  the  program.  The  program  data  op¬ 
tion  SMOOTH  =  .TRUE,  will  apply  this  surface  fitting  technique  to  the 
velocity  potentials.  An  accuracy  evaluation  was  made  on  this  smooth¬ 
ing  technique  and  it  was  found  that  it  produced  reasonable  results 
for  those  planforms  that  did  not  have  many  irregularities  in  the  box 
pattern  that  approximates  the  planform  leading  and  trailing  edges  as 
well  as  meeting  the  condition  of  having  a  small  angle  between  the  fore¬ 
most  Mach  line  and  the  planform  leading  edge.  However,  questionable 
results  were  achieved  whenever  the  above  conditions  were  not  met.  It 
was  also  determined  that  a  least-squares-error  smoothing  process  ap¬ 
plied  directly  to  the  velocity  potentials  produced  unreliable  loading 
results  at  the  ends  of  the  intervals  where  there  may  be  appreciable 
variations  in  the  original  velocity  potential  data.  Although  tne 
least-squares-error  does  fit  the  potentials  reasonably  well  over  the 
interval,  the  first  derivative  of  the  resulting  fitting  function 
(wnich  defines  the  major  portion  of  the  box  loadings)  is  not  well  de¬ 
fined  at  the  ends  of  the  intervals.  Consequently,  it  is  advisable  to 
exercise  caution  when  using  this  smoothing  option.  An  alternate 
smootning  option  is  available  to  the  user  by  specifying  CRDFIT  =  .TRUE, 
within  the  input  data.  This  smoothing  technique  was  devised  to  elimi¬ 
nate  the  need  of  taking  a  derivative  of  a  least-squares-error  poly¬ 
nomial  and  is  obtained  by  an  integration  process  that  provides  smooth¬ 
ing  on  the  average  over  the  interval. 

This  process,  as  shown  in  figure  36,  smooths  the  velocity  potentials 
of  individual  chordwise  strips  by  (l)  formulating  a  finite  difference 
derivative  distribution  from  the  basic  velocity  potential  data,  (2) 
performs  a  least-squares-error  fit  of  the  derivatives,  (3)  integrates 
the  least-squares-error  polynomial  to  obtain  the  smoothed  chordwise 
velocity  potential  distribution. 

It  should  be  noted  that  the  only  basis  for  applying  a  smoothing  opera¬ 
tion  to  the  velocity  potentials  lies  in  the  fact  that  the  velocity 
potentials  do  indeed  become  very  smooth  as  the  basic  box  size  is  in¬ 
creasingly  subdivided  into  many  small  sub-areas. 
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Comparisons  of  chordwise  loadings  obtained  using  smoothed  and  non- 
smoothed  velocity  potentials  are  presented  in  the  flutter  model 
correlation  studies  and  these  comparisons  indicate  that  reasonable 
load  distributions  may  be  obtained  using  the  program  option 
CRDFIT  =  .TRUE. 

PRELIMINARY  CHECKOUT  OF  THE  PROGRAM 

A  preliminary  checkout  of  the  program  has  been  accomplished  to 
evaluate  the  validity  of  the  program  output.  Since  there  are  no 
unsteady  pressure  data  available  for  the  supersonic  cases  then  only  a 
qualitative  evaluation  can  be  made. 

a.  Interaction  Loadings 

A  comparison  of  how  the  chordwise  loadings  on  mutual  interacting 
lifting  surfaces  are  affected  by  horizontal  separation  spacings 
is  presented  by  White  and  Landahl  in  reference  14.  Figure  37  re¬ 
produced  from  reference  14  displays  the  variation  on  chordwise 
loadings  for  various  separation  spacings  in  subsonic  flow.  For  a 
zero  gap  spacing  the  loadings  produced  at  the  wing-tail  junction 
exhibits  the  typical  logarithmic  singularity  identified  for  a  dis¬ 
continuous  downwash  distribution  encountered  on  a  wing-control 
surface  configuration  in  subsc'.i-  flow.  As  separation  begins, 
the  loadings  on  the  wing  near  the  wing  trailing  edge  change  from 
a  singularity  characteristic  to  a  loading  that  becomes  equal  to 
zero.  As  the  gap  spacing  becomes  large  the  loadings  at  the  wing 
trailing  edge  fall  off  as  the  square  root  of  the  distance  from 
the  trailing  edge,  meeting  the  Kutta  condition.  However,  the 
loadings  near  the  leading  edge  of  the  tail  change  from  a  logarith¬ 
mic  character  to  that  of  having  an  inverse  square  root  singularity 
that  is  more  highly  loaded  for  the  large  gap  than  it  is  for  the 
zero  gap  condition.  Although  the  analysis  of  reference  14  is 
applicable  to  subsonic  flows,  the  trends  should  exhibit  the  same 
behavior  in  three  dimensional  supersonic  flow  provided  that  the 
span  station  being  examined  lies  well  within  the  subsonic  commu¬ 
nication  region  where  the  lower  and  upper  surfaces  may  mutually 
affect  each  other. 

Numerical  results  were  obtained  to  evaluate  the  effect  of  gap 
spacing  in  supersonic  flow  for  the  45°  swept  wing-bail  planform 
used  in  the  Cornell  flutter  model  tests  shown  in  figure  38.  The 
gap  spacing  between  the  surfaces  was  obtained  for  integer  numbers 
of  boxes  within  wing  wake  diaphram.  The  fourth  coupled  mode  shape 
of  the  flutter  analysis  which  contains  a  large  amount  of  wing-tail 
angle  of  attack  in  a  streamwise  direction  was  selected  for  the 
data  comparison  of  White  and  Landahl  of  reference  l4. 

Figure  39  shows  the  chordwise  loadings  obtained  for  various  gap 
sizes  on  the  span  station  n  =  0,85  at  M  =  1.238.  The  analysis 
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chord  lies  well  wi'thin  the  subsonic  communication  region  and  the 
results  obtained  exhibit  the  same  trends  as  predicted  by  reference 
lb.  The  loading  in  the  trailing  edge  of  the  wing  change  from  a 
singular  characteristic  to  one  of  tending  to  zero.  Also  the  load¬ 
ings  on  the  tail  near  the  tail  leading  edge  increase  in  magnitude 
with  an  increase  in  gap  spacing  as  is  predicted  by  reference  lb. 

Figure  bO  is  a  plot  of  the  chordwise  loadings  for  the  same  mode 
shape  but  at  the  inboard  station  n  =  0.35.  Tlie  loadings  on  the 
wing  near  the  wing  trailing  edge  do  not  change  appreciably  with 
an  increase  in  gap  spacing.  This  is  to  be  expected  since  there 
is  only  a  very  small  region  where  subsonic  communication  can  exist 
between  the  lower  and  upper  wing  surfaces.  However  the  tail  lead¬ 
ing  edge  loads  do  increase  with  gap  spacing  as  more  of  the  wing- 
wake  region  becomes  available  for  upper  and  lower  surface  communica¬ 
tion. 

Velocity  Potential  Smoothing 

The  velocity  potential  smoothing  option  CRDFIT  =  .TRUE,  was  ap¬ 
plied  to  various  configurations  to  evaluate  its  accuracy  and  be¬ 
havior  characteristics.  For  velocity  potential  distributions 
that  were  initially  quite  smooth  the  smoothing  operation  provided 
a  slight  increase  in  smoothness  and  consequently  provided  a  chord- 
wise  loadings  that  were  more  smooth  on  the  average  than  was  ob¬ 
tained  without  the  smoothing  process.  No  large  variations  in  box 
loads  wore  observed  at  the  ends  of  the  intervals  as  had  been  the 
case  for  some  of  the  configurations  when  the  original  smoothing 
option  SMOOTH  *  .TRUE,  was  used.  The  largest  deviations  were  ob¬ 
served  to  occur  for  those  chordwise  strips  where  large  grad¬ 
ients  exist  in  the  original  velocity  potential  data. 

Figure  bl  shows  the  results  of  applying  smoothing  to  the  velocity 
potentials  on  a  highly  loaded  section  of  the  b5°  wing-tail  analysis 
configuration  of  the  Cornell  flutter  model  test.  The  velocity  po¬ 
tentials  and  resulting  loadings  are  shown  for  a  chord  that  is 
located  at  the  n  «  0,85  span  station.  The  fourth  coupled  mode 
shape  of  the  Cornell  modal  data  is  used  in  the  analysis.  The 
fourth  mode  contains  a  large  amount  of  wing  and  tail  angle  of 
attack  that  generates  the  large  variation  in  chordwise  distribu¬ 
tion  of  the  velocity  potentials. 

Also,  the  second  mode  shape  (primarily  a  differential  bending  mode 
between  wing  and  tail)  was  used  to  evaluate  the  smoothing  process 
at  the  same  span  station  and  the  results  for  this  condition  are 
shown  in  figure  b2.  It  is  apparent  that  no  large  deviations  re¬ 
sult  (even  for  a  highly  load  chord)  in  applying  the  smoothing  op¬ 
tion  and  also  provides  more  realistic  load  distribution  while 
retaining  the  original  overall  loading  characteristics. 
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c.  Normalwash  Examination 


Normalwash  distributions  in  the  wake  of  the  wing  were  examined 
for  smoothness  and  reasonableness  as  a  preliminary  exercise  prior 
to  performing  a  flutter  analysis  of  the  wing-tail  configurations 
for  various  horizontal  separation  spacings.  The  wing  of  the  45° 
sweep  Cornell  Configuration  15  was  used  in  the  evaluation.  The 
third  and  fourth  coupled  mode  shapes  of  the  Cornell  modes  data 
were  used  and  the  deflection  shapes  are  shown  in  figure  43  and 
figure  44.  The  normalwash  distributions  within  the  wing  wake  are 
shown  in  figure  45  and  figure  46. 

There  are  variations  within  the  wing-wake  normalwash  distributions 
as  would  be  expected,  however  there  are  no  large  excursions  and 
the  trends  are  reasonably  consistent  with  the  mode  shapes  causing 
these  distributions. 

Vertical  wash  distributions  were  obtained  for  various  vertical 
separation  values  and  are  displayed  in  figure  47 .  Again  there 
are  no  large  excursions  in  the  distributions  and  the  expected 
trends  of  having  a  diminishing  normalwash  with  an  increase  in 
vertical  separation  is  well  defined. 

5.  FLUTTER  MODEL  CORRELATION 

Flutter  model  test  data  used  in  this  correlation  study  were  obtained 
from  reference  2.  The  tests  were  conducted  in  the  Cornell  Aeronautical 
Laboratory  8'  x  8'  variable  density  Transonic  Wind  Tunnel.  The  test 
Mach  number  range  was  0.40  (  M  <f  1.24.  Of  the  many  test  points  ob¬ 
tained  there  are  only  five  points  that  lie  well  enough  within  the  super¬ 
sonic  range  to  allow  a  theoretical  analysis  to  be  performed.  The 
models  were  constructed  such  that  the  planform  shapes  were  identical 
for  both  wing  and  tail  configurations.  The  supersonic  tests  were  con¬ 
ducted  having  both  models  contained  in  the  same  plane  with  a  very 
small  horizontal  clearance  gap  between  the  models.  The  models  wei*e 
elastically  coupled  by  means  of  an  inter-connecting  torsion  spring 
simulating  a  fuselage  torsion,  degree  of  freedom. 

Test  flutter  modes  consisted  primarily  of  the  bending  modes  of  the 
ving-tail  surfaces  coupled  with  fuselage  torsion  degree  of  freedom. 

The  construction  details  of  che  models  are  well  discussed  within 
reference  2  and  will  not  be  repeated  here.  The  overall  dimensions  of 
the  45°  sweep  configuration  are  presented  in  figure  48  and  figure  49 
displays  the  geometry  of  the  60°  sweep  configuration  used  in  the 
flutter  model  test. 

a.  Vibration  Analysis  and  Flutter  Solution 

The  vibration  analysis,  based  on  the  elastic  axis  and  lumped  mass 
idealization,  was  carried  out  for  each  configuration  using  the 
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following  data  supplied  by  AFFDL: 

J  )  the  weight,  static  unbalance,  and  moment  of  inertia  of  the 
stations  along  the  model  elastic  axis,  40  percent  chord. 

(Tables  3,4  and  5) 

2)  the  first  two  uncoupled  cantilever  bending  and  the  first  un¬ 
coupled  cantilever  torsion  mode  of  the  wing  and  tail  (Tables 
6  and  7) 

3)  the  uncoupled  rigid  body  roll  mode  and  the  uncoupled  fuse¬ 
lage  torsion  mode  (Table  8) 

h)  the  generalized  coordinate  transformation  matrix  between  the 

eight  uncoupled  modes  and  the  five  coupled  modes, (Tables  9,10  &  11) 

The  five  coupled  modes  ured  for  flutter  analysis  were  calculated 
by  taking  matrix  product  of  the  uncoupled  modes  and  the  generalized 
coordinate  transformation.  These  coupled  mode  shapes  are  expressed 
as  bending  and  torsion  displacements  at  the  model  elastic  axis 
stations  and  as  the  angular  displacement  about  the  roll  axis  at 
the  fuselage  wing  and  tail  stations.  The  generalized  mass  matrix 
associated  with  these  coupled  modes  are  calculated  as  described 
in  reference  15.  As  the  program  input  for  calculating 
generalized  forces,  the  coupled  mode  shapes  were  evaluated  at 
the  arbitrary  points  along  the  chords  normal  to  the  model 
elastic  axis. 


The  generalized  forces  matrices  were  calculated  each  time  for  a 
Mach ‘number  and  a  few  reduced  frequencies.  A  chain  of  cubic  inter¬ 
polation  was  then  used  to  find  the  elements  of  the  generalized 
forces  matrices  at  many  intermediate  reduced  frequencies.  This 
substantially  reduces  the  cost  of  performing  a  flutter  analysis 
without  marked  loss  in  accuracy,  since  the  generalized  forces  vary 


smoothly  with  reduced  frequencies.  The  interpolation  routine  also 
scales  the  generalized  forces  matrix  by  -  * — 


The  flutter  speeds  and  frequencies  were  obtained  from  the  follow¬ 
ing  complex  eigenvalue  equation,  using  the  conventional  V-g 
solution  technique  with  complex  QR  algorithm: 


(  [M]  -  f[c]){q}  =  xTKj{q} 

n  _  *1 
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where 


b. 


[M]  is  the  generalized  mass  matrix 

fKj  is  the  generalized  stiffness  matrix;  the  elements 
are  given  by 


K. .  =  to.  M. . 

11  1  11 

uk  is  the  circular  frequency  of  the  i 
coupled  mode 


th 


[c]  is  the  interpolated  generalized  forces  matrix 
[c]  =  -  [Q],  =  -A- - 


-=r  i3  the  airforce  coefficient  equal  to 


p  b3 

20 


is 


/real 
the  s . 

{.imaginary 


part  of  the  complex 


eigenvalue 


is  the  artificially  added  structural  damping 
is  the  flutter  circular  frequencies  when  g  =  0 
is  the  flutter  speed  when  g  =  0 


Node  Line  Correlation 


The  procedures  described  in  the  previous  section  to  obtain  the 
coupled  mode  surface  deflections  were  applied  and  the  resulting 
node  line  comparisons  are  given  in  figure  50,  figure  51,  and 
figure  52  for  the  three  test-analysis  correlation  configurations. 

Only  the  node  lines  of  coupled  modes  number  three  (3)  and  four  (U) 
are  displayed  for  comparison  purposes.  Although  the  data  of  the 
test  report  (reference  2)  showed  good  correlation  between  the 
theoretical  and  experimental  natural  mode  frequencies,  the* node 
line  correlation  displayed  large  discrepancies  in  the  location  of 
the  modal  node  lines.  The  general  shape  of  the  theoretical  node 
lines  obtained  within  reference  2  compare  favorably  with  the  ex¬ 
perimental  node  lines  however  they  are  displaced  aft  toward  the 
trailing  edge  of  both  the  wing  and  tail  planforms.  Also,  the 
theoretical  node  lines  obtained  at  Boeing  using  the  Cornell  data 
are  displaced  even  further  aft  of  the  Cornell  theoretical  node 
line  positions  even  though  the  same  modal  data  were  used  in  the 
calculations.  It  should  be  noted  that  the  largest  discrepancies 
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occur  on  the  tail  for  each  of  the  configurations  evaluated.  It 
was  determined,  through  a  cross  plotting  process,  that  the  Boeing 
node  line  definitions  could  he  made  to  coincide  with  experimental 
node  lines  simply  hy  adding  a  constant  value  to  each  of  the  mode 
shapes.  The  values  to  be  added  to  the  wing  modes  have  different 
values  than  those  to  be  added  to  the  tail  mode  shapes.  As  a 
consequence  of  the  mismatch  in  the  mode  shape  definition  the 
flutter  correlation  study  vas  accomplished  in  two  parts.  The 
initial  evaluation  of  flutter  speeds  was  accomplished  using  the 
original  theoretical  mode  shapes  as  obtained  from  the  Cornell 
data,  and  the  second  correlation  study  was  accomplished  using  modi¬ 
fied  mode  chapes  that  provided  a  match  of  the  test  node  line  lo¬ 
cations  . 

c.  Flutter  Results  Using  Original  Modal  Data 

The  results  of  the  flutter  analysis  using  the  original  modal 
values  are  presented  in  Table  12  for  configuration  15  and  Table  13  and  ] 
presents  the  results  for  configurations  22  and  23.  Although  no 
test  data  are  available  to  correlate  the  effects  of  horizontal 
and  vertical  separations  in  supersonic  flow  theoretical  investi¬ 
gations  were  performed  to  determine  the  flutter  trends  for  these 
spacings.  The  separation  studies  were  conducted  only  on  the  1+5° 
sweep  planform  of  configuration  15. 

d.  Effects  of  Horizontal  Spacing 

The  effect  of  horizontal  separation  on  flutter  speed  is  shown  in 
figure  53  indicating  that  only  a  very  small  reduction  in  flutter 
velocity  is  obtained  for  relatively  large  values  of  horizontal 
separation.  No  test  data  are  available  for  the  supersonic  separa¬ 
tion  studies,  however  the  subsonic  test  results  of  configurations 
5  and  10  supports  this  trend  by  indicating  that  the  flutter  speeds 
are  only  slightly  affected  by  horizontal  separation.  It  should  be 
noted  that  the  lowest  flutter  speed  shown  within  the  horizontal 
separation  studies  was  obtained  for  the  case  of  =  0. 

This  condition  is  one  that  simulates  the  existence  of  an  aerodynamic 
seal  or  a  no  flow  through  boundary  between  the  wing  and  tail  sur¬ 
faces,  Figure  39  shows  that  although  the  loadings  on  the  tail 
leading  edge  are  reduced  by  this  no  flow  thru  condition,  there  is 
a  large  overturning  moment  developed  on  the  wing  causing  the  wing 
to  respond  more  readily  than  it  does  for  the  open  gap  condition. 

e.  Effects  of  Vertical  Separation 

The  effects  of  vertical  separation  on  flutter  speed  and  frequency 
ratio  are  shown  in  figure  5^  for  a  constant  Mach  number,  tunnel 
density  and  fuselage  torsional  frequency.  It  should  be  noted  that 
the  flutter  mode  shape  is  changing  with  an  increase  in  vertical 
separation  as  indicated  by  the  increase  in  the  flutter  frequency 
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ratio  and  the  analysis  configuration  is  taking  on  the  character¬ 
istics  of  single  surface  bending-torsion  flutter  mode. 

Thickness  Effects 


ir 


J.  • 


The  inclusion  of  airfoil  thickness  effects  within  the  flutter 
studies  resulted  in  only  a  very  small  reduction  in  flutter  speeds 
for  all  three  study  configurations  as  indicated  by  analysis  no.  9, 
ll*,  and  19.  The  chordwise  loadings  on  both  the  wing  and  tail  sur¬ 
faces  increased  in  the  region  of  the  leading  edges  and  decreased 
in  value  over  the  aft  portion  of  the  surfaces.  No  large  changes 
in  flutter  frequency  or  mode  shape  resulted  from  including  thick¬ 
ness  effects  within  the  analysis.  This  was  anticipated  since  the 
streamwise  thickness  ratios  are  quite  small  and  cannot  appreciably 
change  the  loadings  or  center  of  pressure  locations . 

g.  Mach  Number  Effects 


Non-dimensional  values  of  flutter  frequency  and  flutter  velocity 
are  presented  as  a  function  of  Mach  number  in  figure  55,  figure  56 
and  figure  57  for  the  three  study  configurations.  The  analytical 
results  obtained  for  M=1.35  of  configuration  15  and  for  M=1.414  for  con¬ 
figurations  22  and  23  were  obtained  by  using  extrapolated  density  data 
and  may  be  somewhat  in  error.  The  overall  trends  as  shown  in  Figures  84 
and  85  indicate  that  the  flutter  speed  ratios  do  not  change  by  any  large 
**  amount  over  the  Mach  number  range  investigated.  The  resulting 
flutter  modes  consisted  of  a  strong  coupling  between  the  wing-tail 
bending  modes  with  the  body  torsion  mode  for  most  of  the  analysis 
cases  with  the  exception  of  configuration  23  for  M  =  1.232.  Con¬ 
figuration  23  at  M  =  1.232  had  a  flutter  mode  that  contained  an 
unusual  amount  of  the  fourth  coupled  mode.  Various  parameters  such 
as  stiffnesses  and  densities  were  changed  to  try  to  identify  what 
was  the  cause  for  this  behavior  however,  the  flutter  mode  remained 
essentially  the  same  throughout  the  parameter  variation  study. 

The  error  between  theoretical  and  experimental  flutter  frequencies 
is  quite  small  for  all  of  the  analysis  cases  with  the  exception  of 
configuration  23  at  M  =  1.232.  However,  the  errors  between  theo¬ 
retical  and  experimental  flutter  velocities  are  larger  than  desired 
and  this  error  could  not  be  greatly  reduced  with  normal  changes  in 
various  parameters.  The  errors  in  flutter  velocities  are  of  the 
same  magnitude  for  the  subsonic  correlation  studies  of  reference  15, 
at  least  for  configurations  1 6  and  22,  as  reported  by  Albano, 

Perkinson  and  Rodden. 

Results  of  the  subsonic  studies  are  shown  in  figure  58  and  figure  59 
that  have  been  reproduced  from  reference  15.  Since  rather  large 
errors  are  incurred  for  both  the  subsonic  and  supersonic  studies 
using  the  original  modal  data,  an  attempt  was  made  to  modify  Uio  ■'>,)< 
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shapes  such  that  the  theoretical  node  lines  would  coincide  with  the  experimen¬ 
tal  node  lines  and  thus  evaluate  the  sensitivity  of  mode  shape  variations  on 
flutter  velocity. 

h .  Effect  of  Mode  Shape  Variations 

Modifications  to  the  mode  shapes  were  accomplished  by  allowing  a  rigid  body 
rotation  about  an  axis  parallel  to  the  node  lines.  The  rotation  axis  was 
located  at  the  attachment  fittings  at  wing  and  tail  root  sections.  The  direc¬ 
tions  of  the  axes  were  essentially  perpendicular  to  the  air  stream  and  could 
be  thought  of  as  allowing  a  pitching  flexibility  of  the  support  system  to 
take  place.  The  amount  of  rotation  allowed  in  each  of  the  modes  was  just 
enough  to  shift  the  theoretical  node  line  position  to  coincide  with  the  experi¬ 
mental  node  line.  A  typical  value  of  rigid  body  rotation  that  was  applied 
within  the  mode  modification  amounted  to  a  rotation  value  of  0.0024  radians 
which  would  change  the  deflection  at  the  theoretical  node  line  position  by 
0.05  inches.  It  should  be  noted  that  all  mode  shapes  were  normalized  to  a 
maximum  value  of  1.0.  Consequently  the  mode  shape  changes  that  were  applied 
amounted  to  values  that  were  approximately  equal  to  5%  of  the  maximum  deflec¬ 
tion.  A  new  set  of  coupled  modes,  generalized  masses,  generalized  stiffness, 
and  airforces  were  produced  for  use  within  the  flutter  solution  program. 

Table  15  presents  the  results  obtained  for  configuration  15,  configuration  22 
and  configuration  23.  Use  of  the  modified  mode  shapes  produce  results  that 
are  closer  to  the  experimental  results  and  may  indicate  that  the  original 
modal  data  and  encompassing  assumptions  may  not  be  completely  correct. 

There  are  no  means  available,  within  the  modal  data,  to  verify  how  much  the 
support  system  flexibility  could  have  contributed  to  the  coupled  modes,  how¬ 
ever,  the  analysis  indicates  that  the  prediction  of  flutter  velocities  are 
highly  influenced  by  small  changes  in  the  mode  shape  distribution. 


i .  Discussion  of  Results 


The  theoretical  trends  obtained  in  the  flutter  analysis  are  similar  to  those 
obtained  during  the  experimental  flutter  tests.  Horizontal  separation  has 
very  little  effect  on  flutter  speeds  for  the  coplanar  cases  and  there  is  an 
Increase  in  flutter  speed  with  vertical  separation.  The  flutter  modes  were 
found  to  be  composed  of  wing  bending  coupled  with  fuselage  torsion  for  all  of 
the  analysis  cases  with  the  exception  of  configuration  23  at  M  «  1.232  where 
the  flutter  mode  contained  a  high  frequency  torsion  mode.  Although  the  theo¬ 
retical-experimental  trends  are  similar,  it  appears  that  the  theory  is  conser¬ 
vative  and  may  be  due  to  unknown  effects  such  as  shocks  attached  on  the  wing 
and  tail  leading-trailing  edges,  or  may  be  due  to  flow  blockage  over  the  tail 
due  to  the  wings  presence. 
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There  are  other  possible  sources  of  error  that  may  contribute  to 
the  conservatism,  one  of  which  is  due  to  the  numerical  solution 
technique  itself.  Probably  the  most  critical  approximation  made 
within  the  "Mach  Box"  (or  any  other  Box  pattern  program)  is  the 
assumption  that  normal  washes  are  of  constant  value  over  the 
invididual  box  areas.  This  approximation  may  apply  very  well  over 
the  lifting  surface  planform  where  the  normal  wash  distributions 
are  relatively  small  and  are  only  functions  of  the  modal  deflec¬ 
tions,  However,  relatively  large  errors  may  result  in  predicting 
the  normal  wash  distributions  over  the  "off-wing"  diaphram  region 
near  the  side  edges  of  the  planform  where  the  normal  wash  varies 
in  proportion  to  the  square  root  of  the  distance  from  the  planform 
edge. 

The  resulting  predicted  normal  washes  will  invariably  be  larger 
than  exist  in  reality  and  will  result  in  predicting  velocity  po¬ 
tentials  (or  loadings)  over  the  lifting  surface  that  are  larger 
than  those  existing  in  the  true  flow  field. 

There  are  possible  sources  of  error  in  the  basic  mathematical  model 
of  the  coupled  modes  used  in  the  flutter  analysis  that  may  con¬ 
tribute  to  the  conservative  velocities.  Five  coupled  mode  shapes 
were  obtained  for  flutter  analysis  using  only  eight  primative 
"uncoupled"  modes  as  a  basic  set  of  raodeB  to  establish  the  final 
coupled  mode  definition. 

The  eight  primitive  modes  consisted  of  two  cantilevered  "uncoupled" 
bending  modes  and  one  "uncoupled"  torsion  for  each  surface  plus 
one  rigid  body  roll  mode  and  one  fuselage  differential  torsion  mode. 
No  allowance  was  made  for  pitch  flexibility  of  the  support  system. 

The  cantilever  analysis  considered  each  lifting  surface  as  a  straight 
beam  of  variable  stiffness  and  inertia  distribution  that  were  al¬ 
lowed  to  deform  in  bending  and  torsion.  This  "elastic  axis"  as¬ 
sumption  requires  all  sections  perpendicular  to  the  elastic  axis  to 
deform  as  a  rigid  body  which  may  be  questionable  when  applied  to 
the  lifting  surface!  used  in  the  present  analysis  that  do  not  have  a 
well  defined  structural  axis  nor  have  rigid  ribs  to  enforce  the 
rigid  deformation  requirements  of  the  cantilever  analysis. 

There  is  no  method  available  to  determine  if  any  errors  (or  even 
less,  the  modal  distribution  of  errors)  exist  within  the  original 
modal  data,  however,  the  flutter  analysis  does  indicate  that  the 
accuracy  of  predicting  flutter  velocities  is  quite  sensitive  to 
small  changes  in  mode  shapes  of  highly  swept  low  aspect  ratio  con¬ 
figurations. 
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TABLE  2  MEASURED  CANTILEVER  COUPLED-MODE  FREQUENCIES 


CONFIGURATION 

NUMBER 

MODEL 

1st  COUPLED 

MODE 

(CPS) 

2nd  COUPLED 

MODE 

(CPS) 

3rd  COUPLED 

MODE 

(CPS) 

U)  2 

U)3 

CO  2 

15 

45P3 

(WING) 

17.6 

67.3 

87.5 

.262 

1.300 

45P6 

(TAIL) 

21.7 

83.5 

106.8 

.260 

1.279 

22  &  23 

60P4 

(WING) 

16.1 

60.7 

94.7 

.265 

1.560 

60P8 

(TAIL) 

19.0 

69.4 

116.0 

.274 

1.671 
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TABLE  3  MASS  PROPERTIES  OF  CONFIGURATION  15 
(  45°  SWEPT  MODEL  ) 


ELASTIC 

AXIS 

STATION 

NO. 

H  =LENGTH 

INCHES  FROM 

ROLL  AXIS 

WEIGHT  (LBS) 

STATIC  UNBALANCE 

(LB-IN) 

MOMENT  OF  INERTIA 
(LB-IN2  ) 

WING 

(45P3) 

TAIL 

(45P6) 

WING 

(45P3) 

TAIL 

(45P6) 

WING 

(45P3) 

TAIL 

(45P6) 

1* 

0. 

.01 

.01 

0. 

0. 

0. 

0. 

2* 

2.83 

.01 

.01 

0. 

0. 

0. 

0. 

3 

5.09 

4.31 

4.33 

-4.11 

-4.18 

12.84 

13.06 

4 

7.26 

3.370 

3.39 

-  .164 

-  .158 

9.86 

9.06 

5 

12.25 

1.003 

.994 

.68 

.799 

13.98 

14.79 

6 

17.44 

.436 

.498 

.560 

.680 

9.79 

11.68 

7 

21.72 

.304 

.325 

.335 

.367 

6.14 

6.66 

8 

25.77 

.304 

.325 

.335 

.367 

6.14 

6.66 

9 

29.82 

.304 

.325 

.335 

.367 

6.14 

6.66 

10 

33.87 

.304 

.325 

.335 

.367 

6.14 

6.66 

11 

38.19 

.313 

.334 

.601 

.655 

5.85 

6.37 

12 

42.69 

.186 

.201 

.782 

.847 

4.84 

5.27 

13 

47.56 

.110 

.118 

.718 

.785 

5.3P 

5.89 

*  STATIONS  ADDED  IN  THIS  STUDY 
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TABLE  4  MASS  PROPERTIES  OF  CONFIGURATION  22  &  23 


(  60°  SWEPT  MODEL  ) 


ELASTIC 

AXIS 

STATION 

NO. 

H=  LENGTH 

INCHES  FROM 

WEIGHT  (LBS) 

STATIC  UNBALANCE 
(LB-IN) 

MOMENT  OF  JNERTIA 
(LB-IN*) 

ROLL  AXIS 

i 

WING 

(60P4) 

TAIL 

(60P8) 

WING 

(60P4) 

TAIL 

(60P8) 

WING 

(60P4) 

TAIL 

(60P8) 

1* 

0. 

.01 

.01 

0. 

0. 

0. 

0. 

2 

5.0 

1.43 

1.5 

-2.35 

-2.70 

18.3 

20.0 

3 

8.34 

2.72 

2.77 

-2.43 

-2.67 

8.74 

10.06 

4 

13.45 

2.32 

2.26 

0.254 

.211 

8.21 

8.82 

5 

18.55 

.871 

.916 

.616 

.686 

11.73 

13.55 

6 

23.75 

.384 

.421 

.538 

.597 

8.19 

10.00 

7 

28.98 

.393 

.426 

.396 

.467 

7.35 

8.29 

8 

34.93 

,393 

.426 

.396 

.467 

7.35 

8.29 

9 

40.45 

.312 

.336 

.481 

.558 

5.55 

6.28 

10 

45.5 

.241 

.261 

.693 

.784 

4.50 

5.14 

11 

50.5 

.170 

.185 

.725 

.820 

4.27 

4.93 

12 

54.9 

.0746 

.0843 

.447 

.513 

2.94 

3.43 

13 

59.74 

.0559 

.0665 

.426 

.509 

3.41 

4.10 

*  ADDED  STATION  IN  THIS  STUDY 
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ABLE  5  ROLL  INERTIA  OF  FUSELAGE  OF  CORNELL  WING-HORIZONTAL  TAIL 
FLUTTER  MODELS 


ROLL  INERTIA  OF 

ROLL  INERTIA  OF 

CONFIGURATION 

FUSELAGE  AT  WING 
LOCATION  ( lb-1 nz) 

FUSELAGE  AT  TAIL 

NUMBER 

LOCATION  (lb-in2) 

TABLE  8  UNCOUPLED  RIGID  BODY  ROLL  AND  FUSELAGE 
TORSION  MODE  SHAPES  OF  THE  CORNELL  WING- 
HORIZONTAL  TAIL  FLUTTER  MODELS 


CONFIGURATION 

NUMBER 

RIGID  BODY 

WING  TAIL 

ROLL 

FREQ. 

(CPS) 

FUSELAGE  TORSION 

WING  TAIL 

15 

1.0 

1.0 

4.0 

-.972 

1.0 

33-0 

22 

1.0 

1.0 

4.7 

-.954 

1.0 

24.7 

23 

1.0 

1.0 

4.7 

-.954 

1.0 

30.7 

TABLE  7  UNCOUPLED  CANTILEVER  MODE  SHAPES  OF  CONFIGURATIONS 


22  and  23  (60°  SWEPT  MODELS) 


ELASTIC 

AXIS 

STATION 

NUMBER 

11  -LENGTH 
(INCHES 

FROM 

ROLL  AXIS) 

1 

0. 

2 

5.0 

3 

8.34 

4 

13.45 

5 

18.55 

6 

23.75 

7 

28.98 

8 

34.93 

9 

40.45 

10 

45.49 

11 

50.5 

12 

54.9 

13 

59.74 

WING  AND  TAIL 


1st  CANTILEVER 

BENDING 

2nd  CANTILEVER 

BENDING 

1st  CANTILEVER 

TORSION 

FREQ. (Hz) 

WING  TAIL 

FREQ. (Hz) 

WING  TAIL 

FREQ. (Hz) 

WING  TAIL 

17.5  20.4 

70.0  80.4 

58.7  71.8 

0. 

0. 

0. 

0. 

0. 

0. 

0.001 

-0.006 

0.002 

0.004 

-0.024 

0.004 

0.011 

-0.055 

0.007 

0.024 

-0.108 

0.022 

0,069 

-0.209 

0.103 

0.178 

-0.325 

0.215 

0.323 

-O.330 

0.304 

0.480 

-0.207 

0.383 

0.650 

0.054 

0.484 

0.811 

0.420 

0.630 

1.00 

1.00 

1.00 

COUPLED  RIGID  BODY  ROLL  AND  FUSELAGE 
RSION  MODE  SHAPES  OF  THE  CORNELL  WING- 
RIZONTAL  TAIL  FLUTTER  MODELS 


CONFIGURATION 

NUMBER 

RIGID  BODY 

WING  TAIL 

ROLL 

FREQ. 

(CPS) 

FUSELAGE  TORSION 

KING  TAIL  (Qp5) 

15 

1.0 

1.0 

4.0 

-.972 

1.0 

33-0 

22 

1.0 

1.0 

4.7 

-.954 

1.0 

24.7 

23 

1.0 

1.0 

4.7 

-.954 

1.0 

30.7 

TABLE  9 


GENERALIZED  COORDINATE  TRANSFORMATION 
MATRIX  FOR  CONFIGURATION  15 


GENERALIZED  COORDINATES 


COUPLED 

MODES 

RIGID  BODY 
ROLL 

FUSELAGE 

TORSION 

Wing  Cantilever 

Tail 

Cantilever 

FREQ. 

CPS 

FIRST* 

BENDING 

SECOND* 

BENDING 

FIRST 

TORSION 

FIRST* 

BENDING 

SECOND 

BENDING 

FIRST 

TORSION 

1 

6.39 

0 

1.598 

-0.04 

0.0800 

1.0 

-0.023 

0.043 

3.9 

2 

0.0571 

0.0571 

-2.671 

0.002 

-0.1962 

1.0 

-0.009 

0.0562 

17.6 

3 

-0.0685 

0.0186 

0.810 

0.025 

0.0852 

1.0 

0.007 

0.0721 

25.3 

4 

-0.0165 

-0.1407 

-0.637 

-0.136 

-0.2034 

1.0 

0.085 

0.1511 

42.8 

5 

0.765 

-1.045 

-29.86 

15.56 

44.92 

1.0 

0.886 

1.219  , 

72.2 

TABLE  10  GENERALIZED  COORDINATE  TRANSFORMATION 
MATRIX  FOR  CONFIGURATION  22 


GENERALIZED  COORDINATES 


CulinrrU|RIGIU  BODY  FUSELAGE^ - --ng  Canti1ever  I  Tail  Cantilever 


FIRST 

FIRST* 

SECOND* 

TORSION 

BENDING 

BENDING 

3.731 

0 

1.446 

-0.0605 

0.106 

1.0 

-0.0424 

0.0676 

4.7 

0.0348 

0.0718 

-2.277 

0.00308 

-0.24 

1.0 

-0.0143 

0.0874 

16.3 

-0.0564 

0.0225 

0.901 

0.0339 

0.123 

1.0 

0.0115 

0.106 

21.6 

-0.0128 

-0.122 

-0.634 

-0.108 

-0.156 

1.0 

0.0921 

0.164 

31.4 

-4.238 

4.705 

128.401 

-171.801 

-159.55 

1.0 

-8.234 

-6.33 

67.2 

*  Denotes  the  Coordinate  Has  Been  Scaled  to  ^ —  ,  br  =  8.0  in. 


TABLE  11  GENERALIZED  COORDINATE  TRANSFORMATION 
MATRIX  FOR  CONFIGURATION  23 


GENERALIZED  COORDINATES 


fail  Cantilever 


D  RIGID  Body  FUSELAGE  Wing  Cantilever _ 

ROLL  TORSION  FIRST*  SECOND*  FIRST  FIRST*|SECOND*|  FIRST 
BENDING  BENDING  TORSION  BENDING  BENDING  TORSION 


3.7306  0.0  1.4456  -0.0606  0.1062  1.0  -0.0424  0.0676  4.7 


2  0.0422  I  0.0432  I  -2.53581  -0.0051  I  -0.2741  I  1.0  -0.0121  0.0890  16.9 


-0.0507  0.0126  0.7009  0.0280  0.0971  1.0  0.013  0.106721.9 


-0.0110  -0.1544-0.7032  -0.2223-0.261  1.0  0.1682  0.219137.0 


-2.1008  2.4649  63.0240  -82.4125  -76.6715  1.0  -6.9564-5.341767.6 


Denotes  the  Coordinate  Has  Been  Scaled  to 
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TABLE  12  ANALYTICAL  FLUTTER  RESULTS  ON  CONFIGURATION  15 
(ORIGINAL  MODES, 45°  SWEPT  MODELS , U>0/U>h=  1  .765) 


■ 

P 

THEORETICAL  RESULTS 

TEST  RESULTS 

ANALYSIS 

M 

H 

U 

»lug/ft3 

THICKNESS 

Vf 

Vf 

u> 

Vf 

n 

NUMBER 

H 

19 

(xIO4) 

CORRECTION 

ft/see 

•f 

cps 

•s* 

ft/l*c 

'f 

cps 

b  ,«UdVM~ 

■ 

1.121 

.021 

0.0 

44.5 

4.0 

1150. 

22.7 

.830 

.688 

■ 

1.121 

.166 

0.0 

NO 

*858.2 

873.3 

24.6 

24.08 

.6190 

.6043 

.7455 

.7297 

1.238 

.021 

0.0 

50.9 

3.5 

1254. 

26. 

.846 

.788 

2 

0.0 

0.0 

NO 

976. 

974.9 

25.70 

25.49 

.6592 

.6584 

.7789 

.7724 

3 

.167 

NO 

1006.8 

1001.1 

26.5 

26.27 

.6799 

.6760 

.8035 

.7959 

4 

.505 

H|l|| 

NO 

1001.2 

1001.9 

26.08 

25.91 

.6762 

.6766 

.7903 

.7851 

5 

1.0 

B 

NO 

991. 

990. 

25.58 

25.37 

.6693 

.6686 

.7752 

.7688 

B 

1.238 

.505 

0.0 

50.9 

3.5 

NO 

1001.2 

1001.9 

26.08 

25.91 

.6762 

.6766 

.7903 

.7851 

n 

.505 

0.0835 

NO 

1061.4 
1057. E 

26.6 

26.34 

.7168 

.7142 

.8042 

.7982 

B 

.505 

.167 

NO 

1128.8 

1123.2 

27.5 

27.23 

.7623 

.7585 

.8325 

.8252 

8 

.505 

.505 

NO 

1430.8 

1436.7 

35.09 

34.69 

.9662 

.9702 

1.0635 

1.0512 

1.238 

.167 

m 

50.9 

3.5 

NO 

1006.8 

1001.1 

26.5 

26.27 

.6799 

.6760 

.8035 

.7955 

.167 

YES 

964. 1 
959.5 

26.1 

25.91 

.6511 

.6480 

.7916 

.7851 

_ 

10 

1.35 

.167 

** 

59.3 

** 

3.0 

NO 

1021.3 

1023.7 

26.4 

26.32 

.6369 

.6403 

.8013 

.7975 

NO  TEST  RESULTS 

AVAILABLE 

*  THE  FIRST  LINE  OF  THE  THEORETICAL  RESULTS  FOR  EACH  CASE  DENOTES  THE 
RESULTS  IN  WHICH  THE  SMOOTHING  TECHNIQUE  (REFINEMENT)  IS  APPLIED 

**  ESTIMATED  VALUE 
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TABLE  13  ANALYTICAL  8.  TEST  FLUTTER  RESULTS  ON  CONFIGURATION  22 


{ORIGINAL  MOOES,  60°  SWEPT  MODELS, Wq/O)^!  .428) 


ANALYSIS 

NWBER 

■ 

■ 

■ 

P 

THEORETICAL  RESULTS 

TEST  RESULTS 

1 

B 

PM 

V 

*lug/n3 

THICKNESS 

Vf 

ft/»ec 

Vf 

JL 

we 

V, 

(U 

■ 

(xio4) 

CORRECTION 

nf 

ep» 

nf 

ep* 

11 

1.151 

.19 

«* 

53.9 

*3.3 

NO 

9)6.6 

21.67 

.6022 

.8773 

NO  TEST  RESULTS 

882.7 

21.14 

.5806 

.8558 

AVA1 

LABLE 

■ 

1.24 

.0257 

D 

65.3 

2.7 

1317 

20.8 

.708 

.84? 

D 

0.0 

Moj 

NO 

1089.0 

22.19 

.6513 

■  89B4 

H 

1073.7 

21.96 

.6422 

.8892 

D 

.19 

BIS 

NO 

1139.8 

23.70 

.6817 

.9594 

H 

1058.7 

22.46 

.6332 

.9091 

13 

1.24 

.19 

0.0 

65. 3 

2,7 

NO 

1139.8 

23.70 

.6817 

.9594 

1058.7 

22.46 

.6332 

.9091 

14 

.19 

0.0 

YES 

1101.9 

1035.2 

23.23 

22.22 

.6591 

.6191 

.9406 

.8994 

15 

1.414 

.197 

Ell 

*116. 

** 

1.52 

NO 

1282.6 

22.02 

.5755 

.8916 

NO  TEST  RESULTS 

■ 

1238.0 

21.61 

.5555 

.8764 

AVAILABLE 

*  THE  FIRST  LINE  OF  THE  THEORETICAL  RESULTS  FOR  EACH  CASE  DENOTES  THE  RESULTS 
IN  WHICH  THE  SMOOTHING  TECHNIQUE  (REFINEMENT)  IS  APPLIED 

**  ESTIMATED  VALUE 
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TABLE  14  ANALYTICAL  &  TEST  FLUTTER  RESULTS  ON  CONFIGURATION  23 


(ORIGINAL  MODES,  60°  SWEPT  MODELS, U)0/U)h  =  1,?75^ 


■ 

■ 

H 

■ 

E 

THEORETICAL  RESULTS 

TEST  RESULTS 

n 

FVVV.  1  ->  i  0 

% 

'MHI 

l 

MUMER 

l 

■ 

iiSH 

THICKNESS 

v, 

H 

Mi 

U 

Cf 

vf 

I  1 

■ 

313  I 

H 

■ 

CORRECTION 

ft/»»c 

ft/m 

Hi 

DPI 

1.161 

H 

ifg 

32.7 

n 

HR 

1214. 

23.5 

.826 

.765 

16 

in 

Bra 

I 

NO 

.6879 

.9388 

■ 

■ 

■ 

909.0 

ti) 

.6179 

.7956 

1.232 

m 

36.7 

4.8 

1292. 

24.2 

.829 

.788 

17 

^l§§ 

0.0 

NO 

1009.3 

2S.<3 

.64/6 

.8283 

1 1 

950.2 

24.26 

.6096 

.7903 

18 

NO 

1063.2 

m 

.6821 

Hfl 

1094.8 

eh 

.7024 

Bfl 

Mi 

m 

EM 

Ml 

18 

1.232 

m 

36.7 

a 

NO 

1063.2 

30.95 

1.008 

■ 

■ 

ME 

■ 

1.011 

■ 

■ 

19 

.19 

HqI 

■g 

YES 

1053.9 

.6763 

1  V 

1 

■ 

.. 

B 

1086.9 

H 

.6970 

jl  .003 

1 

1 

Ml 

20 

1.414 

.19 

[U 

m 

NO 

1094.8 

20.77 

.5703 

m 

NO  TEST  RESULTS 

■ 

II 

■ 

1014.4 

23.50 

.5138 

fi 

AVAILABLE 

*  THE  FIRST  LINE  OF  THE  THEORETICAL  RESULTS  FOR  EACH  CASE  DENOTES  THE  RESULTS 
IN  WHICH  THE  SMOOTHING  TECHNIQUE  (REFINEMENT)  IS  APPLIED 

**  ESTIMATED  VALUE 
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TABLE  15  ANALYTICAL  &  TEST  FLUTTER  RESULTS  ON  CONFIGURATIONS  15,22,23 


(MODIFIED  MODES,  45°  and  50°  SWEPT  MODELS) 


ANA. 

NO. 


13 


22 


CONFIG 

NO. 


18 


23 


15 


22 


23 


1.238 


1.24 


1.232 


«h 

X 

Z 

7s 

u 

o 

THEORETICAL  RESULTS 

TEST  R. 

.SULTS 

«1og/ft3 

(xlO4) 

MOO€S 

modified 

Vf 

ft/»ec 

cp» 

>>,VP 

U 

<*>• 

ft/lK 

'V 

ep» 

Vf 

U 

Jt 

1.765 

.021 

0.0 

50.9 

3.5 

1254. 

26.0 

.846 

.788 

.167 

0.0 

NO 

*1006.8 

26.5 

.6799 

.8035 

1001.1 

26.27 

.6760 

.7959 

.167 

0.0 

YES 

1103.6 

27.6 

.7454 

.8363 

1094.3 

27.26 

.7389 

.8260 

1.428 

0.0 

65.3 

2.7 

1317. 

20.8 

.788 

.042 

0.0 

HO 

1139.8 

23.70 

.6817 

.9594 

1058.7 

22.46 

.6332 

.9091 

i^| 

0.0 

YES 

1205.4 

24.259 

.7210 

.9821 

1 

1138.0' 

23.284 

.6806 

.9427 

IBHH 

0.0 

36.7 

4.8 

1292. 

24.2 

.829 

.788 

1 

0.0 

NO 

1063.2 

30.95 

.6821 

1.008 

I 

1094.8 

31.03 

.7024 

1.011 

I 

0.0 

YES 

1083. 5 

31.17 

.8951 

1.015 

ii 

1123.5 

30.98 

.7209 

1.009 

MM* 

mmm 

THE  FIRST  LINE  OF  THE  THEORETICAL  RESULTS  FOR  EACH  CASE  DENOTES  THE  RESULTS 
IN  WHICH  THE  SMOOTHING  TECHNIQUE  (REFINEMENT)  IS  APPLIED 


6  UNSTEADY  AERODYNAMIC  APPLICATIONS 

A  typical  wing-tail  configuration  has  been  specified  by  the  Advisory 
Group  for  Aeronautical  Research  and  Development  (AGARD)  of  NATO  Cor  use  in 
comparison  of  various  aerodynamic  theories  and  experimental  results.  The 
basic  configuration  is  shown  in  Figure  60. 

As  shown  in  the  figure,  all  dimensions  have  been  normalized  on  Che  basis  of 
semi-span.  In  order  to  test  and  demonstrate  the  various  capabilities  of  the 
computer  program,  three  variations  of  the  basic  configuration  have  been 
used.  These  may  be  summarized  as  follows: 


Configuration  Characteristics 

Number 

L  Wing  alone (with  flow  field  sampling) 

2  Wing  and  tail  with  very  small  longitudinal  separation 
(xWTE  “  XTLE^Root  =  0,25 

3  Wing  and  tail  with  moderate  longitudinal  and  vertical 
separation, -  r„B>Root  -  0.45  and(sT  -  z„)Root 

=  0.5,  and  tail  dihedral, ^  =  30°. 

For  the  case  of  a  wing  alone,  the  perturbation  velocity  components  (longi¬ 
tudinal  wash,  sidewash  and  upwash)  have  been  sampled  in  the  flow  field  for 
some  local-chord  lengths  behind  the  wing.  The  first  sampling  chord  is 
located  at  a  spanwise  distance  of  0.056s  from  the  centerline  and  a  vertical 
distance  above  the  wing  plane  of  0.5s.  The  second  is  located  at  a  spanwise 
distance  of  0.72sand  a  vertical  distance  of  0.1s. 

In  addition  to  the  perturbation  velocities,  the  chordwise  pressure  distri¬ 
butions,  generalized  forces,  and  generalized  aerodynamic  coefficients  have 
been  determined  for  various  reduced  frequency-Mach  number  combinations. 
Table  16  summarizes  the  cases  analyzed. 

Figures  61  through  66  give  the  values  of  longitudinal  wash,  sidewash  and 
upwash  plotted  along  the  sample  chord  for  some  local-chord  lengths  behind 
the  wing.  These  values  are  computed  by  applying  the  smoothing  technique 
for  two  reduced  frequencies,  kg,  equal  to  0.75  and  1.5  based  on  the  wing 

semi-span  for  antisymmetric  wing  bending  and  torsion  modes.  These  mode 
shapes  are  also  given  in  Table 16. 

Sample  plots  of  the  pressure  distributions  obtained  on  the  wing  for  M  = 
1.56  are  given  in  Figures  67  through  74.  The  pressures  are  plotted  at 
Chord  1,  y/s  =  0.056  and  Chord  7,  y/s  =  0.75.  Values  are  given  for 
reduced  frequencies  k  =  0.75  and  1.5  based  on  a  wing  semi-span, s,  and  are 
computed  by  applying  the  smoothing  technique. 
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TABLE  16  SUMMARY  OF  AGARD  CONFIGURATION  CASES 


Case 

Config. 

1 

M 

k 

s 

Surface 

Mode  1 

Mode  2 

Sym1 .  or  Anti  . 

1-1 

1 

1.56 

0.75 

1.5 

Wing 

2 

z  *  y 

z  *  xy 

Symmetric 

1-2 

1 

1.56 

0.75 

1.5 

Wing 

2 

z  *  y 

z  ■  xy 

Antisym. 

2-1 

2 

1.414 

0.75 

1.5 

Wing  & 
Tail 

Wing 

2 

z  =  y 

Tail 

2 

z  *  y 

Wing 

z  =  xy 

Tail 

*  i 

z  =  x  y 

2-2 

2 

3.162 

0.75 

1.5 

Wing  & 
Tail 

Symmetric 

2-3 

3 

1.56 

0.75 

1.5 

Wing  & 
Tail 

Antisym , 

*  The  x'~  reference  location  for  the  tail  torsion  mode  i9  taken  to  be  the 
root  of  the  tail  leading  edge. 


The  generalized  unsteady  aerodynamic  coefficients  at  kg  =  0.75  and  1.5 
for  these  antisymmetiric  wing  bending  and  torsion  modes  are  given  in 
Table  17.  The  lower  values  which  are  superscripted  by  *  show  the  quantities 
in  which  the  smoothing  technique  is  applied. 

The  pressure  distributions  obtained  for  the  wing-horizontal  tail  configuration 
with  very  small  longitudinal  separation  which  is  specified  by  Configuration 
Number  2  are  given  in  Figures  75  through  78  for  M  =  1.414  (Subsonic  edgd)  and 
M  =  3.162  (Supersonic  edge)  at  two  reduced  frequencies  equal  to  0.75  and  1.5 
based  on  wing  semi-span  for  wing  bending,  wing  torsion,  tail  bending  and  tail 
torsion  modes  (symmetric).  The  generalized  unsteady  aerodynamic  coefficients 

at  k  “0.75  and  1.5  for  these  modes  are  given  in  TablelS  for  M  "  1.414  and 
s 

in  Table  19  for  M  =  3.162. 

The  pressure  distributions  obtained  for  the  wing-horizontal  tail  .onf iguration 
with  a  combined  moderate  longitudinal  separation,  moderate  vertical  separation 
and  tail  dihedral  which  is  specified  by  Configuration  3  are  given  in  Figures 
79  and  80  for  M  =  1.56.  The  modes  associated  with  this  case  are  symmetric 
wing  bending  and  torsion  and  symmetric  tail  bending  and  torsion  modes  as 
listed  in  Table  16  for  two  reduced  frequencies  equal  to  0.75  and  1.5  based 
on  wing  semi-span.  The  generalized  unsteady  aerodynamic  coefficients  in 
this  case  are  given  in  Table  20. 
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TABLE  17  GENERALIZED  UNSTEADY  AERODYNAMIC  COEFFICIENTS 
FOR  AGARD  CASE  1-1  (M  =  1.56) 


Note  *  denotes  a  value  in  which  the  smoothing  technique  is  applied. 
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TABLE  18  GENERALIZED  UNSTEADY  AERODYNAMIC  COEFFICIENTS 
FOR  AGARD  CASE  2-1  (M  =  1.414) 


REPUCED  FREQUENCY  kc  «0.75 

w 

\ 

WING 

TAIL 

BENDING 

TORSION 

BENDING 

TORSION 

REAL 

1MAG. 

REAL 

IMAG 

REAL 

IMAG. 

REAL 

IMAG. 

WING 

BEND 

-.061454 

-.060726 

.251265 

.252038 

.206619 

,210765 

1.25459 

1.25701 

.000342 

.000286 

.000448 

.000388 

.001930 

.001622 

-.001304 

-.001217 

TORS. 

-.242002 

-.238927 

.861716 

.866631 

.604679 

.625681 

4.62587 

4.64090 

.001303 

.001078 

.001715 

.001471 

.007390 

.006149 

-.004928 

-.004558 

TAIL 

BEND 

-.067101 

-.067969 

-.169840 

-.169749 

-.547404 

-.551913 

-.52169 

-.51881 

.004768 

.005781 

.226054 

.f>25533 

.360195 

.360229 

.369973 

..366786 

TORS. 

-.110388 

-.112006 

-.234990 

-.236638 

-.842003 

-.852978 

-.65508 

-.65769 

.007463 

.009059 

.  1,45316 
.3^6430 

.565537 

.568114 

.615344 

.611553 

}  REDUCED  FREQUENCY  k$  = 

1.5 

\ 

WING 

_ , _ 

“  H 

TAIL 

BENDING 

TORSION 

BENDING 

TORSION 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

WING 

BEND 

-.249145 

-.245253 

.289408 

.289800 

-.318119 

-.304208 

1.38183 

1.38083 

.000905 

.000768 

.000002 

.000021 

.001277 

.000993 

-.001442 

-.001168 

TORS. 

-.961342 

-.946629 

1.019511 

1.022599 

-1.45263 

-1.39449 

5.15053 

5.15137 

.003461 

.002908 

.000020  ' 
.000088 

.004947 

.003818 

-.005501 

-.004410 

TAIL 

BEND 

-.211750 

-.214975 

-.102803 

-.101748 

-1.03469 

-1.04998 

-.27909 

-.27055 

-.012036 

-.003297 

.214599 

.213108 

.342121 

.343239 

.371472 

.367647 

TORS. 

-.319626 

-.325863 

-.116276 

-.115784 

-1.49341 

-1.52411 

-.25048 

-.24191 

-,029u08 

-.022935 

.325291 

.324948 

- . 

,524819 

.529459 

.620377 

.615706 

The  lower  values  are  ones  in  which  the  smoothing  technique  is  applied. 


I 
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TABLE  19  GENERALIZED  UNSTEADY  AERODYNAMIC  COEFFICIENTS 
FOR  AGARD  CASE  2-2  (M  =  3.162) 


REDUCED  FREQUENCY  k$  = 

0.75 

\ 

WING 

TAIL 

\ 

\ 

x 

BENDING 

TORSION 

BENDING 

TORSION 

\ 

REAL 

IMAG. 

REAL 

— 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

BEND. 

-.018634 

.212557 

.278447 

.841141 

0.00 

O.OC 

0.00 

0.00 

UTNG 

-.017578 

.211982 

.279757 

.837123 

0.00 

0.00 

0.00 

0.00* 

TORS. 

-.082534 

.728195 

.924261 

3.082659 

0.00 

0.00 

0.00 

0.00 

-.078735 

.728178 

.931307 

3.072145 

0.00 

0.00 

0.00 

0.00 

BEND. 

-.041406 

-.070962 

-.233794 

-.162768 

.001027 

.125955 

.176227 

.197210 

TAIL 

-.041299 

-.071173 

-.233216 

-.164091 

.000954 

.126363 

.176360 

.197609 

TORS. 

-.061729 

-.098218 

-.349047 

-.214711 

.001443 

.197056 

.288144 

.327728 

-.061322 

-.098806 

-.347626 

-.218341 

.001168 

.197903 

.288335 

.328579 

REDUCED  FREQUENCY  k$  = 

1.5 

\ 

WING 

TAIL 

\ 

BENDING 

TROSION 

BENDING 

TORSION 

\ 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

BEND. 

-.062751 

.228501 

.175233 

.866024 

0.00 

0.00 

0.00 

0.00 

1.1  TM£ 

-.059223 

.227037 

.181251 

.861009 

0.00 

0.00 

0.00] 

0.00* 

WlNb 

TORS. 

-.271330 

.803960 

.490654 

3.212103 

0.00 

0.00 

0.00 

0.00 

-.258011 

.800790 

.515263 

3.198675 

0.00 

0.00 

a.oo 

0.00 

BEND. 

-.106197 

-.020038 

-.404513 

-.030997 

.003007 

.124253 

.176189 

.195748 

TAIL 

-.106506 

-.020820 

-.405956 

-.034804 

.003091 

.125091 

.176565 

.196603 

TORS. 

-.154867 

-.019555 

-.592626 

-.005787 

.003838 

.194329 

.286812 

.325569 

-.154701 

-.021477 

’-.593312 

-.014707 

.003516 

.396067 

.287027 

.  $77  166 

*  The  lower  values  are  ones  in  which  the  smoothing  technique  is  applied. 
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TABLE  20  GENERALIZED  UNSTEADY  AERODYNAMIC  COEFFICIENTS 
FOR  AGARD  CASE  2-3  (  M=1 .56  ) 


REDUCED  FREQUENCY  k§  =  0.75 


\ 

WING 

TAIL 

\ 

\ 

BENDING 

TORSION 

BENDING 

TORSION 

\ 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

BEND. 

-.056069 

.244731 

.225831 

1.194187 

0.00 

0.00 

0.00 

0.00 

ING 

-.054522 

.243512 

.229380 

1.187548 

0.00 

0.00 

0.00 

0.00* 

TORS. 

-.218139 

.837243 

.697713 

4.392965 

0.00 

0.00 

0.00 

0.00 

-.212479 

.836331 

.716314 

4.378940 

0.00 

0.00 

0.00 

0  00 

BEND. 

-.005720 

-.026525 

-.078585 

0.004701 

.218944 

.344831 

.359459 

ML 

-.006964 

-.026167 

SB 

-.068292 

0.007039 

.226998 

.353614 

.350774 

TORS. 

-.000592 

-.045929 

-.091397 

.004576 

.331840 

.529935 

.602395 

-.003525 

-.045710 

.008575 

. 348615 f 

.549506 

.594668 

REDUCED  FREQUENCY  k$  =  1.5 


\ 

WING 

TAIL 

_ _ , _ 

\ 

BEN 

DING 

TORSION 

BENDING 

TROSION 

\ 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

REAL 

IMAG. 

BEND. 

-.233814 

.284202 

-.24294 

1.31985 

0.00 

0.00 

0.00 

0.00 

ING 

-.228539 

.281925 

-.23006 

1.30945 

0.00 

0.00 

0.00 

0/-' 

TORS. 

-.902674 

.996337 

-1.1513 

4.90717 

0.00 

0.00 

0.00 

0.00 

-.882784 

.991482 

-1.0965 

4,87823 

0.00 

0.00 

0.00 

0.00 

BEND. 

-.057175 

-.00831 

-.24617 

.02726 

.000374 

.20908 

.342106 

.  553397 

ATI 

-.054135 

-.00592 

-.22307 

.03691 

.009689 

.21507 

.352259 

.347620 

Ml  L 

TORS. 

-.096174 

-.02787 

-.43834 

-.02555 

-.012916 

.31745 

.514428 

.394510 

-.091677 

-.02285 

-.39995 

-.00283 

.002798 

.33083 

.539510 

.586211 

The  lower  values  are  ones  in  which  the  smoothing  technique  is  applied. 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  RECOMMENDATIONS  FOR  FUTURE  PROGRAM  DEVELOPMENTS 

The  computer  program  developed  for  evaluating  the  unsteady  loadings  on 
wing-tail  combinations  in  supersonic  flow  appears  to  provide  reasonable 
results  when  applied  to  analysis  configurations  that  meet  the  restric¬ 
tions  of  the  "Mach-Box"  method. 

Analysis  configurations  that  lie  within  the  program  limitations  are 
shown  in  figure  25  through  figure  27.  The  analysis  configurations  that 
do  not  meet  the  program  limitations  are  shown  in  figure  28  and  figure  30. 

Aerodynamic  influence  effects  of  the  tail  on  the  wing  is  allowed  only 
for  coplanar  configurations  and  is  not  allowed  for  vertical  separation 
cases.  It  is  recommended  that  the  program  be  extended  to  include  this 
capability  which  would  involve  a  reorganization  of  the  procedures  used 
to  evaluate  the  normal  washes  and  velocity  potentials  of  the  individual 
surfaces.  Modification  of  the  program  would  allow  analysis  of  a  wide 
variety  of  configurations  such  as  "T"  tails,  horizontal-vertical  tails, 
and  bi-plane  configurations  that  would  be  applicable  to  the  space  shuttle 
configurations  now  under  study. 

Unreasonable  load  distributions  may  result  within  an  analysis  on  extremely 
low  Mach  number  for  those  configurations  that  have  large  leading  edge 
sweep  angles  such  as  shown  in  figure  8l.  The  chordwise  loadings  may  have 
large  fluctuations  that  are  caused  by  the  basic  box  pattern  that  defines 
the  leading  edge  and  is  also  affected  by  the  large  normalwash  values  de¬ 
fined  in  the  "off-wing"  diaphragm  region.  Mach  lines  proceeding  from 
the  "outside"  corners  of  the  planform  box  pattern  will  cause  box  loads  to 
be  overestimated  and  the  loadings  on  those  boxes  affected  by  "inside" 
corners  will  be  underestimated.  The  fluctuations  in  the  loadings  may  be 
reduced  by  using  the  subdivision  options  within  the  present  program. 
However,  the  computer  time  requirements  are  severely  affected  by  the  pre¬ 
sent  subdivision  process,  especially  when  the  program  limits  are  being 
used.  The  present  program  subdivides  all  boxes  within  a  specified  stream- 
wise  band  width  from  the  left  hand  Mach  line  to  the  right  ’  \nd  Mach  line. 
This  procedure  does  improve  the  chordwise  load  distributions,  however  it 
is  costly  to  use  on  highly  swept  configurations  at  very  low  Mach  numbers. 
It  is  therefore  recommended  that  the  subdivision  procedure  be  reorganized 
such  that  the  boxes  to  be  subdivided  are  only  those  boxes  in  the  diaphragm 
that  are  cut  by  the  planform  leading-trailing  edges  shown  in  figure  82. 

One  further  item  that  needs  attention  is  the  basic  solution  method  itself. 
The  present  solution  is  modal  dependent  in  that  the  normalwashes  in  the 
diaphragm  and  the  velocity  potentials  on  the  planform  are  evaluated  for  a 
given  mode  shape,  reduced  frequency,  and  Mach  number.  The  process  is 
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repeated  over  again  for  every  mode  shape  to  be  included  in  the  analyses. 
Recent  advances  in  computer  technology  now  make  ib  practical  to  remove 
this  modal  dependency  by  performing  a  matrix  manipulation  on  the  known 
aerodynamic  influence  coefficients  and  known  diaphragm  velocity  po¬ 
tentials  as  shown  in  figure  83.  Application  of  this  procedure  will  re¬ 
move  the  modal  dependency  from  the  solution  and  will  provide  a  program 
that  will  be  an  order  of  magnitude  more  economical  than  presently 
available  in  evaluating  interaction  loadings  in  supersonic  flow. 

RECOMMENDATIONS  FOR  FUTURE  FLUTTER  MODEL  TESTS 


Flutter  analyses  of  Cornell  configuration  15  and  configuration  22  were 
extended  to  evaluate  the  flutter  speed  variations  at  Mach  numbers 
higher  than  those  tested  within  the  Cornell  flutter  model  tests  of  re¬ 
ference  2.  Although  no  tables  of  density  versus  Mach  number  are  avail¬ 
able  for  the  analysis,  an  extrapolated  density-Mach  number  variation 


was  devised  for  use  in  the  analysis  to  represent  the  possible  upper  and 
lower  limits  of  density  as  a  function  of  Mach  number.  The  assumed  den¬ 
sity  variations  are  shown  in  figure  84  and  figure  85  where  the  results 


of  the  flutter  analysis  given  in  terms  of 
are  presented  in  table  21. 


V  , 

W9' 


Numerical  results 


It  should  be  noted  that  there  is  no  large  change  in  the  velocity  parameter 

or  in  the  frequency  ratio  m _  for  configuration  15,  indicating  that 

UJJP 

adverse  coupling  persists  throughout  the  Mach  range  for  the  particular 
value  of  wing  bending  to  fuselage  torsion  ratio  used  in  the  test.  There 
is  a  small  change  in  the  velocity  parameter  with  Mach  number  for  con¬ 
figuration  22,  however,  adverse  coupling  still  remains  for  most  of  the 
Mach  range  investigated.  The  change  in  V/b^W^/PTis  only  slight  when 

compared  with  the  results  of  a  single  surface  analysis  which  displays  a 
dramatic  increase  in  V  _  with  Mach  number. 

The  increased  values  for  single  surface  tests  are  usually  observed  to 
start  within  the  transonic  range  and  are  attributed  to  an  aft  shift  in 
chordwise  center  of  pressure  that  are  accompanied  by  a  reduction  in  CL 
values  with  increasing  Mach  number.  The  source  of  the  adverse  a 

coupling  problem  associated  with  wing-tail  configurations  in  close 
proximity  appears  to  be  identified  with  the  fuselage  flexibility  that 
allows  differential  bending  (or  flapping)  to  take  place  that  causes  aero¬ 
dynamic  driving  force  to  exist  on  the  tail  that  is  out  of  phase  with 
the  wing  motion.  Within  the  differential  flapping  mode,  the  tail  will 
be  aerodynamical ly  driven  in  the  opposite  direction  to  the  wing  motion 
(caused  by  the  normalwash  of  the  wing  wake) .  The  highly  loaded  tail 
will  then  tend  to  mechanically  drive  the  wing  to  larger  amplitudes  (than 
would  be  experienced  by  a  single  surface  alone)  by  means  of  the  fuselage 
torsion  spring.  The  amount  of  circular  coupling  that  is  present  in 
wing-tail  configurations  appears  to  be  a  direct  function  of  relative 
stiffnesses  of  the  +wo  surfaces,  the  span  length  ratios  of  the  surfaces, 
and  the  fuselage  torsional  stiffness  parameter.  Mach  number  effects 
appear  to  have  a  negligible  influence  in  changing  the  velocity  perimeter 
V 


It  appears  that  the  most  critically  coupled  wing-tail  combination  is 
one  that  has  (l)  the  wing  and  tail  in  the  same  plane  (coplanar), (2) 
leading  surface  more  flexible  than  t^e  trailing  surface,  (3)  planform 
taper  ratios  of  1.0  or  larger,  (U)  span  lengths  are  equal,  (5)  and 
fuselage  torsion  flexibility  is  allowed. 

The  maximum  aerodynamic  loading  induced  on  the  tail  (tending  to  aero- 
dynamically  drive  the  tail)  is  obtained  when  the  tail  is  fully  exposed 
to  the  wing  wake.  It  is  expected  that  tall  induced  loadings  will  be 
smaller  for  tail  spans  that  are  less  than  or  greater  than  the  span  of 
the  leading  surface  and  this  in  turn  will  proportionally  reduce  the 
coupling  effects  of  the  system. 

Therefore,  it  is  recommended  that  series  of  supersonic  flutter  model 
tests  be  conducted  on  variable  sweep  wing-tail  models  to  evaluate  the 
effects  on  flutter  speed  due  to  the  critical  parameters  of  (l)  wing-tail 
flexibility  ratio,  (2)  ratio  of  wing-tail  span  lengths,  (3)  fuselage 
torsional  flexibility,  (!<)  and  Mach  number. 

It  is  suggested  that  the  variable  sweep  capabilities  be  incorporated 
only  on  the  wing  (or  .leading  surface).  Planforms  of  the  wing  and  tail 
should  be  representative  of  a  typical  aircraft  configuration  (such  as 
the  F-lll)  to  provide  a  realistic  evaluation  of  the  tail  induced  load 
coupling  effects  of  modern  aircraft.  Inexpensive  model  construction 
techniques  may  be  applied  to  the  model  design,  however,  a  very  accurate 
description  of  the  modal  data,  inertia  and  mass  data,  and  structural 
influence  coefficients  are  required  and  must  be  carefully  defined  and 
documented.  Vibration  tests  and  theoretical  vibration  analyses  must  be 
conducted  to  correlate  the  quality  of  the  mode  shape  definition  prior 
to  tunnel  testing  and  flutter  analyses. 

The  wing  models  to  be  constructed  should  provide  a  large  variation  in 
wing  flexibility  to  evaluate  the  effect  of  wing-tail  flexibility  ratios 
on  flutter  speed. 

The  model  support  system  should  allow  flutter  testing  .  wing-tail  com¬ 
binations  having  horizontal  and  vertical  separations  as  well  as  being 
able  to  test  coplanar  configurations.  The  support  system  should  be  de¬ 
signed  to  allow  testing  over  a  wide  range  of  fuselage  torsion  spring 
rates  to  determine  spring  rates  that  may  preclude  any  adverse  coupling 
effects. 

Finally,  the  Mach  number  range  to  be  covered  within  the  tests  should 
extend  from  very  low  supersonic  Mach  numbers  to  Mach  numbers  that  are 
slightly  higher  than  the  sonic  edge  Mach  number  of  the  planform  having 
the  largest  sweep  angle. 
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FIGURE  25  PLANFORM  EDGE  DEFINITION,  WING  OR  TAIL 
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Wing,  tali.  Arbitrary  dihedral 
and  vertical  reparation. 


FIGURE  26  ALLOWABLE  WING/TAIL  CONFIGURATIONS 


Overlapped  wing  and  tall  allowed  If  mutual  interaction 
does  not  occur.  (  Tall  Mach  cone  does  not  Intersect  wing  ). 


FIGURE  29  TAIL  OVERLAP  RESTRICTION 
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Dihedral  *45  * 
except  vertical  tall 

y  Multiple  dihedral  anpl 

FIGURE  30 

OTHER  CONFIGURATIONS  MOT  ALLOWED 
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(  Alto  tee  figure  33) 
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FIGURE  34-  EFFECT  OF  SUBDIVISION  ON  CHORD  WISE  LOADING. 


CHORDWISE  MACH  BOX  NO. 


FIGURE  35  EFFECT  OF  SUBDIVISION  ON  VELOCITY  POTENTIALS. 

(Reproduced  from  reference  10) 
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Figure  36  VELOCITY  POTENTIAL  SMOOTHING  PROCESS 


FIGURE  37  NON  -  DIMENSIONAL  LIFT  DISTRIBUTION 
(  Reproduced  from  reference  14  ) 


FIGURE  33  ANALYSIS  PLANFORM  OF  CHORDWISE  LOADINGS  -  GAP  STUDY 


SPAN  STATION 


y/» 


0.85 


-0.2  • 

-0.4  - 

-0.6  - 

-0.8  ' 

-1.0  • 

-1.2  • 

-1.4  • 


g 

g 

e 


=  0 
=0.166 


=0.5,  1,0 


FIGURE  39  CHORDWISE  LOADINGS  FOR  VARIOUS  WING  -  TAIL  SEPARATIONS 
IN  WING  TIP  REGION. 


RESULTING  CHORDWISE 


FIGURE  43  WING  DEFLECTION  SHAPE,  MODE  NO.  3,  CORNELL  CONFIGURATION  15 


FIGURE  48  4«?SWEPT  MODEL  ASSEMBLY.  (Reproduced  from  reference  2) 


Cornell  experiment 


Cornell  calculated  nodes 
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Figure  53.  Flutter  Speed  and  Frequency  vs  Horizontal  Separation 

For  Configuration  15,  45°  Swept  Model,  Mach  Number  =  1.238,  z/b„  =  0 
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Figure  54.  Flutter  Speed  and  Frequency  vs  Vertical  Separation  for  Configuration  15, 
45° Swept  Model,  Mach  Number  =  1.238,  x/b.  =  0.505 
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MODIFIED  MODES 
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Figure  55.  Flutter  Speed  and  Frequency  vs  Mach  Number  for  Configuration  15, 
45°  Swept  Model,  z/bs  =  0,  x/b$  =  0.166 
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Figure  56.  Fluttor  Speed  and  Frequency  vs  Mach  Number  for  Configuration  22, 
60°  Swept  Model,  z/bs  =  0,  x/bs  =  0.19 
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Figure  57  Flutter  Speed  and  Frequency  v*  Mach  Number  for  Configuration  23, 
60°  Swept  Model,  z/b,  «  0,  x/b,  =  0 .19 
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FIGURE  59  FLUTTER  SPEED  VS  MACH  NUMBER,  60°  SWEPT  MODELS 

(Reproduced  from  reference  15) 
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Figure  61  .  Longitudinal  Wash  on  A  Sampling  Cho'd,  y/«  ■  0.72,  zjt  ■  0.2 
For  k,  -  0.76  and  M  -  1.66  (Antisymmetric) 


Figure  62.  Longitudinal  Wash  on  A  Sampling  Chord,  y/s  ■  0.72,  ill  *  0.2 
For  k.  *1.5  and  M  -  1.56  (Antisymmetric) 
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Figure  65  .  Upwaah  on  A  Sampling  Chord,  y/»  ■  0.066,  i/$  «  0.5 
For  k,  *  0.76  end  M  *  1.66  (Anttaymmatric) 


Figure  66  .  Upwreh  on  A  Sampling  Chord,  y/s  =  0.056,  zh  -  0.5 
For  k$  *  1.5  and  M  *  1.56  (Antiiymmatric) 
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BOX  PRESSURE  DIFFERENCE  BOX  PRESSURE  DIFFERENCE 


Modi  1  k,=  0.75 


chordwise  oistance 


Figure  6  7 .  Pressure  Distribution  Along  Chord  No.  1,  y/s  »  0.066  For  k,  -  0.75 
M  -  1.56  and  Mode  1  (Antisymmetric  Bending  Mode) 
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Figure  6  8 ,  Pressure  Distribution  Aong  Chord  No.  1,  y/s  •  0.066  For  k,  ■  0.75, 
M  *  1.56  and  Mode  2  (Antisymmetric  Torsion  Mode) 
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Figure  71.  Pressure  Distribution  Along  Chord  IMo.  7,  y/s  =  0.72  For  ks  =  0.75 
VI  =  1.56  and  Mode  1  (Antisymmetric  Bending  Mode) 


Figure  72.  Pressure  Distribution  Along  Chord  No.  7,  y/s  =  0.72  For  k$  =  O.i 
M  =  1.56  and  Mode  2  (Antisymmetric  Torsion  Mode) 
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PRESSURE  DIFFERENCE 


Figure  73.  Premure  Dirtribution  Along  Chord  No.  7,  y/»  -  0.72  For  ks  =  1.5  M  *  1.56 
and  Mod*  1  (Antiiymmatric  Banding  Moda) 
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FIGURE  75  CHORDWISE  PRESSURE  DISTRIBUTION,  AGARD  PLANFORM, 
SMALL  SEPARATION,  WING  PITCH  MODE.  M  «  1.414,  k  -  0.75 
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FIGURE  76  CHORDWISE  PRESSURE  DISTRIBUTION,  AGARD  PLANFORM, 
SMALL  SEPARATION,  WING  PITCH  MODE,  M  =  1.414,  k  =  1.5 


No  smoothing 
Smooth!  ng 


SMALL  SEPARATION,WING  PITCH  MODE,  M  -  3.162,  k,  -  0.75 


MODERATE  SEPARATION, WING  PITCH  MODE,  M  =  1.56,  k. 


FIGURE  81  SOURCE  OF  LOAD  IRREGULARITIES  DUE  TO  L.E.  DEFINITION. 


FIGURE  8  2  RECOMMENDED  REVISED  SUBDIVISION  PATTERN. 
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CORNELL  CONFIGURATION  NO.  15 
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CORNELL  CONFIGURATION  NO.  22 


TABLE  21  ANALYTICAL  FLUTTER  RESULTS  v.s.  FLOW  MACH  NUMBER 


CONF. 

NO. 

1 

t 

M 

P 

slug/ft^ 

(xlO4) 

- ! 

V 

vf 

(ft/sec) 

«f 

wm 

0) 

i>sU0Vy 

U) . 

0 

15 

1.121 

4.0* 

44.5 

858.15 

24.6 

.6192 

.7455 

3.8 

46.84 

869.81 

24.4 

.6119 

.7394 

1.238 

3.5* 

50.9 

1006.78 

26.51 

.6794 

.8035 

3.2 

55.63 

, 

1039.55 

26.39 

.6710 

.7997 

1.350 

3.2 

55.63 

993.77 

26.46 

.6417 

.8018 

2.7 

65.93 

1069.95 

26.43 

.6348 

.8008 

1.414 

3.0 

59.33 

1000.87 

26.17 

.6255 

.7929 

2.5 

71.20 

1083.46 

26.16 

.6180 

.7928 

1.60 

2.6 

68.44 

1093.25 

25.60 

.6363 

.7758 

2.0 

89.00 

1237.66 

25.58 

.6318 

.7751  ' 

1.80 

2.15 

82.79 

1217.05 

25.05 

.6443 

.7592 

1.50 

118.67 

1460.09 

25.42 

.6456 

!  .7702 

2.40 

1.0 

178.00 

1889.93 

24.28 

.6819 

.7358 

0.5 

356.00 

2702.31 

24.50 

.68% 

.7424 

22 

1.151 

3.3 

53.43 

916.59 

21.67 

.6060 

.8773 

1.240 

2.7* 

65.30 

1139.70 

23.70 

.6815 

.9594 

1.414 

1.85 

95.27 

1180.40 

22.13 

.5844 

.8959 

1.70 

103.70 

1224.82 

22.09 

.5813 

.8944 

1.750 

1.1 

160.31 

2004.95 

25.06 

.7653 

1 .0147 

0.7 

251.85 

2465.35 

24.87 

.7506 

1.007 

2.100 

0.72 

244.88 

2431.23 

24.38 

..7508 

.9870 

0.30 

587.70 

3742.02 

24.28 

.7460 

.9829 

2.400 

0.50 

351.62 

i_  . 

2786.76 

23.22 

.7173 

.9401 

*  denotes  the  test  values. 
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